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Currently,  there  is  a concerted  effort  being  made  by  the  Navy 
and  other  Governmental  agencies  in  exploring  the  feasibility 
of  using  alternate  concepts  to  present  day  naval  ship  design 
for  the  Navy  of  the  future.  These  investigations  have  led  to 
the  consideration  of  Air  Cushion  Vehicles  (ACV)  and  Surface 
Effects  Ships  (SES)  as  viable  candidates.  These  vehicles  offer 
the  potential  for  much  greater  versatility  and  higher  operational 
speeds  than  hithertofore  possible  with  conventional  ship  design. 
The  ACV  with  its  totally  flexible  skirt  system  presents  an 
amphibious  capability  most  attractive  for  coastal  and  near- 
shore operations,  assault  landing  operations,  and  for  arctic 
environmental  use.  The  SES,  on  the  otherhand,  while  not  of  an 
amphibious  nature,  provides  an  ocean  going  vehicle  capable  of 
very  high  speed  performance  in  reasonable  sea  states  and  weather 
conditions . 

Interest  in  these  concepts  has  led  the  Navy  into  a development 
program  in  which  two  air  cushion  assault  vehicles  are  presently 
being  constructed  for  evaluation  purposes.  Additionally,  two 
100-ton  surface  effect  ships  have  been  built  and  tested,  under 
Navy  contract,  with  sufficiently  encouraging  results  that  the 
Navy  is  currently  conducting  a detail  design  of  a 3000-ton 
class  SES.  It  is  apparent  from  this  activity  that  more  than 
just  casual  interest  is  being  given  to  these  vehicles  and  in- 
deed, dependent  on  the  results  of  the  above  programs,  they 
may  prove  to  be  the  forerunners  of  a completely  new  class  of 
fighting  ship  for  the  Navy  of  tomorrow. 

The  advent  of  the  Surface  Effect  Vehicle  as  a serious  contender 
for  Naval  applications  has  led  to  the  need  for  an  evaluation  of 
the  vulnerability  of  this  type  of  craft  under  typical  tactical 
situations.  As  presently  envisaged  the  role  these  vehicles  are 
to  play  in  naval  operations  is  one  of  antisubmarine  warfare  (ASW) , 
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escort  duties  and  near  or  offshore  patrol  and  rescue,  which 
operations  require  a dash  or  high  speed  capability  coupled  with 
maneuverability,  a feature  characteristic  of  air  cushion  vehicles 
(ACV)  and  surface  effect  ships  (SES)  alike. 

Due  to  this  mounting  interest  it  is  appropriate,  at  this  time, 
to  obtain  an  assessment  of  the  vulnerability  of  such  craft  to 
possible  threats.  In  identifying  possible  threat  areas  one 
outstanding  possibility  is  that  due  to  explosion  generated  waves. 
Past  experience  in  this  field.  Reference  1 and  2,  has  shown  the 
great  damage  potential  such  a phenomenon  can  have  on  submarines 
and  conventional  ships.  The  effects  on  ACV's  and  SES  are  ex- 
pected to  be  of  greater  significance  since  the  unique  features 
of  these  vehicles  make  them  particularly  susceptible  to  sudden 
and  anomalous  changes  in  sea  surface  topography,  such  as  are 
known  to  be  produced  by  nuclear  detonations. 

Past  studies  have  been  primarily  concerned  with  the  behavior 
of  ships  and  submarines  within  the  transient  surf  zone  produced 
by  high  yield  explosions  at  the  continental  margins  (Van  Dorn 
Effect) . However,  because  of  their  dynamic  response  we  expect 
that  the  damage  potential  on  SES  and  ACV's  cannot  only  be 
restricted  to  these  conditions  but  must  be  extended  to  include 
the  effects  of  small  and  moderate  yield  devices  and  operations 
in  deep  water.  It  is  evident  that  even  under  these  latter 
conditions  waves  can  be  produced  that  are  capable  of  limiting 
the  performance  of  these  craft. 


The  radical  differences  between  the  design  of  these  craft  and 
those  of  present  naval  ships  makes  it  impossible  to  extrapolate 
the  results  obtained  in  past  studies  to  the  present  case.  It 
is  only  by  conducting  an  investigation,  wherein  the  features 
f of  these  vehicles  are  faithfully  modeled,  that  the  vulnerability 

of  these  craft  can  be  determined. 
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In  light  of  the  above  discussion,  it  is  deemed  imperative  that 
such  a study  be  conducted  with  the  objective  of  defining  the 
operational  limits  of  ACV's  and  SES  under  explosion  generated 
waves,  and  to  ascertain,  where  possible,  the  survival  potential 
of  these  vehicles  when  subjected  to  tactical  situations  of  this 
nature . 

The  criteria  used  in  defining  the  structural  design  and  stability 
characteristics  of  SES  are  derived  principally  from  the  desired 
operational  envelopes.  The  envelope  defines  the  speed-wave 
height  domain  over  which  the  craft  will  operate.  Typically, 
such  an  envelope  is  shown  in  Figure  1. 

Two  factors  which  greatly  affect  the  basic  structural  design 
of  the  SES  are  the  highest  wave  environment  to  be  encountered 
when  operating  on-cushion  and  the  maximum  impact  loading  to  be 
seen  by  the  hull  during  operation.  The  former  factor  is  of 
prime  importance  in  selecting  the  height  of  the  flexible  skirt 
system  and  thus  impacts  hull  design.  The  latter  determines  plating 
thickness  and  consequently  weight.  From  Figure  1 it  will  be  seen 
that  point  A on  the  chart  determines  maximum  wave  height  on 
cushion.  The  worst  combination  of  sea  state  and  speed  will  be 
determined  by  line  AB  along  which  maximum  impact  loads  are  likely 
to  occur.  If  such  an  operating  envelope  is  determined  without 
due  consideration  for  potential  threats  as  outlined  above  grave 
consequences  can  arise.  It  is  easily  conceivable  that  a wave 
environment  outside  the  typical  boundaries  now  being  considered 
in  the  SES  field  can  be  generated  by  low  to  moderate  yield 
devices.  Such  circumstances  could  cause  structural  and  operational 
failures . 

In  addition  to  the  above  impacts  the  question  of  craft  stability 
and  survival  are  of  equal  importance.  The  response  of  an  SES  to 
a typical  explosion  generated  wave  profile  could  lead  to  con- 
ditions of  craft  plow-in,  pitch  poling  and  capsizing.  Such 
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extreme  motions  are  indeed  possible  under  certain  conditions  of 
speed,  water  depth  and  yield.  This  aspect  of  vulnerability  is 
therefore  of  equal  importance  in  analyzing  SES  operational 
characteristics . 

The  problem  at  hand  can  be  divided  into  two  basic  sub  tasks: 

(a)  The  analytical  description  and  modeling  of  explosion 
generated  water  waves , and 

(b)  The  analytical  treatment  of  the  craft  dynamics  and 
motions  when  subjected  to  a disturbing  functions  as 
defined  in  (a)  above. 

Whereas  previously  conducted  work  by  Tetra  Tech,  References  1 
and  2,  is  directly  applicable  to  the  first  of  these  areas,  the 
second  provides  a new  and  added  dimension  due  to  the  radical 
difference  between  ACV  and  SES  and  conventional  ships.  Analytical 
modeling  of  SES  motions  and  maneuvering  however,  have  also  been 
conducted  by  Tetra  Tech,  Reference  3 and  has  been  used  as  a 
basis  of  departure  for  the  present  program. 

The  present  report  deals  with  the  first  phase  in  the  investigation 
of  the  response  of  a typical  SES  to  an  explosion  wave  environment. 
This  initial  phase  has  been  directed  to  the  formulation  and 
development  of  the  analytical  model  describing  the  dynamics  of 
a surface  effect  vehicle  and  the  description  of  the  explosion 
generated  wave  environment.  This  analytic  model  has  been  used 
to  assess  the  effects  of  a chosen  explosion  condition  on  an 
SES.  Exercise  of  the  program  in  this  area  has  been  limited, 
pending  the  start  of  the  second  phase  which  will  investigate 
various  parameters  of  the  problem  such  as  the  effects  of  yield, 
standoff  distance,  vehicle  size,  water  depth  and  tactical 
maneuvers  to  enhance  survival. 
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In  order  to  fully  exercise  the  analytic  program  and  ensure  its 
validity  several  cases  of  sinusoidal  waves  and  solitary  waves 
were  also  run.  These  latter  waves  are  representative  of 
waves  in  the  shallow  water  environment  and  consequently  are 
worthy  of  investigation  in  their  own  right. 


The  work  described  in  this  report  was  conducted  for  the  Office 
of  Naval  Research  under  contract  N00014-76-C-0261.  This  report, 
as  already  mentioned,  covers  work  during  the  first  phase  of  the 
contract  and  is  consequently  an  interim  report. 


2.  FORMULATION  OF  PROBLEM 
2.1  Coordinate  System 

The  motion  of  the  craft  will  be  described  in  terms  of  the 
relationship  between  a body  fixed  reference  frame  and  a co- 
ordinate system  fixed  in  space.  The  initial  coordinates 
(xQ  ,y0  / z0  ) and  the  body  coordinates  (x,y,z)  are  both  designated 
according  to  a right  hand  convention  with  zQ  and  z positive 
downward  as  shown  in  Figure  2.  The  origin  of  the  body  frame  is 
kept  fixed  at  the  center  of  gravity  of  the  craft  for  all  time,  t. 
The  x -axis  is  parallel  to  the  baseline  of  the  craft,  positive 
forward,  and  positive  y is  therefore  pointing  to  starboard. 

The  two  coordinate  systems  coincide  initially  at  time  zero. 


The  body  moves  with  six  degrees  of  freedom;  at  time  t,  there  will 
be  three  linear  displacements  (surge  x,  sway  y,  and  heave  z) 
and  three  angular  displacements  (roll  4> , pitch  0,  yaw  f)  . 

2.2  Kinematic  Relations 

As  the  body  moves  with  six  degrees  of  freedom,  forces  and 
moments  are  generated  and  act  on  the  body.  For  the  convenience 
of  analysis,  we  shall  resolve  the  total  force  into  three  components 
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about  the  body  axes.  Definitions  and  symbols  of  the  six  com- 
ponents of  force,  displacement  and  velocity  are  given  by  Table  1. 


TABLE  1 


Motion 

Force  or  Movement 

Displacement 

Velocity 

Longitudinal 

X 

X1 

u 

Lateral 

Y 

n 

V 

Normal 

Z 

*i 

w 

Roll 

K 

CD 

p 

Pitch 

M 

9 

q 

Yaw 

N 

♦ 

r 

It  is  noted  that  the  linear  displacements  along  the  inertial 
axes  have  been  defined  by  x,  y,  and  z.  These  are  related  to 
the  displacements  along  the  body  axes  by  the  following  equations: 


x 

y = 


cos  9 cos  y sin  9 sin  j>cosf  - cos^sinf  sin9cos<f>cosy+sin4isiny 
cos  9 sin  y sin  9 sin  $siny  + cosucosy  sin9cos ^siny-sin^cosy 
- sin  9 cos9  sin*  cos  9 cos  $ 


In  Table  1,  u,  v and  w are  the  linear  velocities  along  the  body  axes 
x,  y and  z,  and  p,  q and  r are  the  angular  velocities  about  these 
axes  respectively. 
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2.3  Equations  of  Motion 

Consider  the  craft  in  six  degrees  of  freedom  and  let  u,  v,  and 
w be  the  linear  velocity  components  of  the  craft  center  of  gravity 
along  the  body  axes  x,  y and  z and  p,  q and  r be  the  angular 
velocities  about  these  axes,  respectively.  The  equations  of 
motion  in  this  body  coordinate  system  are  then  given  by 

m (u  + qw  - rv)  = x 
m (v  + ru  - pw)  = y 
m (w  + pv  - qu)  = z 
lx  p + (Iz  - Iy)  qr  =K 
Iy  q + (lx  - Iz)  rp  =M 
Iz  r + (Iy  - lx)  pq  =N 


Where  m is  the  mass  and  lx,  Iy  and  Iz  are  the  moments  of  inertia 
of  the  craft  about  the  respective  axis.  Terms  on  the  lefthand 
side  represent  the  rigid  body  inertial  reactions  and  the  cen- 
trifugal effects  acting  at  the  origin  with  respect  to  the 
moving  coordinate  system.  The  terms  on  the  righthand  side  refer 
to  the  total  forces  and  moments  applied  to  the  craft,  including 
the  hydrodynamic  effects  arising  from  the  overall  motions  of  the 
craft  as  well  as  the  results  of  propulsion  and  control  forces 
which  may  affect  the  craft  maneuvers.  In  a functional  form, 
these  components  can  be  expressed  generally  as: 


z > = f (u, v,w,p,q, r ,u, v,w,p ,q,r ,xQ ,yQ , Zq , $ ,0 , f , 6 , e)  (3) 
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In  the  above  equation,  x , y , and  z are  the  position  components 
of  the  linear  displacements  of  the  craft  and  <f> , 0 , and  ip  are  the 
angular  displacements.  The  parameter  <5  represents  a general 
description  of  the  effect  of  various  propulsion  and  control 
schemes,  and  the  parameter  c represents  the  effect  due  to  the 
environmental  disturbances  such  as  waves.  The  functional  form 
of  the  equation  shows  clearly  the  dependence  of  the  external 
force  and  moment  on  the  various  variables.  To  reduce  the 
above  functional  relationship  into  a useful  mathematical  form, 
a Taylor  expansion  is  usually  applied  provided  that  the  linear 
and  non-dimensional  proportionality  constants,  are  known  or  deter- 
minable. By  keeping  a sufficient  number  of  terms  for  each 
variable,  the  forces  and  moments  can  be  expressed  in  any 
desired  order  of  these  variables  to  account  for  the  non-linear 
effects . 

The  determination  of  the  proportionality  constants,  or  the 
hydrodynamic  derivatives,  by  analytical  methods  is  generally 
limited  only  to  the  linear  terms.  The  non-linear  coefficients 
are  normally  determined  experimentally  by  means  of  captive  model 
tests.  In  the  present  analysis,  however,  external  forces  ana  mo- 
ments are  determined  analytically  on  the  basis  of  physical 
concepts.  While  the  non-linear  coefficients  are  not  explicitly 
identifiable  by  this  approach,  this  method  is  more  convenient 
to  include  various  non-linear  features  without  the  backup  of 
experimental  information.  The  general  representation  of  the 
total  force  (or  moment)  acting  on  an  SES  is  assumed  to  be  com- 
posed of  various  components  as  follows: 


where  i = 1 to  6,  represents  a particular  mode  or  direction  of 
motion.  The  calculation  of  each  of  the  component  forces  is 
discussed  in  the  following  section. 

3.  FORCES  AND  MOMENTS  - CRAFT  DYNAMICS 

3.1  Sidehull  Forces 

The  calculation  of  the  forces  acting  on  the  sidehulls  assumes 
that  these  forces  fall  into  two  major  catagories,  namely  viscous 
and  non-viscous  components.  The  non-viscous  components  of  forces 
and  moments  are  those  directly  related  to  dynamic  fluid  pressure 
resulting  from  the  sidehull  motion.  These  forces  are  intimately 
associated  with  the  energy  exchanges  between  the  fluid  and  the 
moving  sidehull  and  can  be  deduced  from  the  fundamental  principles 
of  classical  mechanics.  Consequently  all  non-viscous  terms,  both 
linear  and  non-linear,  can  be  analytically  identified  as  functions 
of  the  body  added  inertia,  provided  that  the  non-viscous  dissi- 
pative damping  is  negligible.  The  viscous  components  are  drags 
created  through  various  origins.  The  term  drag  customarily 
refers  to  the  total  resistance  of  the  craft  in  its  axial  direction 
which  consists  of  many  components  from  many  different  items,  and 
will  be  considered  in  detail  separately,  in  a later  section.  In 
the  present  section,  only  contributions  due  to  viscous  cross 
flows  on  the  sidehulls  are  considered.  These  contributions  are 
normally  treated  as  dependent  on  the  square  of  the  velocity 
through  proportional  empirical  constants.  Some  details  for  the 
calculation  of  both  the  viscous  and  non-viscous  forces  on  the 
sidehull  are  given  in  the  following: 

(a)  Hydrodynamic  pressure  on  sidewall 

Because  of  the  narrow  thin  geometry,  the  calculation  of  the 
hydrodynamic  forces  on  the  sidehull  can  be  performed  according 
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to  the  fundamental  concept  of  slender  body  theory.  For  a 
slender  body  of  constant  speed  U in  an  inviscid,  incompressible 
fluid  the  linearized  free  surface  condition  is  given  by: 


2 

U i> 


+ 9 *z  = o 


(4) 


Where  $ is  velocity  potential.  Since  the  sidehull  immersion  is 
normally  small  in  comparison  with  the  craft  length,  a normalized 
equation  for  the  above  condition  can  be  written  as  follows: 

2 


d * 

l X ' X 


+ $ 


Z ' 


= 0 


(5) 


Here  F is  the  Froude  number  based  on  craft  speed  and  sidehull 
length.  x'  is  the  non-dimensional  axial  coordinate  referenced 
to  the  craft  length  i,  and  z'  is  the  non-dimensional  vertical 
coordinate  referenced  to  craft  immersion  d.  From  the  above  equation, 
it  is  clear  that  even  though  F may  be  large  (typically  of  the 


order  of  1 or  2 for  an  SES)  the  fact  that  the  immersion  ratio 
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d/j i is  small  (of  the  order  of  10  ) makes  the  second  term  dominate. 


Consequently,  the  free  surface  can  be  regarded  as  a reflection 
boundary  where: 

(6) 


$ , =o 

z ' 


which  is  equivalent  to  the  condition  for  a positive  reflection 
in  the  free  surface. 


The  derivation  of  the  boundary  condition  suggests  that  the  problem 
can  be  treated  as  a body  moving  in  an  infinite  medium,  in  which 
the  dissipative  damping  is  negligible.  Consequently,  as  shown  by 
Lamb,  reference  4,  the  hydrodynamic  effect  on  the  body  is  entirely 
determinable  as  a function  of  its  added  mass  along  the  principal 
axes  of  the  body.  Following  the  procedure  of  classical  mechanics, 
the  effects  of  the  hydrodynamic  pressure  on  the  craft  can  be 
obtained . 
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In  deriving  the  force  relations,  we  shall  break  the  three 
dimensional  sidehull  into  a number  of  segments  along  the  longi- 
tudinal axis.  Each  segment  will  be  considered  individually  as  a 
two-dimensional  problem;  interferences  between  segments  will  be 
ignored.  This  is  the  basic  approach  of  the  slender  body  technique 
which  is  adopted  in  the  present  analysis. 

In  so  doing,  relative  fluid  velocities  at  the  center  of  a segment  Z, 
from  the  normal  plane  are  given  by 


u 

r 

(5/ 1)  = u 

V 

r 

+ 

> 

n 

UJ> 

£r  - 

f 

P 

(7) 

w 

r 

u,t)  = w - 

Cq  + 

b 

P 

where 

P = 

\ 

q = 

0 

r = 

i 

and 

2 

u 

2 2 
+ v + w = 

u2 

where  U is  the  resultant  velocity  of  the  craft  o,i,b,  and  f 
are  the  lateral  and  normal  moment  arms  about  the  center  of 
gravity  respectively.  The  above  relations  are  applicable  to 
both  the  starboard  and  port  sidehulls;  for  the  port  sidehull 
however,  a negative  value  of  b should  be  used. 

We  shall  first  consider  the  segment  as  being  axially  symmetric 
and  having  component  added  masses  m and  mzz  along  the  lateral 
and  normal  axes,  respectively.  For  asymmetrical  segments  with 
respect  to  the  axial  axis,  additional  treatment  will  be  con- 
sidered later. 
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Specifically , m 

and 

m can  be  written  as  follows: 
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= k 

TT 
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yy 
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IT 
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(8) 

n2(C) 
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(O 

= k 

zz 

zz 

where 

k 

yy 

and  k 

zz 

are 

the  added  mass  coefficients  which  are 

generally 

functions 

of  geometry 

and  frequency;  n(£)  is  the  local  beam 

at  the 

water 

line 

and  c 

(C) 

is  the  local  draft. 

The  added  mass  component  along  the  axial  direction  is  ignored 
according  to  the  slender  body  approach. 


The  kinetic  energy  of  a unit  slice  of  the  fluid  can  then  be  written 
as 


T (§.t)  ( 
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yy 
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Neglecting  the  second  order  terms,  the  hydrodynamic  forces  and 
moments  per  unit  axial  length  are  given  by ■ 
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The  kinetic  energy  T at  a fixed  normal  plane  is  a function  of  g 
and  t.  The  total  derivative  therefore  must  reflect  the 
changing  coordinate  g^  of  the  normal  plane  with  time.  Thus 

d _ _9_  + Sen  5 _ u 5 

dt  " at  at  ag  ' st  a§ 

Substituting  (9)  into  (10)  to  (14),  carrying  out  the  differentiation, 
and  then  integrating  over  the  sidewall  length,  gives  the  following 
expressions  for  the  normal  and  lateral  forces  and  moments 
acting  on  the  body.  In  each  of  these  expressions  the  integrals 
are  taken  from  stern  to  bow. 
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The  first  subscript  designates  the  contributing  force  component 
and  the  second  subscript  stands  for  the  motion  direction.  For 
instance,  represents  the  rolling  moment  generated  by  lateral 

forces  which  are  induced  by  lateral  motions  vr  and  vr. 


It  has  been  mentioned  earlier  that  m and  m are  for  axially 

yy  zz  J 

symmetric  sections.  More  often,  however,  the  sidewall  sections 
are  asymmetrical.  These  asymmetries  give  rise  to  cross  coupling 
effects  which  are  estimated  as  follows: 


m 

(£) 

= k . m ( £) 

yz 

y zz 
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(£) 

= k . m ( £) 

zy 

z yy 

m 

yz 

(?) 

represents  the 

relating  the  fluid  momentum  in  the  lateral  direction  y to  the 
local  normal  velocity  in  the  direction  z.  Similarly,  m (£) 
can  be  interpreted  as  the  added  mass  relating  vertical  fluid 
momentum  to  the  local  lateral  velocity.  The  coefficients  k 

y 

and  kz  are  estimated  using: 


k = N (£) 

y y 


where  N (£)  and  Nz  (£)  are  average  values  of  the  horizontal 
and  vertical  unit  normal  components  of  the  hull  cross-section 
at  any  station.  The  average  is  taken  with  respect  to  the  wetted 
length  of  the  hull  cross  sectional  area.  Note  that  for  a side 
hull  with  axial  symmetry  N (£)  is  zero  and  hence  there  would  be 
no  cross  coupling  forces  or  moments. 

The  forces  and  moments  due  to  cross  coupling  are  then  given  by: 
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(b)  Hydrostatic  Forces  and  Moments 

The  hydrostatic  force  acting  on  the  body  is  obtained  by  inte- 
grating the  hydrostatic  pressure  over  the  entire  wetted  body 
surface  and  is  numerically  equal  to  pgA,  where  p is  the  density 
of  the  fluid,  g is  the  acceleration  of  gravity  and  A is  the 
volume  of  the  displaced  fluid.  Let  the  sectional  area  at  station 
£ be  s(£),  which  is  a function  of  local  draft  5 defined  as 

?(C)  = D (£)  + z - £sine  + Bsin<J>  (29) 

where  D (£)  is  the  initial  local  draft  of  the  body  at  station  £ . 

The  total  buoyancy  force  is  then  given  by: 
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FB  = 0&fs  s(b)d*  (30] 

where  the  integration  is  carried  from  stern  to  bow.  The  force 

component  along  the  body  normal  axis  z is  then  given  by 

7 — „ /*b  /_%  


= -p  g cos  o 


?y*bs  (5)  d§ 


and  the  component  along  the  longitudinal  axis  is 
xs  = p g sin  ey*b  s(c  ) 


The  hydrostatic  pitching  moment  is  given  by: 


Ms  = P 8 

Similarly,  the  rolling  moment  is  given  by: 


rb 

= ~ P S Js  3 (?  > £ ) • b {5  ) d§ 


If  b-B  is  small  as  compared  to  B,  then  approximation  of  K 


can  be  simply  obtained  by 


= B . Zg  = - pgB  f %(£)<% 


(c)  Axial  Drag 


The  axial  drag  on  the  sidehulls  arise  from  two  basic  sources. 
Firstly,  the  frictional  drag  caused  by  the  viscous  effects  of 
the  fluid  over  the  body,  and  secondly  the  base  pressure  drag 
which  arises  due  to  the  separated  wake  existing  aft  of  the 
transom. 
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The  sidehull  viscous  drag  is  primarily  a function  of  the  Reynolds 
number  and  surface  finish  of  the  body  and  is  determined  from: 


Drag  = y p u2  Sw 


where  Sw  is  the  wetted  surface  given  by 


and  G (?)  is  the  girder  length  at  section  5. 


The  frictional  drag  coefficient  may  be  obtained  by  using 


the  standard  ITTC  relationship  given  by 


_ 0-075  _ 

f " (Dog10  Rn  -2 V 


where 


Rn  = 


. the  Reynolds  number, 


in  which 


L = sidewall  length 


v = kinematic  viscosity  of  fluid. 


In  addition,  a pressure  drag  component  exists  on  the  sidewall; 
because  of  the  usual  design  of  this  type  of  body  being  long  and 
slender  and  having  a sharp  transom.  This  can  be  estimated  by 
an  additional  drag  coefficient  given  by  the  following  expression: 


where 
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(38) 


= drag  coefficient  based  on  base  area. 


= frictional  coefficient  based  on  base  area. 


This  value  of  drag  holds  provided  the  base  area  is  not  fully 
ventilated.  At  high  speeds  and/or  at  shallow  immersions  the 
likelihood  of  ventilation  is  almost  certain.  When  this  occurs, 
the  base  drag  coefficient  is  given  by 


B 


whete 


H 


fh2 


= Froude  number  based  on  transom  immersion 


(39) 


U . 


The  transition  between  a fully  vented  and  viscous  wake  is 
determined  from  the  following  empirically  established  relation: 

if  >_  3.2  the  base  is  fully  vented. 

The  force  and  moment  components  due  to  axial  drag  are  thus 
given  by 


X = - 
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u2  [Sw  C + S(l)  . C l 
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= - 2 P u2[Cfy'c(5)h(’  )d*  + CB  S(l)f(l)] 


N = - B X 
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(40) 

(41) 

(42) 
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where  h(?)  is  the  vertical  moment  arm  of  the  girder  at  section  5 
and  S ( 1)  is  the  immersed  transom  area. 


In  addition  to  the  normal  components  of  drag,  viz  the  skin 
friction  drag  and  the  pressure  drag,  there  exists  expecially  at 
high  speed,  a significant  spray  drag.  Unfortunately  very  little 
information  exists  regarding  this  drag  component.  However,  using 
the  results  of  some  experimental  work,  references  3 and  5 the 
spray  drag  can  be  expressed  in  a general  form  as: 

Dspray  = f(<3'c't}  (43) 

where  q is  the  dynamic  pressure,  c is  the  characteristic  length 
from  the  point  of  generation  of  the  spray  to  the  maximum  thickness 
point  and  t is  the  maximum  thickness  of  the  body. 


Based  on  the  results  of  reference  3 the  following  formulae  were 
obtained  for  the  spray  drag  caused  by  a typical  SES  sidehull 
conf igeration : 

D = 0.75  C.  qct 

spray  f n (44) 


In  this  formula  the  value  of  t is  taken  to  be  the  maximum 
thickness  in  the  waterline  plane  and  the  friction  coefficient  Cf 
is  evaluated  at  the  appropriate  Reynolds  number. 

This  result  has  shown  excellent  agreement  with  the  test  results. 


(d)  Viscous  cross-flow  effect 


As  mentioned  this  component  of  the  sidehull  forces  and  moments 
is  a non-linear  term  arising  due  to  the  real  fluid  effects 
occurring  on  the  sidehull.  The  contribution  of  this  term  to  the 

overall  force  on  the  sidehull  is  small  for  small  excursions  of 
the  hull  but  becomes  the  dominant  term  as  craft  motions  become 
greater.  This  force  is  usually  cast  in  the  form: 
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Cross-flow  forces  = ^pc^S  [ | v (45) 

where  CD  = cross  flow  drag  coefficient 

S = projected  area  of  the  sidehull 
V = relative  flow  velocity 

The  coefficient  CD  is  a function  of  geometrical  shape  of  the 
body  and  Reynolds  number.  It  is  usually  obtained  from  experi- 
mental data  by  judicial  interpretation  of  the  results  from 
tests  done  on  idealized  geometric  shapes. 

Accordingly,  the  forces  and  moments  due  to  this  effect  are 
given  by: 
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In  addition  to  these  terms  there  are  also  cross  coupling  terms 
arising  due  to  these  cross  flow  drag  forces.  These  forces  are 
given  by: 
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where  the  superscript  C designates  cross  flow,  and  the  subscripts 
have  the  same  meaning  as  defined  previously. 

The  value  of  the  cross  coupling  coefficients  k and  k in 

c 3 yc  zc 

this  case  are  taken  to  be  proportional  to  the  cotangent  of  the 
local  deadrise  angle  of  the  sidehull  at  any  given  section. 


3.2  Cushion  Pressure  Forces 


In  addition  to  the  forces  imparted  to  the  craft  through  the 
sidehulls,  the  cushion  pressure  supporting  the  craft  has  a 
significant  effect  on  the  craft  dynamics.  For  the  present 
investigation,  since  a general  type  of  craft  is  being  considered, 
the  supporting  air  cushion  is  considered  as  basically  a rectangular 
box  bounded  by  the  sidehulls  and  the  forward  and  aft  seals.  The 
plenum  is  fed  by  a fan,  or  system  of  fans,  with  a specified 
fan  characteristic.  Details  of  the  fan  ducting  and  heave 
alleviation  devices  which  are  usually  used  on  such  craft  have 
not  been  included  in  the  analysis  as  this  would  require  a more 
detailed  definition  of  the  lift  fan  system,  a task  beyond  the 
scope  of  the  present  preliminary  study. 

The  basic  equation  governing  the  air  flow  into  and  out  from  the 
cushion  is  the  conservation  of  mass  which  states  that 


ih  = p (Q . - Q 

in  out 


(53) 
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where  m = rate  of  change  of  mass  in  the  plenum 
Q^n  = total  flow  into  the  plenum 

Qout  = leakage  flow  out  under  the  seals  and  sidehulls 

The  flow  into  the  air  plenum  is  governed  by  the  lift  fan 
characteristic  which  is  represented  by: 

Qin  = Qf  = + *1  Pc  + ^2PC  (59) 

where  $ , i $2  are  constants 
p = cushion  pressure 

The  leakage  flow  is  considered  to  be  governed  by  an  orifice 
type  flow  equation  given  by: 


^out  C0  AL  ^ Pc  Pa 

P 


where  CQ  = discharge  coefficient 
p = density 

p = atmospheric  pressure 

Cl 

Pc  = cushion  pressure 
A = leakage  area 

Li 


(60) 


The  leakage  area  in  this  equation  is  comprised  of  several 
components.  These  can  be  represented  as 


AL  = 


+ A + 
sw 


where  Aq  = equilibrium  leakage  flow 

A = leakage  flow  under  the  sidehull 
sw 

As  = leakage  flow  under  the  seals 


(61) 


25 


The  above  equilibrium  leakage  flow  is  that  leakage  required  to 
maintain  the  craft  at  a given  equilibrium  condition  when  not 
disturbed  by  any  waves.  This  leakage  area  can  be  adjusted  by 
changing  the  setting  of  the  seals  under  actual  conditions  and 

it  determines  the  equilibrium  immersion  of  the  craft.  The 
equilibrium  pressure  is  obviously  given  by 

(pc  ‘ pa}  Ac  = W - FB  (62) 

where  W = craft  weight 

FB  = buoyancy  force 

A = plenum  area 
c 

The  areas  Agw  and  Ag  are  obtained  at  each  instant  in  time  by 
integrating  the  clearance  between  the  sidehull  and  seals  with 
respect  to  the  local  water  elevation. 

This  leakage  area  obviously  changes  as  a function  of  time 
depending  on  the  craft  motions  and  the  free  surface  elevation. 

The  pressure  in  the  plenum  is  assumed  to  vary  according  to  an 
adiabatic  compression  law,  viz 

(Pc  + Pa>  v'  = constant  (63) 

where  V = plenum  volume. 

From  the  conservation  of  mass  equation  the  following  equation 
can  be  obtained, 

V - °in  - °out  <««> 

on  substituting,  m =pV. 
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Using  the  above  equations  the  pressure  and  flows  into  the 
plenum  can  be  found  and  the  resulting  forces  and  moments  on  the 
craft  can  be  calculated  as  follows: 
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where  VCG  = vertical  height  of  CG  above  keel 
B = Width  of  the  plenum 

(j>  = roll  angle 

1 = length  of  the  plenum 

In  addition  to  the  above  forces,  the  pressure  acting  on  the  free 

surface  causes  a wave  drag  effect  which  has  to  be  accounted  for 
in  the  computations.  The  wave  drag  has  been  discussed  extensively 
in  the  literature. 


The  wave  drag  referred  to  here  includes  that  due  to  the sidehulls , 
the  pressure  planform,  and  their  interactions.  The  calculations 
for  the  wave  resistance  of  a pressure  patch  and  a pair  of  thin 
walls  is  straight  forward  provided  that  the  geometrical  form  is 
simple.  Following  the  method  of  reference  6,  the  total  wave 
resistance  for  a combination  of  a pressure  planform  and  two 
sidehulls  in  a channel  of  width  W can  be  written  as  follows: 
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Considering  that  the  pressure  planform  is  rectangular  and  the  sidehulls 


» are  of  parabolic  shape,  the  above  integrals  can  be  evaluated  easily 

and  the  results  have  been  given  by: 
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where 


W. 
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total  weight  of  the  craft  = PQBL 
.bubble  width 

bubble  length 
sidewall  width 


3.3  Seal  Forces 

For  purposes  of  the  present  study  a very  simplified  seal  config- 
uration has  been  adopted.  Again,  this  has  been  done  in  order  to 
avoid  too  many  detail  points  which  would  reflect  a given  design 
rather  than  a general  craft. 

The  present  seals  are  assumed  to  be  flexible  fabric  seals,  such 
as  a bag  and  finger  design,  which  when  immersed  in  the  water 
simply  deflect  and  lie  on  the  water  surface.  Hence  they  do  not 
contribute  any  forces  or  moments  to  the  craft  except  for  their 
axial  drag  and  the  forces  and  moments  arising  due  to  the  shift 
in  the  center  of  pressure  of  the  air  in  the  plenum  caused  by 
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the  changing  imprint  length  on  the  water.  Referring  to  figure  3, 
which  shows  a simple  bow  seal,  the  following  equations  can  be 
derived : 


Z seal  = (p^  - p ) lw  B 

C a 

M seal  = Z seal  Is 


where  lw  = Is  tan  9 
sin  9b 

Is  = distance  of  seal  tip  to  C.G. 
9_  = sheer  angle  of  seal 

D 

9 = trim  of  craft 

the  axial  drag  due  to  the  seal  is: 


X seal  = Cf  - v “ IwB 

where  is  the  friction  factor,  derived  from  the  Reynolds 
number  as  follows: 


C£  = 0.044  (75) 

f ~R r'6 

e 

Rg  is  the  Reynolds  number  based  on  the  seal  wetted  length  lw. 

3.4  Aerodynamic  Forces 

The  aerodynamic  forces  and  moments  acting  on  the  craft  have 
been  simplified  and  are  represented  by  an  overall  drag  coefficient 
based  on  frontal  area  of  the  craft.  This  drag  coefficient  has 
been  selected  to  correspond  to  test  results  on  SES  type  con- 
figurations. The  force  is  therefore  simple: 

2 

X aero  = CD  }_  p V~Af  (76) 

2 

where  A^  = frontal  area  of  craft. 
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No  aerodynamic  lift  or  moments  have  been  used  in  the  present  study 
as  the  craft  is  operated  in  a straight  line  course  at  all  times 
and  no  wind  conditions  are  considered,  consequently  only  an  axial 
aerodynamic  force  exists. 

3.5  Propulsion  and  Control  Forces 

Various  methods  of  propelling  and  controlling  the  craft  exist. 
Current  emphasis  for  SES  propulsion  is  a water jet.  These  devices 
allow  for  thrust  vectoring  or  differential  thrust  for  maneuvering 
and  turning.  Since  we  are  only  considering  straight  line  oper- 
ation in  the  present  study  no  such  devices  will  be  considered. 

The  methodology  used  is  to  calculate  the  drag  for  a given  immersion 
and  trim  at  a specified  speed  and  allow  the  thrust  to  be  equal  to 
the  calculated  drag.  As  the  craft  responds  to  the  wave  environ- 
ment its  attitude  and  draft  will  be  varying  causinq  chances  in 

drag  and  cushion  pressure.  The  thrust,  however  is  kept  constant 
to  the  initially  calculated  value  thus  slight  surge  motions  will 
occur  during  the  time  the  craft  encounters  the  waves. 

3.6  Appendages 


Usually,  especially  in  the  case  of  an  SES,  directional  stabilizers 
or  fins  are  fitted  in  order  to  ensure  directional  stability. 
Standard  representations  of  these  appendages  exist  in  the  program 
code  developed,  which  account  for  drag  and  lift  forces.  In  the 
present  study  a nominal  fin  has  been  assumed. 

These  items  are  considered  as  base  vented  parabolic  sections 
designed  to  produce  the  required  lateral  stiffness  to  the  craft  to 
ensure  stability.  For  present  purposes,  therefore,  two  items 
attributing  to  the  total  drag  of  these  fins,  namely  pressure  drag 
and  frictional  drag,  are  considered. 


Since  the  quality  of  these  surfaces  has  to  be  kept  smooth  and 
constantly  clean  to  ensure  cavitation  free  operation,  it  is 
assumed  that  for  all  intents  and  purposes  the  surface  is  close 
to  being  hydrodynamically  smooth  and  consequently  the  frictional 
drag  is  computed  on  this  basis.  The  total  drag  of  the  stabilizer 
surface  can  then  be  written  as: 

Df in  = *5  pV2A[Cd  + 2 Cf]  (77) 

where 

2 

A = fin  surface  area,  ft 

= pressure  drag  coefficient 
= frictional  drag  coefficient 

For  a base  venting  parabolic  section  we  have 


where 

t = thickness-chord  ratio 

c 


and,  for  a smooth  surface  the  frictional  coefficient  can  again  be 
approximated  by  the  formula: 


Cf  = 0.044/Rc1//6 


(78) 


where 


Rc  = Reynolds  number  based  on  mean  chord. 

The  lift  force  from  the  fin  is  calculated  using  the  following 
classical  lift  equation. 
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(79) 


where 


AR 

V 


2 7T 


AR 

AR+3 


Aspect  ratio  of  the  fin. 


Relative  velocity  of  the  fin  in  the  water. 


4.  WAVE  ENVIRONMENT 

4 . 1 Wave  Representation 

The  computation  of  explosion  generated  waves  can  be  divided  into 
three  parts;  they  are  the  modeling  of  the  source  condition,  the 
calculation  of  propogation  and  transformation  of  waves  over  a 
given  bottom  topography,  and  the  determination  of  breaking  inception 
and  wave  run-up  according  to  some  acceptable  criteria.  The  last 
two  parts  would  involve  tedious  bookkeeping  of  propagation  history 
from  point  to  point,  should  the  bottom  topography  be  irregular. 

Since  the  study  in  the  current  phase  emphasizes  specifically  the 
mathematical  modeling  of  the  craft,  the  details  of  the  bottom 
irregularities  are  not  considered.  If  the  continental  shelf  is 
assumed  two-dimensional  and  to  have  a constant  mild  slope,  the 
wave  environment  can  simply  be  classified  into  two  characteristically 
different  groups;  deep  water  and  shallow  water  waves. 


4 . 2 Deep  Water  Wave  Generation 


The  deep  water  waves  theoretically  can  be  represented  by  sinusoids 
of  various  frequencies.  While  the  craft  responses  in  sinusoidal 
waves  are  to  provide  a general  indication  of  the  craft  characteristics 
as  a function  of  wave  period,  they  provide  little  information  as  to 
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how  the  craft  responds  when  it  is  sufficiently  close  to  the 
source  region,  as  the  wave  amplitudes  are  normally  very  large 
such  that  the  linear  superposition  technique  is  not  valid  and 
applicable.  The  present  model  is  capable  of  simulating  either 
a sinusoid  wave  system  or  an  idealized  explosion-generated  wave 
system  at  a given  stand-off  distance  from  the  source  at  any 
time  after  detonation.  Whereas  the  sinusoidal  wave  form  is 
simpler  and  well-known,  only  modeling  of  the  explosion-generated 
waves  is  discussed  in  the  following: 

The  problem  concerning  waves  generated  by  an  arbitrary  but 
localized  disturbance  on  a free  surface  has  been  investigated  by 
Kajiura,  reference  7.  In  analyzing  the  explosion-generated 
waves,  the  initial  disturbance  is  usually  assumed  as  being  of  a 
parabolic  crater-like  shape  with  radial  symmetry  such  that: 


n (r)  = 


nQ  [2(r/RQ)  2-H  for  r < Rq 


for  r > R 


where  n0  = crater  height 

R = crater  radius 

o 

r = radial  distance 


The  waves  resulting  from  this  disturbance  at  a distance  r from 
the  center  are  then  given  by,  Van  Dorn  et  al,  reference  8,  as 


, . . n R 

n (r,t)  = o o 


dv/dk 


] J3  (kRo)  cos  (kr-wt) 


where  k = wave  number,  determinable  from  the  relationship, 
between  the  group  velocity  V and  the  arrival  time 
t,  such  that 


V (k)  = H k (1+  sink  2kd)  = ^ 
w = / gk  tanh  kd 
d = water  depth 

J 3 = Bessel  function  of  the  1st  kind  of  order  3. 
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The  above  equation  shows  that  the  traveling  wave  train  possesses 
a series  of  amplitude  peaks  primarily  governed  by  the  moderating 
Bessel  function  J3.  The  problem  that  remains  is  to  relate  the 
crater  dimension  nQ  and  Rq  to  the  yield  of  a given  explosion  so 
that  prediction  of  waves  at  a given  location  r and  time  t can  be 
made . 


It  is  noted  that  both  n and  R are  not  easily  measurable.  What 
one  can  measure  are  the  wave  height  and  period  at  a large  distance 
from  the  source  disturbance.  It  is  in  fact  more  convenient  to 

measure  the  amplitude  peak  ’max  in  the  first  wave  envelop  at  a 

• k. 

given  range  r,  and  the  corresponding  wave  number  max  can  be 

evaluated  by  knowing  the  arrival  time  t from  the  above  equations. 
Analytically,  one  can  show  that,  for  a particular  source  distur- 
bance n (r) , the  amplitude  of  the  maximum  wave  nmax  is  inversely 
proportional  to  r,  and  the  corresponding  wave  number  max  is 
independent  of  the  crater  height  r)Q.  For  an  explosion  in 
sufficiently  deep  water,  ^max  can  be  determined  from  the  first 
stationary  value  of  J3  as: 


k 


max 


R 


0 


4.2 


Once  the  measurement  of  max  is  obtained,  the  crater  radius  can 
be  readily  estimated.  From  equation  (81),  one  finds: 


n R =1.63  nmax  r 

0 0 

when  k = max.  Consequently,  the  crater  height  can  also  be 
estimated  from  the  measurement  of  wave  height. 


Empirical  correlations  of  measurements  of  nmax  with  the  explosion 

yield  W and  the  detonation  depth  Z show  that  there  is  a certain 

„ 0 54  03 

trend  between  the  parameter  max  r/W  ' and  the  parameter  Z/W 

(W  in  lbs  of  TNT  equivalent) ; this  is  best  presented  graphically 
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by  plotting  the  experimental  data  points  as  shown  in  Figure  5. 

It  is  noted  that  there  are  two  peaks  appearing  in  the  former 

parameter  over  a range  of  the  latter.  One  of  these  peaks  occurs 

0 3 

at  Z/W  ’ = -0.05  and  is  commonly  termed  as  the  upper  critical  depth. 

Detonation  at  this  depth  is  seen  to  produce  the  highest  responses. 
The  other  peak  occurs  at  Z/w° . 3 = -2.7  and  is  usually  called  as 
the  lower  critical  depth. 

}c 

As  discussed  before,  the  parameter  max  can  be  determined  by 
measuring  the  arrival  time  of  the  first  wave  at  a given  distance. 

By  analyzing  the  wave  profiles  obtained  from  the  measurements, 
an  empirical  relationship  between  the  parameter  max  and  yield 
might  W has  also  been  established,  namely: 


max  =0.44  W 
= 0.39  w* 


_0  3 


for  0 > yw  ' 


- 0.25 


-0.25  > Z/W  * 5-7.5 


Using  these  empirical  relations  together  with  the  measured  results 
as  shown  in  Figure  5,  the  source  parameters  no  and  can  be 
determined  for  any  yield  at  any  water  depth  and  detonation  depth. 
Consequently,  the  wave  history  at  any  point  r and  time  t can  be 
calculated  according  to  equation  ( 81) . 

4.3  Shallow  Water  Waves 

Two  types  of  waves  should  be  considered  with  regard  to  shallow 
water  waves.  Firstly,  those  waves  produced  in  deep  water  as  a 
result  of  an  offshore  explosion  which  transform  their  height,  shape 
and  internal  characteristics  through  the  process  of  shoaling, 
refraction  and  reflection  when  they  propagate  shoreward  into 
shallower  water.  Secondly,  waves  directly  generated  by  explosions 
in  shallow  water  on  the  continental  shelf.  As  far  as  the  wave 
characteristics  are  concerned,  these  waves  can  be  considered 
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identical  and  be  treated  in  a similar  manner.  Before  entering 
into  the  discussion  of  how  to  model  these  waves  mathematically, 
however,  correlations  of  yield  with  wave  generation  in  shallow 
water  are  briefly  outlined  below. 


’I 


The  method  of  correlation  between  wave  heights  and  yields  dis- 
cussed in  the  previous  section  is  limited  to  deep  water  wave 
generation  such  that  d > 6 W0.3.  For  explosions  in  water  of 
depth  such  that  1 < /W° . 3 <6,  Le  Mehaute,  reference  9,  proposed 
a simple  interpolation  rule  to  fit  the  experimental  data,  as 
follows : 


n 


ndeep 


[H  + — ( 


l o 


d 

W°  . 3 


-1)  I 


(83) 


This  shows  that  the  generation  efficiency  is  reduced  by  half 
when  the  parameter  d/W0.3  approaches  unity.  In  the  case  of  very 
shallow  water  where  d/W0.3  <<1,  the  linear  model  is  no  longer 
valid  and  different  correlations  must  be  used.  Unfortunately, 
there  are  very  few  data  collected  of  shallow  water  explosions. 
Among  the  available  data  as  listed  in  Table  2,  only  the  WES 
test  data,  reference  10,  provide  systematic  variations  of  charge 
weight  and  water  depth. 


By  means  of  small-scale  charges  (0.5  -2048  lbs.)  the  WES  program 
was  designed  to  estimate  wave  effects  due  to  a 20  KT  explosion 
in  water  of  30  to  200  feet  deep.  The  charge  position  varied 
from  beneath  the  bottom  to  above  the  free  surface.  The  results 
showed  that  variations  of  Z/d  from  -1.0  to  0 had  little  effect 
on  wave  height.  In  contrast  to  deep  water  explosions,  the  most 
significant  parameter  for  wave  generation  in  shallow  water  is 
water  depth,  instead  of  charge  position. 

The  other  significant  feature  is  that  the  dispersion  law  is 
different  for  waves  propagating  in  deep  and  shallow  water.  In 
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deep  water,  wave  height  varies  inversely  with  radial  distance  r 
as  a combined  result  of  frequency  and  radial  dispersions.  In 
extremely  shallow  water,  the  large  leading  wave  is  expected  to 
behave  like  a solitary  wave  so  that  its  height  varies  inversely 
as  r2/3  instead  of  r.  In  moderately  shallow  water,  the  relation 
below  should  hold 

n r“  = constant  2/3  £ 8(d)  £ 1 (84) 


In  correlating  the  WES  test  data,  the  following  empirical  formula 
is  derived 


nmax  r 


W6/ 3+0.25 


= 1-44  (£i/3)°.93 


where  3 = 0. 83  (^ /3) 


0 07 


(85) 


It  is  noted  that  the  power  8 varies  as  a function  of  the  depth 
parameter  d/W1/3;  for  the  very  shallow  case,  8 approaches  2/3 
as  a limit.  While  the  derivation  of  the  above  realtionship  has 
assumed  that  reasonable  extrapolation  of  the  WES  data  is  valid, 
it  must  be  noted  that  the  correlation  is  based  upon  the  experi- 
mental data  covering  d/W1 /3  up  to  0.585.  There  is  no  indication 
that  it  will  approach  the  empirical  relation  (83)  as  d increases. 


Equation  (85)  provides  an  empirical  relationship  for  predicting 
the  maximum  wave  height  at  any  distance  r from  a shallow  water 
explosion.  After  the  wave  height  is  determined  for  a given 
explosion,  the  important  procedure  required  for  numerical  simu- 
lation is  a mathematical  representation  of  the  wave  history  as 
a function  of  time.  As  mentioned  earlier,  disregarding  whether 
the  waves  are  generated  in  shallow  water  or  are  propagated  into 
shallow  water  from  offshore,  their  internal  characteristics  are 
approximated  identical  if  their  height  and  period  are  the  same. 
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The  most  important  parameter  which  affects  these  waves  in  this 
case  is  the  local  water  depth.  As  is  well  known,  when  waves 
propagate  into  shallower  water,  their  crests  become  more  peaked 
through  shoaling.  When  the  local  depth  d becomes  so  shallow 
that  the  wave  height  h = 0.67  d to  0.78  d,  waves  start  to  break. 
Analytical  and  experimental  studies  of  wave  propagation  and 
transformation  have  been  discussed  in  detail  by  Le  Mehaute  et  al, 
reference  13,  and  Divoky  et  al,  reference  14.  Their  analyses 
show  that,  among  many  existing  wave  theories,  the  cnoidal  wave 
theory  is  good  for  describing  the  transition  from  deep  water 
waves  to  shallow  water  waves  but  the  solitary  wave  theory  best 
describes  the  long,  shallow  water  waves  including  the  spilling 
type  breakers.  In  the  present  study,  the  solitary  wave  form  is 
used  for  numerical  modeling  of  the  long  period  waves  on  the 
continental  shelf.  After  the  wave  height  and  period  is  determined 
according  to  the  yield  weight,  the  mathematical  representation 
of  the  waves  in  water  of  depth  d is  given  by 

n (x,t)  = h sech2  a(x-ct)  (86) 

where 

h = wave  height 

a = V 3h/  4d 
c = /gd  (1+  h/2d) 

5.  NUMERICAL  TECHNIQUES 

The  computer  program  developed  to  integrate  the  equations  of 
motions  under  the  influence  of  the  forces  and  moments  inputed  to 
the  craft  by  the  wave  environments  described  previously  will 
now  be  briefly  discussed. 

Initiation  of  the  computation  is  made  by  entering  into  the  program 
the  initial  conditions  of  the  craft  such  as  altitude,  speed  and 
craft  weight.  Overall  craft  dimensions  and  the  geometry  of  the 


sidehulls  and  seals  are  also  required.  With  the  above  information 
the  submerged  geometry  of  the  sidehulls  and  seals  are  calculated 
and  the  forces  and  moments  from  all  sources  described  in  section  3 
are  calculated.  The  initial  values  of  all  the  variables  are 
then  used  as  starting  values  at  time  t=0,  to  initiate  integration 
of  the  equations  of  motion. 


A fourth  order  Runge  Kutta  scheme  is  used  for  the  integration 
and  all  variables  are  updated  at  each  time  step.  Time  steps  of 
the  order  of  0.1  seconds  are  normally  used. 

The  wave  field  is  also  started  at  t=0  and,  depending  on  the  par- 
ticular case  being  considered,  propagates  towards  the  craft  in 
a predetermined  direction  (heading)  . Three  options  are  currently 
available  for  the  wave  field,  as  previously  discussed.  Experienced 
gained  with  the  program  indicates  that  runs  vary  typically  from 
about  15  seconds  to  70  seconds  when  using  a CDC  6600. 

An  overview  flow  chart  illustrating  the  general  operations  per- 
formed in  the  computer  is  shown  in  figure  5. 

6 . RESULTS 

The  program  was  exercised  under  various  wave  conditions  and  craft 
headings  to  investigate,  on  a preliminary  basis,  the  response  of 
a typical  SES . In  order  to  conduct  this  study  several  assumptions 
had  to  be  made  regarding  the  craft  size  and  dimensions.  In  order 
to  make  the  results  relevant  to  current  interests  an  SES  having 
characteristics  similar  to  the  2000  ton  class  was  chosen.  Some 
of  the  salient  features  of  this  craft  are  listed  below: 


Craft  weight  = 2000  lng  tons 
Cushion  length  = 240  feet 
Cushion  beam  = 88  feet 
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Center  of  Gravity  location  = 130  ft.  forward  of  transom 

24  ft.  above  keel 

Lift  fan  characteristic  : 

Qf  = 75,497  - 121  (pc  - pa)  cfs 

Bow  seal  angle  = 30° 

Stern  seal  angle  = 60° 

Initial  air  leakage  area  = 49  ft.2 

The  definitions  of  inputs  and  a sample  input  case  are  shown  in 
Appendix  A.  This  input  provides  further  information,  including 
details  of  the  sidehull  shapes  chosen.  This  shape  is  also 
representative  of  typical  sidehull  designs  for  SES . 

The  results  of  these  runs  are  shown  in  the  following  figures. 

6.1  Sinusoidal  Wave  Response 

In  order  to  exercise  the  program  and  obtain  a reference  base 
of  craft  response  a series  of  runs  was  conducted  using  a 
sinusoidal  wave  excitation.  This  wave  was  chosen  to  correspond 
to  the  significant  wave  height  and  period  of  a Sea  State  3 
Pierson  Moskowitz  spectrum,  namely: 

Wave  height  = 5 feet 
Period  = 6 seconds 

1 

Figures  6 through  9 illustrate  the  results  of  these  runs  for 
heading  angles  of  0°,  45°,  165°  and  180°,  respectively.  In 
, these  figures  the  cushion  pressure  and  wave  profile  are  shown 

in  the  upper  figure;  the  pitch  and  yaw  in  the  middle  and  the 
heave  response  and  roll  in  the  lower  diagram.  The  curves  are 
shown  as  a function  of  a non-dimensional  time,  t/T.  The  re- 
quired conversion  factor  to  real  time  is  given  in  each  caption. 
For  these  specific  runs  the  craft  immersion  at  the  center  of 
gravity  was  1.9  feet  at  a initial  trim  of  0.9  degrees. 
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Figure  6 Craft  Response  in  Sea-State  3 - average  wave  height  = 3 ft 

average  wave  period  = 6 sec,  craft  speed  = 80  knots,  craft 
heading  = 0 deg  * 1>10  sec 


Figure  7 Craft  Response  in  Sea-State  3 - average  wave  height 

average  wave  period  = 6 sec,  craft  speed  = 80  knots 
craft  heading  = 45  deg  1 = 0.89  sec 


It  will  be  seen  from  these  results  that  generally  the  craft  is 
well  behaved  in  this  sea  state.  However,  in  the  cases  of  165° 
and  180°  heading  fairly  large  oscillations  occur  in  cushion 
pressure  with  corresponding  heave  excursions.  These  two  runs, 
which  deal  with  head  seas,  illustrate  that  the  cushion 
pressure  is  certainly  more  responsive  to  head  seas  than  to 
the  following  seas,  shown  in  figures  6 and  7. 

It  should  be  emphasized,  however,  that  no  heave  alleviation 
devices  have  been  modelled  in  the  present  program  and  conse- 
quently this  type  of  behavior  is  not  unexpected.  It  is  antici- 
pated that  installation  of  such  a device  would  alleviate  this 
situation  considerably. 

The  cases  with  quartering  seas,  namely  figures  7 and  8 show  roll 
and  yaw  responses.  In  particular  figure  7 shows  a roll  amplitude 
of  approximately  2°.  This  case,  run  for  a quartering  following  sea, 
also  shows  an  increased  pitch  response  which  is  to  be  expected. 

The  wave  length  to  cushion  length  ratio  for  these  runs  is  1.29 
which  is  removed  from  the  wave  pumping  value  of  2.0.  It  should 
also  be  pointed  out  that  the  natural  frequency  of  the  craft  under 
the  above  initial  trim  and  heave  conditions  is  approximately: 

Pitch  natural  frequency  = 1.45  Hertz 

Heave  natural  frequency  = 1.67  Hertz 

6.2  Solitary  Wave  Response 

As  discussed  in  section  5 , within  the  continental  margins  in 
shallow  water  , the  waves  caused  by  deep  water  explosion  can  be 
represented  by  solitary  waves.  In  this  section  we  have  con- 
ducted a series  of  runs  wherein  the  craft  response  to  a solitary 
wave  at  various  headings  has  been  investigated.  Furthermore, 
the  effect  of  varying  water  depth  and  wave  height  is  also  shown. 
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For  the  present  runs  the  initial  trim  and  center  of  gravity- 
immersion  was  taken  to  be  1 degree  and  2 feet,  respectively. 
The  runs  were  performed  at  a craft  speed  of  50  knots,  except 


for  one  run  which  was  conducted  in  a near  hovering  mode.  (actual 
speed  was  5 knots.)  The  wave  height  and  period  were  varied 
during  this  series  and  are  identified  in  the  caption  of  each 
figure . 

The  hovering  condition  is  shown  in  figure  10.  As  will  be  seen 
in  this  run  the  water  depth  is  taken  as  60  feet  with  a wave 
period  of  15  seconds  and  wave  height  of  6 feet.  With  these 
conditions  the  ratio  of  wave  length  to  cushion  length  is  2.91. 
Behavior  of  the  craft  is  quite  acceptable  with  the  maximum 
pitch  and  heave  excursions  shown  in  Table  3. 

The  effects  of  varying  heading  angle  for  the  conditions  described 
in  the  above  case  are  shown  in  figures  11,  12  and  13.  As  will 
be  seen  from  these  runs,  all  at  50  knots,  the  pitch  excursions 
increase  as  the  heading  varies  from  a beam  sea  condition  to  a 
head  sea.  Attendant  with  this  change  in  heading  the  roll 
and  yaw  decreases.  In  the  case  of  90°  heading  or  beam  seas  the 
roll  motion  is  excited  at  a natural  frequency  of  about  1.34  Hertz. 

It  is  apparent  from  these  curves  that  the  craft  will  survive  this 
wave  environment  without  undue  difficulty.  It  should  be  pointed 
out  that  whereas  it  may  appear  desirable  to  head  away  from  a 
blast  situation,  should  one  occur,  the  results  would  depend  on 
the  craft  relative  velocity  to  the  wave.  This  approach  for  reducing 
wave  induced  damage  may  or  may  not  be  appropriate. 

Figures  14,  15  and  16  illustrate  the  behavior  of  the  craft  under 
differing  combinations  of  wave  height,  water  depth  and  wave  period 
for  a head  sea  i.e.,  heading  of  180°.  The  exact  conditions 
are  given  in  Table  3 . In  the  first  two  cases,  figures  14  and 
15,  craft  response  is  reasonable  although  some  relatively  large 
heave  and  pitch  excursions  occur  in  the  10  foot  wave  condition. 
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Figure  10  Craft  Response  In  Shallow  Water  - water  depth  = 60  ft 
wave  period  = 15  sec,  wave  height  = 6 ft,  craft  speed 
(hovering),  craft  heading  = 0 deg,  JL  = 3.82  sec 
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Solitary  Wave  Results 


Period  = 15 

sec . 

Wave  Height 

= 6 ft.  Depth  = 60  ft. 

X/L  = 2.91 

Heading 

(Deg) 



Speed 

(Knots) 

Maximum 

BUM 

Trim  (Deg) 

Roll  (Deg) 

Yaw  (Deg) 

0 

0 

-3.61 

0 

0 

-4.7 

90 

50 

-0.31 

+5.10 

-2.7 

-5.0 

135 

-1.97 

+4.19 

-0.89 

-4.48 

180 



so 

-3.69 

0 

0 

-4.72 

Period  - 30 

sec. 

Wave  Height 

= 10  ft. Depth  = 60  ft. 

X/L  = 6.01 

180 

— 

50 

-5.37 

0 

0 

r 

-8.08 

Period  = 30 

sec . 

Wave  Height 

= 5 ft.  Depth  = 30  ft. 
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= 6.5  ft. Depth  = 30  ft. 

X/L  = 2.18 
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0 
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Note:  Heading  angle  defined  as  0°  for  following  seas, 

180°  for  head  seas. 
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craft  heading  = 90  deg,  t/T  = 3.12  sec 


Figure  12  Craft  Response  in  Shallow  Water,  water  depth  = 60  ft, 

wave  period  = 15  sec,  wave  height  = 6 ft,  craft  speed  = 50  knots, 
craft  heading  = 135  deg,  iL  = j.36  sec 
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The  normal  reaction  of  the  craft  on  encountering  the  wave  front 
is  to  increase  its  trim  with  a simultaneous  decrease  in  immersion. 
This  behavior  is  to  be  anticipated  since  the  craft  is  impacting 
the  wave.  In  accordance  with  this  reaction  the  cushion  pressure 
decreases  due  to  increase  leakage.  After  the  wave  crest,  trim 
decreases  and  immersion  increases. 

i 

In  figure  16  it  is  apparent  that  large  excursions  in  pitch  and 
heave  are  occurring  and  furthermore  these  motions  are  diverging. 
This  particular  run  condition  however,  has  been  taken  at  a ratio 
of  wave  length  to  cushion  length  of  2.18,  which  is  very  close  to 
the  wave  pumping  condition  of  2.  Therefore  it  is  expected  that 
severe  conditions  will  arise,  as  seen  in  Table  3.  The  maximum 
heave  and  pitch  are  larger  in  relation  to  the  wave  amplitude  than 
in  all  other  cases.  It  is  apparent  that  under  this  condition  the 
craft  is  not  likely  to  survive  without  evasive  action. 

6.3  Deep  Water  Explosion  Wave  Response 

A deep  water  explosion  wave  environment  was  generated  by  con- 
sidering the  wave  caused  by  a yield  device  of  1 kiloton. 

It  was  assumed  that  the  stand-off  distance  of  the  craft  from 
the  center  of  the  blast  was  7,500  feet.  A device  having  this 
yield  and  exploding  at  the  upper  critical  depth  would  cause  an 
initial  water  disturbance  having  an  initial  radius  of  740  feet 
and  initial  water  height  of  55  feet.  Using  these  conditions 
the  craft  response  was  calculated  for  a hovering  situation  at 
0 degree  heading,  and  at  50  knots  for  headings  of  90°,  (beam 
, sea),  45°  and  135°.  The  results  of  these  runs  are  shown  in 

figures  17  through  20. 

From  figure  17  it  is  seen  that  at  this  distance  from  the  blast, 
the  maximum  wave  height  encountered  is  approximately  6 feet. 

The  graph  shows  the  arrival  of  the  first  wave  group  and  the 
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subsequent  response  of  the  craft.  As  seen  all  the  variables  are 
within  normal  excursions  with  a maximum  pitch  of  -1.55°  and  heave 
of  2.85  feet.  The  wave  envelop  shown  in  this  figure  is  typical 
of  the  explosion  generated  wave  envelopes. 

Should  a blast  occur  off  the  beam  of  the  craft  when  operating 
at  50  knots  the  results  indicate  that  the  craft  will  barely 
survive  the  waves.  As  seen  in  figure  18  large  excursions  in 
yaw,  roll  and  heave  are  experienced.  The  maximum  excursion  in 
these  variables  are: 

maximum  heave  = 4.85  feet 
maximum  roll  = 8.98° 
maximum  yaw  = -11.91° 

The  maximum  pitch  angle  experienced  is  -1.52°  which  is  quite 
nominal.  It  is  apparent  from  this  response  that  the  craft  is 
quite  vulnerable  to  beam  explosions. 

Should  the  craft  be  operating  at  50  knots  and  an  explosion  occur 
the  question  arises  as  to  what  evasive  action  it  should  take. 

As  a preliminary  maneuver  it  was  assumed  that  a reaction  time  of 
80  seconds  is  required  for  the  craft  to  either  alter  its  course 
to  another  heading  or  head  up  into  the  blast  and  kill  its  engines. 
We  have  seen  that  in  this  latter  mode  it  can  survive  the  present 
explosion.  The  question  arises  as  to  whether  an  alternative 
course  heading  is  preferable.  To  investigate  this  possibility 
two  headings  of  45°  and  135°  were  investigated. 

Figure  19  shows  the  response  of  the  craft  to  the  waves  environ- 
ment on  a heading  of  45°  assuming  such  a heading  is  achieved 
80  seconds  after  the  blast.  As  will  be  seen  little  if  any 
motion  occurs  to  the  craft  since  the  craft  is  heading  away  from 
the  wave  front  and  is  apparently  in  small,  long  period  waves 
ahead  of  the  main  group  of  waves.  Provided  sufficient  clear  sea 
is  available  the  craft  could  outrun  the  wave  until  the  waves 
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had  decayed  sufficiently  to  allow  a change  in  heading.  This 
situation  would  obviously  be  changed  if  the  continental  margin 
were  reached,  since  in  that  case  the  waves  would  begin  to 
experience  bottom  effects  and  become  solitary  waves. 

Should  the  craft  head  into  the  blast  on  a 135°  course,  an  un- 
likely situation,  unless  it  was  already  on  this  course  when  the 
blast  occurred  the  response  is  shown  in  figure  20.  Here  it  will 
be  seen  that  the  motions  are  diverging  and  indeed,  based  on 
the  present  analysis,  the  craft  will  not  survive.  As  will  be 
seen  the  run  was  actually  terminated  before  the  motions  become 
excessive.  The  maximum  excursions  of  the  craft  at  this  point 
are : 

maximum  heave  = 4.526 
maximum  pitch  = 1.548 
maximum  roll  = 2.325 
maximum  yaw  = -0.331 

It  is  apparent  from  this  brief  survey  that  dependent  on  the 
location  of  the  blast  relative  to  the  craft  and  the  available 
response  time  several  possible  scenarios  exist  for  evasive 
action  subsequent  to  a blast.  This  evasive  action,  however 
depends  on  a great  many  variables  and  will  be  the  topic  of 
further  investigation  during  the  next  phase  of  the  present 
program.  It  it  clear  however,  that  relatively  moderate  yields 
can  cause  an  SES  considerable  difficulty  if  cognizance  of  the 
seriousness  of  the  situation  is  not  realized. 

It  should  be  reiterated  at  this  point  that  the  above  analyses 
were  performed  without  heave  alleviation  devices  on  the  craft, 
results  should  therefore  be  viewed  in  this  light. 


The 
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7. 


CONCLUSIONS 


We  have  presented  the  preliminary  results  derived  from  modelling 
the  response  of  a typical  SES  to  an  explosion  type  wave  environ- 
ment. This  effort  represents  the  status  of  the  project  at  the 
end  of  its  first  phase  and  is  by  no  means  complete.  Results 
obtained  to  date  however,  indicate  that  some  consideration  is 
needed  to  developing  suitable  tactical  maneuvers  for  minimizing 
damage  potential  from  explosion  generated  wave  environments. 

It  is  further  apparent  that  the  unique  high  speed  capability  of 
these  vehicles  provides  them  with  a great  deal  of  versatility 
in  generating  evasive  action. 

It  is  presently  planned  in  the  next  phase  of  the  project  to 
take  the  procedure  further  with  the  development  of  the  computer 
code  to  enhance  its  versatility  and  to  perform  a parametric 
analysis  on  the  various  variables  of  the  problem,  in  order  to 
develop  a meaningful  procedure  for  action  under  possible  threats 
of  this  nature.  Such  scenarios  would  be  most  valuable  in  future 
planning  of  the  operational  fleet. 
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APPENDIX  A 


COMPUTER  INPUT/OUTPUT  FORMAT 


f 


Card  2:  Format  (F10.0,  1415) 

1)  DT  - Number  of  intervals  per  wave  period. 

2)  NSTEP  - Number  of  integration  steps. 

3)  NPRNT  - 1 

4)  IP  - Debug  flag  for  component  forces  and  moments,  printed 

in  main  program. 

If  IP  = 0,  debug  not  printed. 

If  IP  ^ 0,  debug  is  printed. 

5)  IFIN  - Flag  on  inclusion  of  stabilizer. 

If  IFIN  = 0 , do  not  include  stabilizer. 

If  IFIN  ^ 0,  include  stabilizer. 

6)  IPLOT  - Flag  on  plotting. 

If  IPLOT  = 0,  call  PLOTT. 

If  IPLOT  = 1,  call  PLOTXY . 

If  IPLOT  = 2,  call  PLOTT  and  PLOTXY. 

If  IPLOT  > 2 , do  not  plot. 

7)  IPT  - Number  of  points  to  plot  for  PLOTT. 

8)  NJET  - Number  of  jets  for  thrust  vector  control. 

9)  INT  - Flag  for  printing  cumulative  integrals  and  geometrical 

variables . 

If  INT  = 0,  do  not  print. 

If  INT  ^ 0,  print. 

10)  IBUG  - Flag  on  debug  for  subroutine  BUOY. 

If  IBUG  = 0,  do  not  print  debug. 

If  IBUG  f-  0,  print  debug. 

11)  IW  - Flag  on  wave  type. 

If  IW  = 1,  sinusoidal  wave. 

If  IW  = 2,  solitary  wave 
If  IW  = 3,  explosion  wave 

12)  IPR  - Flag  on  debug  for  pressure  subroutine 

If  IPR  = 0,  do  not  print  debug. 

If  IPR  ? 0,  print  debug. 

13)  ICO  - Flag  for  generating  new  derivatives  when  draft  changes 

by  more  than  ICO  feet. 

If  ICO  = 0,  do  not  change  derivatives. 

If  ICO  > 0,  ICO  equals  the  change  in  draft  required  u 
update  derivatives. 


7) 


Card  3:  Format  (8F10.0) 


1}  THT  - Initial  pitch  angle  of  craft  (deg) . 

2)  PHI  - Initial  roll  angle  of  craft  (deg) . 

3)  PSI  - Initial  yaw  angle  of  craft  (deg) . 

4)  Z - Heave  (set  = 0)  . 

Card  4:  Format  (8F10.0) 


1)  AA  - Distance  from  transom  to  C.G.  (ft)  . 

2)  BB  - Half  spacing  of  side  walls  (ft) . 

3)  CC  - Side  wall  immersion  at  C.G.  (ft). 

4)  DD  - Distance  from  keel  of  craft  to  C.G.  (ft). 

5)  AM  - Craft  weight  (tons) . 

6)  DXDU  - Added  mass  coefficient  of  side  wall  in  axial  flow. 

7)  AIX  - Moment  of  inertia  of  craft  about  the  x-axis  ( ton-ft-sec2) 

8)  AIZ  - Moment  of  inertia  of  craft  about  the  z-axis  (ton-ft-sec2) 

9)  AIY  - Moment  of  inertia  of  craft  about  the  y-axis  (ton-ft-sec2) 

Card  5:  Format  (8F10.0) 

1)  WL  - Reference  length  of  craft  (ft) . 

2)  SP  - Approaching  speed  (knots). 

3)  RHO  - Density  of  water  (lb.  sec  /ft4) 

4)  ANU  - Kinematic  viscosity  of  water  (ft2/sec). 

5)  CDLL  - Drag  coefficient,  lateral  force,  lateral  motion. 

6)  CDNN  - Drag  coefficient,  normal  force,  normal  motion. 


Card  6:  Format  (8F10.0) 

1)  OMEGA  - Dihedral  angle  of  stabilizer  (deg) . 

2)  CR  - Chord  length  of  stabilizer  at  root  (ft) . 

3)  CT  - Chord  length  of  stabilizer  at  tip  (ft). 

4)  S - Stabilizer  span  (ft). 

Card  7:  Format  (8F10.0) 

1)  CCO  - Side  wall  immersion  at  C.G.  before  turning  (ft). 

2)  THTO  - Pitch  angle  before  turning  (deg) . 

3)  SPTURN  - Assigned  speed  at  turn  if  different  from  SP  (knots) 

4)  DFTH  - Control  for  differential  thrust  (set  = 0) . 

Card  8:  Format  (8F10.0) 


XARM  - Longitudinal  distance  of  water  jet  nozzle  location 
from  craft  C.G.  (ft). 

ZARM  - Vertical  distance  of  water  jet  nozzle  location  below 
craft  C.G.  ( ft) . 

BACE  - Vertical  location  of  the  stabilizer  attachment  belo\ 
the  keel  line  (ft). 
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Card  9:  Format  (8F10.0) 


Ik 


1)  YARM(I)  - Transverse  location  of  Ith  water  jet  nozzle  from 
craft  centerline  (ft) 

NJET  values.  Positive  starboard  side. 

Card  10:  Format  (8F10.0) 

1)  DELJET ( I ) - Deflection  angle  of  nozzle  I (deg). 

NJET  values.  Positive  toward  port  side. 

Card  11:  Format  (8F10.0) 

1)  RMCP(I)  - Engine  power  level  delivered  to  nozzle  I. 

NJET  values. 

Card  12:  Format  (8F10.0) 

1)  ALPHA ( I ) - Vertical  tilt  angle  of  nozzle  I (deg). 

NJET  values.  Positive  upward. 

Card  13:  Format  (8F10.0) 

1)  DWET  - Distance  from  keel  to  wet  deck  (ft). 

2)  WAMP  - Wave  amplitude  (ft) . 

3)  WPER  - Wave  period  (sec.). 

4)  BETA  - Heading  angle  (deg) 

BETA  = 0,°  following  or  overtaking  waves. 

BETA  = 180°,  head  waves. 

5)  WDEP  - Water  depth  (ft) . 

6)  XO  - Distance  from  center  of  explosion  to  craft  (ft). 

7)  RO  - Crater  radius  (ft). 

8)  ETAO  - Crater  height  (ft). 

9)  TO  - Reference  time  with  respect  to  time  of  detonation  (sec) . 
Card  14:  Format  (1615) 

1)  NST  - Number  of  sections  along  craft  from  transom  to  bow. 

Card  15:  Format  (8F10.0) 

1)  BUBL  - Air  cushion  bubble  length  (ft). 

2)  BUBB  - Air  cushion  bubble  width  (ft). 

3)  WALB  - Maximum  width  of  side  wall  (ft) . 

4)  DEPTH  - Depth  of  craft  (ft). 

Card  16:  Format  (8F10.0) 

1)  SLBOW  - Length  of  planing  bow  seal  (ft) . 

2)  SLSTRN  - Length  of  planing  stern  seal  (ft). 

3)  THETA  - Angle  of  planing  seal  (deg) . 

I 
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Card  17:  Format  (8F10.0) 


1)  DRISE(I)  - Dead  rise  angle  at  station  I (deg). 

NST  values.  I = 1 at  transom,  I = NST  at  bow. 

Card  18:  Format  (8F10.0) 

1)  ENTRCE ( I ) - Average  entrance  angle  at  station  I (deg). 

NST  values. 

Card  19:  Format  (8F10.0) 

1)  CHINE (I)  - Height  of  chine  above  keel  line  at  station  I (ft). 
Card  20:  Format  (8F10.0) 

1)  NSW ( I ) - Number  of  water  lines  used  for  defining  offsets 
at  station  I. 

Card  21:  Format  (8F10.0) 

1)  XSW(I)  - Distance  from  transom  to  station  I (ft).  NST  values. 
Card  22:  Format  (8F10.0) 

1)  HSW(I)  - Height  of  bottom  profile  above  keel  line  at 

station  I (ft).  If  profile  below  keel  line  HSW(I) 
is  negative. 

Card  Group  23:  Format  (8F10.0) 

1)  D1 ( I , J)  - Height  of  Jth  waterline  above  keel  at  Ith 

station  (ft).  NSW(I)  values  of  Dl  for  each  I. 

All  values  are  positive.  Dl(I,l)  = 0.0.  Refer 
to  figure  A.l. 

Dl  is  input  as  follows: 

Card  1 - Dl(l, I),  Dl(l,2),  ...  Dl(l,  NSW(l)). 

Card  2 - Dl(2,l),  Dl(2,2),  ...  Dl(2,  NSW(2)). 

Card  NST  - D1(NST,1),  Dl(NST,2),  ...  Dl (NST, NSW(NST) ) . 

Card  Group  24:  Format  (8F10.0) 

1)  W1(I,J)  - Horizontal  offset  of  the  starboard  wall,  right 

side  of  vertical  reference  plane,  at  Ith  station 
and  Jth  waterline  (ft).  NSW(I)  values  of  Wl  for 
each  I.  All  values  are  positive.  W1(I,J)  input 
similarly  to  D1(I,J).  Refer  to  figure  A.l. 
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Card  Group  25:  Format  (8F10.0) 


i 

i 


I 


I 


1)  W2(I,J)  - Horizontal  offset  of  the  port  wall,  left 

side  of  vertical  reference  plane,  at  Ith  station 
and  Jth  waterline  (ft).  NSW(I)  values  of  V\2  for 
each  I.  All  values  are  positive.  W2(I,J)  input 
similarly  to  Dl(I,J).  Refer  to  figure  A.l. 

The  following  nondimensional  derivatives  may  be  obtained  by 
calculation  or  from  model  tests. 

Card  26:  Format  (8F10.0) 


1) 

DYV  - 

y ' 

J V 

2) 

DYP  - 

y ' 

1 P 

3) 

DYR  - 

y' 

1 r 

4) 

DYDV 

- y'  . 

5) 

DYDP 

- y ' . 

* P 

6) 

DYDR 

- y 1 . 

1 r 

Card  27:  Format  (8F10.0) 


1)  DKV  - k ' 

2)  DKP  - k ' 

3)  DKR  - k' 

r 

4 ) DKDV  - k ' . 

v 

5)  DKDP  - k' . 

P 

6)  DKDR  - k' . 

r 

Card  28:  Format  (8F10.0) 

1)  DNV  - n' 

v 

2)  DNP  - n ' 

3)  DNR  - n ' 

4)  DNDV  - n ' ^ 

5)  DNDP  - n ' ^ 

6)  DNDR  - n ' £ 


If  CC  ^ CCO  a second  set  of  nondimensional  derivatives  are  read  as 


cards 

29 

, 30  and 

31 

• 

Card 

29: 

Similar 

to 

card 

26 

for 

CC, 

THT 

case 

Card 

30: 

Similar 

to 

card 

27 

for 

CC, 

THT 

case 

Card 

31: 

Similar 

to 

card 

28 

for 

CC, 

THT 

case 

Definition  of  Output 


Input  data  are  reproduced  as  they  appear  on  data  cards, 
with  the  exception  of  Dl,  Wl,  W2  which  are  not  printed 
in  the  order  they  are  read. 


Any  input  data  that  are  converted  in  the  program  are 
printed  in  new  units. 


1)  SP  ( ft/sec)  . 

2)  AM  (non-dimensional) . 

3)  AIX  (non-dimensional) . 

4)  AIY  (non-dimensional) . 

5)  AIZ  (non-dimensional) 

6)  FROUDE  (non-dimensional)  - Froude  number 

3.  Craft  attitude 

1)  Draft  (ft). 

2)  Trim  (deg) . 

4.  Non-dimensional  derivatives  printed  as  read  from  input. 

5.  Stabilizer  coefficients. 

6.  Coefficients  for  ship  plus  stabilizer. 

7.  Stability  criterion  for  ship  only. 

8.  Stability  criterion  for  ship  plus  stabilizer. 

9.  Center  of  pressure  of  sidewall. 

10.  Non-dimensional  cumulative  integrals. 


I , DdF 

jPt4S 


I 


DF2I 


DF3.T 


\ DF  dF 

S 

f 2 

t DF  dF 

J pi:S 
f 3 

! DF  clF 

^ p &s 


DC  dF 
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DC  21 


DC?'dF 


DC  31 


DC  FI 


- i DC3dF 

Jp&s 

- f DCF  dF 

Jp&ts 


DCF2I  - 


DC  2 FI 


B3BI 


DCF2dF 


f DC2FdF 
Jp&c  s 

3 f 

BB  \ B dF 
Jp«s 


where 


p&s  - Integration  limits  over  both  port  and  starboard  sidewalls. 

D - Draft  at  successive  stations . 

F - Distance  from  C.G.  to  successive  stations. 

C - Vertical  moment  arm,  at  successive  stations,  for  submerged 

portions  of  craft  (ft). 

B - Bearn  at  successive  stations. 

BB  - Half  spacing  of  side  walls. 

11.  Non-dimensional  Geometrical  Variables  as  Function  of  Roll. 

GI(I)  - Integral  of  girder. 

SI(.I)  - Integral  of  cross  sectional  area. 

S1(I)  - Cross  sectional  area  at  transom. 

TDRAF(I)  - Draft  at  transom. 

If  CC^  CCO  output  from  (3)  to  (11)  will  be  printed  for  new  craft 
attitude  corresponding  to  CC  and  THT. 

12.  Craft  characteristics . 

13.  Wave  characteristics . 
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14.  Table  of  output  plus  units: 

1)  T - Time  (sec) . 

2)  U - Craft  speed  (knots) 

3)  BETA  - Sideslip  angle  (deg) . 

4)  W - Heave  rate  (ft/sec) . 

5)  X 

6)  Y Location  of  craft  (craft  lengths) . 

7)  Z - Heave  (ft) . 

8)  PHI  - Roll  angle  (deg) . 

9)  THETA  - Pitch  angle  (deg) . 

10)  PSI  - Yaw  angle  (deg) . 

11)  PC  - Cushion  pressure,  gage  pressure  (psf) . 

12)  QF  - Fan  flow  (ft3/sec) 

13)  QO  - Leakage  flow,  difference  from  initial  condition  (ftJ/sec). 

14)  VOLDOT  - Rate  of  cushion  volume  variation  (ft3/sec) . 

15)  WD  - Heave  acceleration  (G's). 

16)  WH  - Wave  elevation  at  C.G.  of  craft  (ft). 

17)  VOL  - Cushion  volume  (ft3). 

15.  Legend  for  computer  plots: 

1)  HEAVE  - Heave  (ft) . 

2)  ROLL  - Roll  angle  (deg) . 

3)  PITCH  - Pitch  angle  (deg) . 

4)  YAW  - Yaw  angle  (deg) . 

5)  WAV  HGT  - Wave  elevation  at  C.G.  of  craft  (ft) . 

6)  PRESS/100  - Cushion  pressure  divided  by  100  (pst) . 
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i 


OflApTs  2.U0 


TRIMs  1.00 


( 


t 


N/JN  JiMfeNSiONAL  nTHIVATIVES 


OYR#OYo,OY«,OYV,nY»»,OTCP»OYUG,OYQP,OYr>V,OYOw 

-.‘;3b/iE-03  0.  . 1517F-02  -.27721-02 

0. 

.23650*3 

0, 

.2343003 

-.9062003 

0, 

02P#O2f3,OZH#O2vfn2«fO20P,GZJf-i»D2OH,l)2OV,O2D« 
J.  ltjj^f-02  0.  0. 

-,2V8uO0<J 

0 

-.4103003 

0, 

0. 

1310E-02 

,j<H,3M,GKW#0'<v,.-)*-',GADR»DKO(;.&KOW,i)KOV,OKO* 
il  0 « 2 e:  - 0 !3  0.  • 12962-0 > -.2404003 

0. 

.21  73005 

0, 

.2410004 

-. /779E-04 

0, 

JMP#QMn,  GHH,l)MV  ^n^^fOMOP/nMOll.UMORfOMOVfOMOx 

).  -./n37f*-G3  o.  a. 

-.3235  0 03 
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0. 

0. 
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0. 
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i-S^TC**  OF  PWLSbi'wr  at  cr.NTCR  OF  GRaviTYs  22.35 


sec 

01 

OF  I 

0F2I 

OF  3 1 

OCI 

OC?I 

DC  5 I 

L'CF  I 

0CF2I 

UC2FI 

6 56  I 

i 

0 • 

0. 

0. 

0. 

0 . 

0 . 

0. 

0. 

0 . 

0 . 

0. 

j 

, 130002 

-.895003 

.449003 

-.22900  3 

•1662-03 

. 1532-04 

. 1 422-05 

-.8252-04 

.4152-04 

-. 7 6 2 1 - 0 5 

.2092-0 

3 

,340002 

- • 152  002 

.70  Li; -0  3 

-.3312 -0 $ 

.3152-03 

•2932—94 

.272005 

-.  1 4 1 1,  — 0 3 

• 64  *E -0  4 

- . 1312-04 

.4 1 62-0 

4 

. H 1 O 0 2 

-.19  3002 

.321003 

-.36HL-93 

.44*0-  3 

| Pt-ou 

.39000  5 

-.17*2-05 

• / 6 jt -04 

-.16/2-04 

. 0 2 or  - " 

t; 

.*02002 

-•21500? 

.360 E— 03 

-.377003 

• 5642-0 3 

•5282-04 

.4*42-05 

-.2012-03 

.8052-04 

- , l 6/2-04 

, / 9 7 F - 0 

6 

. » 0«E-02 

-.223002 

.8  75003 

-. 3732-03 

•6622-0 $ 

.622E-04 

,5842-0*. 

- .2  08003 

.31 52-04 

- . 1 Out -0  4 

.9  Sufe‘-0 

7 

.7372-02 

-.221  00? 

.3  7700  3 

- . 3732-03 

. 7 42003 

. o99t-na 

.659005 

-.2062-05 

.3 l 0 c — 0 4 

- . 1 9 2C. -04 

. 1 022-0 

* 

. o5  1 Ov)2 

• .2 1 3092 

.8*1/00  3 

-.37 62-03 

•5032-03 

. 7592-04 

.71 72-05 

- . 1962-03 

.6252-04 

- . 1 e3c-0  4 

. 1072-* 

9 

.3952-92 

-.20  3002 

.9  1 '4003 

-.3702-03 

• 6470  *'3 

•O0 1 2-0 4 

. 7591 -05 

- . 1 3 6 1 - 0 5 

.3522-04 

-.1752-04 

. 1 1 02-0 

l 0 

.920002 

19*»00? 

.9<i  1 003 

— . 30000  5 

.e/lC-03 

. **262-04 

• 7h  3r -0  5 

- . 1 du  2-0  3 

.37*2-04 

- . 1 6 1 1 -0 4 

. 1 122-'' 

tl 

.9252-02 

-.  19200? 

.95000  3 

-.357003 

.»  7 0 C - 0 3 

• 8 31 E-04 

. 7 8 J 2 - 0 5 

-. 1 7«2-03 

.3382-04 

- . 1 *5c -0  4 

.1  122-- 

l 

.9202-02 

19200? 

.953F-03 

- . 3 5 5 L - 0 3 

.6782-03 

.6322-04 

. 7902-05 

-.  1 7 72.-05 

.3912-04 

- . 16*2-04 

.1122-0 

GEUmCToICAL  ^AKUOLCS 


MfJLL  ‘ Opt*.) 

Gt 

31 

SI 

10HAF 

2 - 0 0 0 

.4>»22t-ol 

.3  502E -1)  3 

• 75U2E-0 3 

. 2 /|  u 4 1 -0  1 

1 A ;)  0 0 

•*045l -0 l 

.266 0 fc - 9 3 

• 64 762-4  $ 

.21212-4 1 

0*000 

.3'5  l 4t  -0  1 

.20  3 '00  3 

.54592-0  7. 

. 1 797^-0 1 

-l * >00 

.2**01  -0  1 

.12/92-0$ 

.40092-05 

.13:  >e-ot 

-2*000 

• 1 O'Sot  -0  1 

.61042 -00 

.30 18E-93 

. 10C7L-0 1 

\ 


82 


APPENDIX  B 


COMPUTER  LISTING 
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PROGRAM  5ES"AVC ( INPUT /OUT RUT, T APE5  = INRUT, 7 APC6=0U THU  T, TAPE 1 
*, TAPL2) 

OlMfNSlON  YC12),YP(12) 

III ME  NS  X UN  RP(70&),VP(700)»HHXn(700),XP(700),YYP(700),ACC(700) 

*,BETARt)(7  0y ) ,UP(700)  , THTD(700),ZP(70u),hP(700),wD(700) 

*,PSXtf(  70  0 ,wAV(700),PCP(700) 

COMMON  /A/  PUUT , QUOT  , RDU  I »PHfDOT, T HTUU I , PS  I OU I * UOOT , VDOT , WOO T , 

*XDUT,YL)OI  ,7UUT 

COMMON  /B/  P,TJ,R,  X , YY  , Z ,U,  V,  W,PHI  , TH  r ,PSI 

COMMON  /NPp/  DYPf DYQ,i)YR,DYV,DYW,DYOP,l)YOU,OYDR, DYl)V,DYPW, 

* nZP,OZQ,OZR,DZV  ,0Z*rn£l)P,0Z0Q,0Zl>n,0£UVrDZDH, 

* ttKP,l)KQ.DKR,l)KV,OKM,DKOP,DKOU,DKDI<,OKDV,DKDW, 

A f>MM,  OMO,  OMR,  OMV,DMYi,I)m|)P,omDQ,  DOOR,  pMOV,  DMO", 

* ni|P,ONQ,ONR,DNV,DNw,r)NOP,r}NOR,ONDI<#ONl)V/ONOW 
COMMON  /IJErY/  NK,  DELTA,  F X , F Y , F K , FN , XUOELU, OK AGY , OR AGN  , 

*nELTAY,OtLTAN,DMACK,BETAR»ur-  TH , PEL T A X , OEL  T AK  , Df-  THI,  T HR ATE 
*,l)ELP,DCLS,KPM,  IFUIL,  IFIN 

COMMON  /In/  AA, \IX,AIZ,AM,BB,CB,CF,Orfi,OXI)U,FO,G,NST,NVAL, 

*Pl ,RHG,SP,  ill),  WL,  x LG,  XHi,  COLL,  CONN,  KRUUOE , CC , UU , AMI,  ALOl),CLD 
*,NC,NG,SHTiiwN,  I plot,  ipt,  aiy 
COMMON  /PRnT/  0 I ,NSTCP»NPRNT,  IP 

COMMON  / IEmP/SX,SY,SK,SN, NAVE0G,AEHGDG,HYDROF,SPHYDG,SEALDG, 

*$K  I HOG»  FI  N()  G 

COMMON  /TEMPI/  THXGH, TLOw.SHJPDG,  TGTLDG.TX 
COMMON  /ThrST/  rci)Nl,TC0N2,TC0N3,I0RAG,CC0,THT0 
COMMON  /U/  Gl(C'>),SI(2S),Sl(25),PHlH20),TDRAF(25) 

COMMON/ x/  l SELT (25) ,DX (25) ,I)F  I (25) ,t)F2I (25) ,DF3I (25) , DC  I (25), 

*0C2I (25) ,DCSI (CS) ,OCF U25),l)CF2I  (25)  ,DC2FJ C 25 ) , B3HI (25) , XS* (25) 

COMMON  /!'/(,(./  X ARM,  Z ARM,  8 ACE,  YARM(R)  , DEL  JET  (4)  ,RMCP(a)  ,NJfc  f 
* , ALPHA (fl) 

COMMON  /HGT/  BOUYAN, TNHGT , HMq, nXO 
COMMON  /FLO*/  PC , OF , DO, VOOf P , AUP , A ( 

COMMON  /VOlM/  VULP 

COMMON  / T U v * / PWH,PWZ,PMC,PZC,PSLZ,PSLM 

COMMON  /"Ay/  OuFT ,RAMP,WPER,CCL,CAY, I BUG, F (25) .BETA, IH.wOEP.OFFStT 

* , "LG , ICn,Xn,WlJ,FTAU,TU 
COMMON  /ASnFC/  CCX 

SlNH(E)=(ExP(E)-CXP(-E) )/2. 

SECH ( XKGJSM,  / O'XP ( ARG) +EXPC-ARG) ) 

oA I A COO, CO  1 »Lo2,CU3»CU4/ U.53R2 4 656,-52.707 16255, 1 07. 18 76202, 

*- l OO.Ru5obi B, 55.25071574/ 

L 

C CALL  INPUT 

C 

(_TAPC  = 1 

CALL  INP((LTAPE) 

c 

C INITIALIZATION 

L 

|-|FTHI  = 0.U 
OELTAX=0.0 
DcLTAYpO.H 
OFLTA/=0. 

DEL  T *K  = 0.<» 

DtLTAMsO. 

nCLtAN=0.lt 

SHIPDG=0. 

TnTLDG=0. 

T = 0 . 0 

» >1 


U=uu 

jnox=o 

COEF=0.5*RhO*W|_**<?*SP**2 

nE=12 

yU  )=U 

Y(2)=V 

y(3)=w 

y(9)=P 
v(ii)=Q 
y(6)=R 
y(7)=X 
y(8)=YY 
Y (9)=Z 
y(tO)=PHl 
Y<n)  = fHr 
Y ( 1 2) =PSI 
C 

L CALCULATE  wflGHt  OF  CRAFT  AT  INITIAL  TRIM  wITH  NO  NAVE 

C 

jNWGTsO 

.•ALL  S L A N A V C rt  X , N Y , NZ  / nK  , NM , wN , VOL  / AQ , Y , T ) 

rUUYAN=-«2 

WM0=WM 

WXO=NX 

TNWGT=1 

call  NUNGS(T,nT,NE,Y,YP,INDX) 

KNT  = 1 

tp=t**l/sp 

mPIKNT)=U*sP/I.G89 

vPIXNT)=-AtAn(v/U)/DTR 

wPCKN  n=«*sR 

WlHKNTJsHDOt 

pp=p 

nP  = QASP/NL/DTf? 

pP  UNT  )=K*sP/NL/DTR 

yPlKNT)=X 

YYPCKN()=Yy 

^PIKNT)=A*wL 

pMIO(KNTJ=PHl/OTR 

THTOCXNT;  = THf/|)TR 

pSIP(KNTJrpS[/l)TR 

WAV(KNT)=0. 

*CC  CKN I )=F  RUunL'A*.‘,N(K*UPVDOT  ) 
pAUIUS  = O.I> 

jFIH.NE.U.j)  RaOIIJS  = U/R»NL 
PFTAPO(KNT)=l3n  AR/OTR 
TXP=TX*CUfcF 

turg=totldg*coef 

L 

c integration  uy  rungs 

c 

10n0  inc=o 

uAMSAM*G/22-'IO.*(0.5*RMO*wL**:S) 

p8t-i  = 2.*Po 

pTNTl)=  IhTo/OTR 

pSP  = SP/ 1 .089 

rPETApoEI A/n I R 

CRAFTL=X5w(N5T)*WL 

wNITLCO.Ji'.jJ  Pam,  CRAFTL/PRO,PTHTU,LC,PSP 
3o0  FORMAT ( 1 M 1 ,2  IHCRAFT  CH AR AC t CR I S T I CS /5X , 7H<* f L GH I = , F6 . 0 , 5H  TONS/ 
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Ai;X,  1 3HCRAFT  LENGTHS  ,Fo.l,i»H  FT  . /S  X , 1 UHCUSH  I ON  h I 0 THs,  F6  . I , 4H  FT. 

INITIAL  THIM=,F6.2,5H  DFG  . /SX , t 4H  I N i T I AL  Oft  AF  T = , F 7 . 3 , 4 H FT. 

*/SX,  lanlNUIAL  SPEED=»F5.0,6H  KNOTS) 
lFUh.LT.  3)  hftl  TC(6,  SOI)  PHETA  , .sAMP,  «Ptft  , > LG.CLL  ,CAY 
3 o 1 fCHMAT1//2[M  h a V C CHA ft ACTLftISTICS/SX, BUM CAOINCs, f 5. 0,5h  DEG./ 
*GX,7HhLI6mt=»FS.i,4H  F T ./SX  , 7HPCII  100=, F5 . 1 , SH  Sl.'C. 

*/GX,7HLCNGTHS,ro.l,/|H  FT./SX,9hCLLt«XTY=,F6.1,7H  FT/SEC 
*/SX,Uh-.AX  nU.=,F7.4) 

tFUh.LO.3)  htll  TCC6,302)  PbElA,*OFP,RO,ETAU,xn,T0 

3n2  FPMMA I 1//2IH  KAV[f  CHAftACTEHlsriGS/f,X,eHHCA01NGs,i'5.0,SH  DEG. 

*/SX,5HhDLP  = ,F6.  J ,4H  FT./5X,3hRO=,F7.I,4H  FT. 

*/5X,5HLTAO=,T6.1,'IH  F I ./5X,3HX0=,F8. 1 ,4H  FT. 

*/SX,3HIO=,F&.l,5H  SEC.) 
wRITLlG,22o) 

2;>0  FORMAT  Uftl,5X,lHT,6X,lMU,3X,ilHbCTA,6X,lh«,6X,lF‘X,6X,lHY,7X,lH/ 

*,CX,3HPfH,-3X,5hThETA,5X,3HPSI,OX,2FIPL,6X,JHGF,6X,2HQO,2X,6HVOLOOT, 
ArX,2H«0,‘jX,2HhFI,7X,3HV0L)  C 

WHITL10,222) 

2^2  FORMAT ( 3X, rhSEcS.PX , 5HK  NOT  S , 9 X , JHD t G , l X, oHF  T/SLC  , a X , 3H/LC , 4 X , 3H/LC 
* , 6 X , 2 H F Tr5x#iH0FG»5X,iFlDtG,5X#3HDEG,5X,3HPSF, 1 X , 7HF T 3/SEC , I X , 

*7HF  T3/SCC, lX,7HFT3/SEC,aX,3HFT2,5X,2HF T , 7 X , 3 HE  I 3/ ) 

IFUP.NE.O)  hniTC(fa,221) 

2t I FORMAT  U7X,2HSX,9X,?HSY,9X,2HSN,5X,6HXUDELU,6X,5HOHAGY, 

*6X,5HDKAGN,5X,6HOELTAY,5X,GHDELTAN,6X,1jHTHIGH,7X,«HTLOW,5X, 
a GHTUTLDG//) 

WftITL(0,2i)o)  TP, UP (KNI ) ,VP(KNT) ,nP (KNT ) ,XP(KNT) , YYP(KNT), ZP(KNT) , 

ApHiDIKNT ) ,TMTD(KNT) ,PS1D(KNT) , PC , OF ,GO, V l)U  TP  , AOP,lsAV(KNT ) 

TFUP.nE.H)  "RITE  (6,206)  SX  , SY  , SN,  XUDELU  , DR  AGY  , DHAGN , DLL  TA  Y ,DELT  AN 
A , In IGH, TLUw, TOTLDG 

*,  WAVLDG,  ALHC1UG,HYDK’IJF,SPRYDG,S£  A LOG,  SK  INDG,  F TNUG 
tNTHY=0 

rcx=cc 

nil  2 I=1,nsTLP 

iFUNTftY.EQ.l.AMD.  tCU.NE.O)  CALL  SELDER  ( LT  Afjfc  ) 

j N I ft  Y so 
T NC= INC  F 1 

fALL  ftUNGS(T,DT,NE,Y,YP,INDX) 
ll  = Y(l) 

V=Y(2) 

X = Y f 3 ) 
p=Y(lt) 

0 = Y f S) 
p = Y ( 6 ) 

X=Y(7) 

VY=Y(8) 

7=Y ( 9 ) 
pHl  = Y (10) 

THTsY  Cl  l) 
pSl=YU2) 

TFUNL.NE.nF'KNT)  go  tu  2 

KNT  = KN  I ♦ 1 

TP=T*"L/SP 

ijPIKNT  )su*sP/l.bB9 

vPIKNMs*atAN(v/U)/OTR 

mPIF.N1 J s s * g P 

wD l KN  T ) = «T)f)T *SP*  aP/nL/G 

PCP(KM)=Pc 

pP=P 

nP=G*Sft/"| /DTR 
pPIKNt ) sft*gH/"L/D TK 


.3^. ...  . .. 


r 


XP(KNT)=X 
YYP(KNI )=YY 

7PIKNT)=A*wL 
PMIDUNT  J=pmI/0TR 
THTOCXNT  J=tHT/oTR 
pSIO(KNTJ=pSl/OfW 
CTsCLL*TP 

llT=CY  (7J*C0SlKCTA)  + Y(a)*SiNf BUT*)  )*^L 
r,0  TO  (H',ll,l2),lw 
[0  WAV(KN  r)=-(wAMP*SIN(CAY*(UT-CT) ) ) 
rD  TO  13 

ii  tt=abscu r-c o/weg 

A1=CAY»(UT-C1+0FFSET+IT**LG) 
wAV(KNT)  = ,iaMP*sECHUI)**2 
r 0 TO  13 

h=roep 

jw=TPtTO 

S=UTtXO 

pF=S/TW/SORT(G*H) 

TFIRF.GE..3)  GO  TO  fl 
«K  = J 
rO  TO  5 
'I  pF2=RF*RF 
bF3=KF2*RF 
pF«=MFj*RF 

HK=COR*»Fa+C03*RF34C02*RF2+C01*RF+C00 
5 rAY=MK/H 

nMEGA=CAY*sURT(G*TANH(HK)/CAY) 

rFL=OHLGA/tAY 

rT=CLL*TW 

mk2=2.*hk 

<;HK2sS1NH  (hK2) 

ARG=HK2/SH<2 

aRGJ=1.+Ahr 

aRG2  = -ARGI./(ARG*(  t.-HK2/TANH(HK2)  ) +0 . 5 * ARG 1 **2- AKG 1) 

P0K=CAY*KO 

rALL  hcssel(3,rok»bj:4J 
4WGU=CAY* ( xU-CTYUT) 

wAV(KM)  = (|tTA0*R0/S)*SGHT(ARG2)*BJ3*C0S(APG4) 
l 3 CONTINUE 

ACC(KNT)=FRlHJPt**2*(R*U+V0UT) 
pAUI US=0 .0 

p(:  TARD(KNT)suf.  tar/dtr 
iFIR.NL.G.'j)  R AO  I USsO/R*  WL 
TXP=fX*CUCF 
XORGsTUTI.OrxCOCF 
pH  T =PWM+Phr  + PSI.M 
bG=PC-21  1 7. 

WRITE  (0,200)  TP,UP(KNT),VP(KNn»l->P(KtNn,XP(KNr),YYP(KNT),ZP(KNT), 

»PHID(KNT  ) ,THlO(XN  r ) ,l‘3IO(KNT)  ,PG,')F,  UU,  VDOTP,  wl)(KN  I ) , »AV(KN I ) , AOP  C 

jFUP.YE.O)  wRl  rE(0,206)  SX,SY,SN,XUI)tLU,ORAGY,ORAGN,DCLlAY,DELTAN 
* , THIGH,  n ow,  TOTLOG 

*,WAVLOG, AFRUt)G,HYOROF , SPR YOG , SE ALDG , SK 1 MDG , F I NOG 
TFUCO.Em.o)  GO  TO  15 
7TESr=ZP(XMl)+wAV(KNT) 

7TST  = ABS(ZTTS  f +CC-CCX) 

XFlZT8l.LT.Itu)  GO  TO  15 
fMIRYSl 

rrx=tc  yztest 

l 5 CONTINUE 


I 
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I N C = •) 

Z fONTXNUE 
pCPCl  )=0. 
r)0  803  I J K = 2 » K N T 

aa3  pCP(IJa)=(pCP<iJk)-21 17,)/loo, 

kNBsO 

oil  707  IKQ=l,KNT,0 
«n8=kN8+ 1 
ZPIKN«)=ZP(IKQ) 

ThTO(KN87  = ThTO(XK(3) 
pHI0(KN8)=pHl[)( IKQ) 
pSI0(K.%8)=pSir>(  IKQ) 
wAV(HN8)=waV(.IkQ) 

767  PCP(K^8)=PcPUKU) 
jptskno 

lFUPLUT.OT.2)  CALI.  EXIT 

TFUPLUT.CQ.O)  CALL  PLUTr<ZP,THTU,"AV,PHID,PSID,PCP,IPT,NC»NG) 
rFUPLLM.EO.l)  CALL  PLOT X Y ( XP , Y YP , KM T ) 

tFIIPLOT.Lq.2)  CALL  PLtITT ( ZP , T HTD , « A V , Pm  Ip , PS  I U , PCP , IPT , NC , NG ) 
f F t IPLuT  • t.0.2 ) CALL  PLO  TXY  C XP , Y YP , KN  f ) 

2nO  pOKMAT C6F 7.2» 0F8.3, OF  8.0 ,F7.3,F7.2, FI  0.0) 

208  fO«mAT18X,11L'U.3//) 

STOP 

pNO 
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L 

L nEFIMIlIUN  OF  INPUT  FLAGS 

L 

L nSTEP=NO.  of  Time  steps  to  execute 

C „,PUN  SPRINTING  INCREMENT 

• c tP.nl. 0-  ° t? I N T DEBUG 

L if-  IN.NE. D-iNCLUDE  PIN 

p iPLUf— — — -FLAG  ON  PLOTTING 

p IT  I PLUT  ~0  , PLOT  PL.UTT 

l ir  iplot  =i  plot  plotxy 

L IT  I PLOT  = 2 PLOT  PLOTT  AND  PLOTXY 

P IT  IPLOT  GT  2 PONT  PLU’ 

L tPT— — — —NUMBER  OF  STEPS  TO  PLOT  PLUTT 

L uJET— — — --NUMBER  OF  JETS  FOR  THRUST  VECTOR  CONTROL 

l rNT— —.—PRINT  FLAG  TOR  CUMULATIVE  INTEGRALS  AND  GEOMETRICAL  VARIABLES 

. IT  TNT=U,  HUNT  PRINT 

P IE  INT.NL.O,  PR/NT 

p fM------- — -flag  for  rave  type 

P IT  t R - 1 , SINUSOIDAL  RAVL 

P IF  IR  = ?,  SOL  I TARY  RAVE 

L IT  IwsjS,  EXPLOSION  rave 

p i CO— FLAG  FOR  GENERATING  NEW  DERIVATIVES  WHEN  uRAF'T 


SUBROUTINE  INPT(L) 

DIMENSION  T I I L I.  ( 2 0 ) 

COMMON  /A/  PDOTf QDOT, ROOT, PHJ DOT. THTDOl , PS  I DOT , UDOT , VDUT , WOUT , 
*XDOT,YDOI,ZDUT 

COMMON  /b/  P,u,H,X, YY»Z.U»V,h,PHI,TMf,PSI 
COMMON  /DERV/  nk.deli  A,FX,FY,FK,Fn,XUDELU,ORAGY,DRAGN, 
i *oELTAY,DtLTAN,ORACK,BETAR,r>(-  TH » DLL  1 A X , DEL  T AK  , Of-  TH  I , THRA  Tt 

*,DfcLP,DELS,HPM, IFOIL, IF  IN 

COMMON  /V COFF / FYNCL/F INYV,FjNYR,FINKV,FlNKR,FiNNV,F'INNH 
COMMON  /F  flyL  / C»ALFA»GAMA,XF 

COMMON  /GtnMM/NSw(25) , Rl (25,25) ,R2( 25, 25) ,01 (25,25) 

COMMON  /IN/  AA, AIX,AU,AM,BH,CB,CF,DfR,DXOU,FU,G,NST,NVAL, 
*PI,RHU,SP,U0, RL »XLG,XEG, CULL, CONN, FRUUDE , CC » OD, ANU, ALUO, CLD 
*,NC,NG,SPTiiRN,  IPLOT ,IPT,AIY 
COMMON  /PRnT/  DT,NSTLP,NPRnT, IP 

COMMON  /PReS/  C0IS,RHUWA ,PHIO,PHI 1 , a TM,PMAX, AC,DEM, IPR 

common  /temp/  sx,sy,sk,sn 

COMMON  /IFMPI/  Thigh, TLOW,SHlPDG,TOTLDG 

COMMON  /IHRST/  TCON1, TC0N2, TCON3, IORAG, CCO, THTO 

COMMON  /U/  Gt (25) ,SI(25) ,Sl(25) ,PHO(25)  , TORAF  (25) 

COMMON/ X/  JSECTC25), 01(25) ,DFI (25) ,0F21 (25), OF  3 1 (25) , DC  I (25) , 
*oC2K2S),Dc21(20),DCFI(25),DCF2I(25),0C2F'I(25),f)3BI(25),XSR(2S) 
COMMON  / 1 N r>  IK  / CR,CT,S, OMEGA 

COMMON  /1VCL/  XARM,ZARM,BACE,YAHM(R),DLLJET(4),RMCP(a),NJET 
*, ALPHA(a) 

COMMUN  /ABt/  Draft  (25)  , RE  I GHT  , UUBB , BU13L  , W AL.B , SLBOW , SLS  TRN  , THE  T A , 
*I5FPTH,SPRAyL 

COMMON  /CDE/  DRISE(23) ,ENTRCE(23) ,CHINE(22) ,HSPRAY(23) 
common  /SOS/  hs'r  (25) , OEL 1 , UEL2 » N l , N2 

COMMON  /NOD/  DYP,DYD,OYR,DYV,DYR,OYDP,DYDG,DYO(<,OYOV,OYDW, 

* nZP,DZa,DZR,DZV,DZW,DZDP,DZDU,nZDR,DZDV,OZDW, 

* PKP,DKQ,OKR,DKVfDKW,UKI)P,DKDO,DKDH,DKOV,OKOfc, 

* DmP,DMD,OMR,DMV,DMw,DMDP,DMDD,D«DR,DMDV,DMDW, 

* DNP , DNQ , UNR , DNV , DNW , DNDP , DNDU , PNDR , OND V , DNDW 

COMMUN  /RAV/  OWFT , RAMP, RPtW,CLL,CAY, IHUG,F (25) , bFTA, I W, WOE P, OF  F SE  T 
R,wLG,ICO,X|j,Kfj,F.  TA(),TU 
COMMON  /PSeAL/  (HTO,ThTS 

common  /asofg/  CCX 
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changes  more  than  icn  feet 
ir  ico=o,  don r change  derivatives 

IF  ICU.NE.O,  CHANGE  IN  DRAFT  RFOUIRFD  TO  CHANGE  DERIVATIVES 
jFU.EQ.2)  GU  to  1000 
BEAD  AND  KRllf  INPUT 
REACH5»IU5)  IT  I TLC ( I ) * I = l » 20 ) 

READ (5,  10  l ) OT,N$TCP,NPRNT , IP , IF  IN , I PEG T , I PT , N JET , INT , I BUG, I W , I PR 
*,ICO 

rEAOCS, 100)  THT»PHI,PSI,Z 

re  AO (5, IOC)  AA,BR,CC,00, AM.DXOU, AIX, AIZ, AIY 

READ  (‘J  / l Oil)  wL,SP,RHC,ANU.C01L,CDNN 

READ (5,  IOC)  UMlGA,CR,CT,S 

r E A D ( 5 # 100)  C CO/ THTO/3PTUHN/DFTH 

BE AO ( 5 , 100)  XARM,ZARM,RACt 

BE AO (5#  100) (YARM( I) , l=l,NJET) 

RE  AO (5 , 100)  lOlLJCTCl) , I=l,NJET) 

BEAD  15, 100)  IRMCPI I ) , I = 1,NJET) 

BEAD (5/ 100)  1 ALPHA (I ), 1 = 1, N JET) 

WHITE (0, 20y ) 

WRITE (0,20 i ) (TITLE(t) ,1=1 ,20) 

wHITE(O,2i'2)0T,NSTEP,hPRNT,IP, IF  IN, I PLOT, IPT , N JET , I NT , I BUG,  I W , I PR 
*,1CU 

wRITL(0,2fl3)  THT,PHI,PSI,Z 

wRITL(0,2iV|) AA,BQ,CG ,00, AM  , 0XDU , A I X , A I Z , AIY 
wRITl(o,2os)«l,sp,rho, anu,cpll,cdnn 
wR  I TC ( 0 » C07 ) UhFga,CR,CT,S 
wRITE(0,20b)CC(),THTO,SPTUHN,OFTH 
wRITL(0,2lo)  XAHM, ZARH,HACE 
WHITEC0,2i 1) 

wHITlCO,  1 o.j)  CYARM(I)  , 1 = 1 ,NJCT) 
wRITE(0,2(2) 

wRITE(0,10<j)  (Of'LJL'T(I),l  = l,NJET) 
wHITLCO, 21 3) 

WHITLCO,  lily)  (RMCP(I)  , I = 1,NJET) 

WHITlCO,  2l/|) 

WRIT  ECO,  10 IJ)  (ALPHA  (T)  ,1  = 1,NJET) 
read  and  wr  l TC  INPUT  FUR  wave 

read (5, 100)  OWE  I , WAMP, wPFR,UET  A.wOEP, XU, HO, E TAG, TO 
WRITE(0,20D  OwFT,WAMP,,iPER,BETA,wOEP,XO,HO,ETAQ,Tl) 

r£AO  AND  WRirc  INPUT  FDR  PRESSURE 

READ  (5, 100)  CDIS.RHUWA,  ATK,  PHI0,PHU,  TFITB,  THTS 
wRI  TE  (0,2/|2)CO  lS,liHOWA,  ATM  , PH  10,  PH  II  , TFITH,  THTS 

READ  AND  WRITE  INPUT  FUR  SPRAT 

rEAD(S,102)  NST 

oFAOIS,  100)  HUBL  , RUllB  , W At  U, DEPTH 

BEAD (S,  100)  SLRUw, SLS ! RN,  (HETA 
rF AO (5, tOO)  (ORISCt I) , 1=1 ,NST) 

READ (5, 100)  ILNIRCECI) ,T=1 ,NST) 
rEADCS,  100)  l EH  I NC ( I ) , I = 1 , NST ) 
rEAO(5,10P)  (NSw(l),  1 = 1, MST) 

rEAOCj,  100)  (XSw(I)  , t = t ,'JST) 
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1 


si 

I S2 

S3 


r 


S'» 

L 
L 

u 


i 

t 

L 


I 

t L 

I. 

C 

I 


rFADC5»  100)  IHSWCX),I=1,NST) 

1)051  I = t»NsT 
9VS=N5*  ( I ) 

bead (5.  too)  truu,j),J=i,NVS) 

1)052  X=l»NsT 
N V3  = NS  * ( I ) 

PEA0(5»  10(1)  i*l(I,J),Jsl,NVS) 

i)05j  1 = 1 > ns  I 

0VS=NSRlI) 

READ  (5 » 100)  Cw2lI,J),J=l,NVS) 
w«ITE(0,2lsJ  NST 

wRITCC0,2to)  BUHL »0UB6,» ALB, DEPTH 

wRI  T L (0,2  17)  SLBOi»,SLSTRN,  THETA 
wPITC(0,2lfl) 

wRITCCO, Toy) (DRISECI) »I*1»NST) 

WRITE (0,2 19) 

wRITLCO,10y)  IENTHCC(U»I  = 1»NST) 
wRXTC(0,22.j) 

WRITE (0, lOy) ILHINCCI), I=1,NST) 
WRITL(0,22i) 

WPITCCO,  1 (12)  (NSWC  1) , 1 = 1 ,NST) 
wR I TE  CO, 222) 

WRITE(0,  l(l.j)  (XSw(l)  , 1 = 1 ,N3T) 
wR I Tt ( 0 , 22j) 

WRITE (0,  Jo.j)  CHSw(X)  , 1 = 1 ,NST) 

1)054  1 = 1 ,NgT 
WVS=NSn ( I ) 

WRITE (0,22/1)  a,  coin, J)  ,J=1,NVS)J 
wR I TE ( O , 225)  Cwl (I,J),J=1,NVS) 
wRITCC0,22(,)  Cw2CI,J).J=l,NVS) 
CONTINUE 
CONSTANTS 

nC=20 

N 0 = 0 

0=32.2 

Pl=3. 1415927 
i)TR=PX/  180. 
p = U = R = V = * = x = W = 0, 

Od=l  . 
o=un 

wEIGHT=AM*2240. 

CONVERT  To  radians 

7HT=  f HT  *DTr 
pHl=PHX*i)TR 
pSI=PSX*DTr 
OHEGA=0MLGa*DTR 
ThTO=  THTU*oTR 
pFTAsbL I A*pTR 
XHT0=TmTd*oTH 
7HTS=TmTS*|)TR 

CONVERT 

a I X = A I x*22/|0  . 

A rr  = AlY*22/|0, 

A T 2 = A IZ*22tO . 

»m=am*2240./G 


C 

C 

C 
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I 

i 


I 


SP=SP*1  .ofiq 
SPTURN=SKTtlKN*i  .689 
pROU&E  = Si'/SURT  ( G*  AL ) 

L 

E CALCULATE  CAT,CfL»F  FUR  SUBROUTINE  SWAVE 

L 

Rn  TO  (4I,n2,«3),lW 
<|1  cAF=‘l.  *PI**2/(G*wPCR**2) 
cEL  = D.i.*GAWPLR/PI 
wLG=0.i,*b*wPER**2/PI 
r,0  Til  'IP 

l|2  cAT  = D.86o*gaKT(  WAMP/WOEP ) /WOEP 

cEL  = 0.U*UumT  lG*WDEP)*(2.  + »AMp/viOCP) 

wLG=CEL*''PeP 

rO  TO  '19 

«3  cAV=9.»PI**2/(G*30.0a*2) 
CEL=O.5*G*~0.0/PI 
wLG=i).U*G*30.0**2/PI 
q'l  CONTINUE 

HIT  19  I = 1 , nS  T 
1 9 p(I)=XSWtI)-AA 

L 

L CALCULATE  TIME  increment 

L 

nT  = C. ILF./ AOS  CCfL-SP*COS  (BETA))  )/DT 
f)T  = Df*bP/*L 

L 

C iqON  OIMENSiONALIZE  input 

L 

„ENn^  = o .5*«hu*wL*a5 

i IX=AIX/DENUM 
A I Y = A I Y /UENUM 

aIZ=aiz/uenom 

am=am/(0.G*KHU*wL**3) 

2=1/ «L 
VLG=AA/KL 
XFG=1  .-XLG 
X ARM=X ARM/wL 
ZARM=ZARM/wL 
no  20  1=1 ,nji t 

TO  vAUMCl JsT  ArIHD/WL 
fjO  21  i = 1 » m3  T 
yS" t I )=XSR(I J/WL 

WRITE (0,227)  SP, AM, Atx, AIY,AJZ,FR0UDE 
C 

u CALCULATE  Mi, N2, DEL  1 / 0EL2 

C 

msti=nst-i 
no  5 i=2,nsti 
rSAVE=I 

0EL1=XG.KI)-XSW(I-1) 
nF.Ll  =UELl  *WL 
0EL2=XSwCl+l)-xSw(I) 
r»EL2=DLL2*wL 

tfiaosu .-nrL2/ocu ) .gt.o.d  go  to  o 

0 CONTINUE 
6 M l = I3AVL 

hj2  = N5T*  ISA  vC  t 1 
E 

L CALCULATE  aC,OL'M  FiiR  PRESSURE 
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r r r.  r-  c-  r. 


aC=2.‘0B*XshINi)awl 
i)FN  = i>.l,*HHn*'*L**2*SP**2 

CALCULATE  Fll  AT  CG 

nO  3 I=2/NST 
F01=AhS(AA-XSw(I-l)*WL) 

Fn2  = ABS(A*-XSw  ( I)*««L) 

IFIF02.L1  .HJ1)  KFOsl 
3 CONTINUE 

READ  AND  WRITE  NUN  DIMENSIONAL  OCR  I VA  T I VCS 


10oO  iFIL.EU.2)  GO  TO  1001 
lPDEK=U 

THTO(’R  = THTn/OTM 
wRITC(0,2i>r)  LC<J,THTOPR 

CALL  DLRUA,ii8,CC0,DD,THT0,PHl,NST,Nl,N2,DELl,DCL2,HSw,NSw,XSR, 
*f)l»Wl*'*2/HHO»KL) 

CALL  GLOUa,BB,CCO,OD,wL,NST,THTU,KFO,FO) 

GO  TO  1002 
pN  TRY  SFLOeR 
1 0 0 1 THTPW  = THf/nTH 

jFUNT.EU.u)  [PDCRsI 

CALL  OLRUa»BB,CCX,DD,  THTO,PM ,NST,N1 , N2, DEL  1 » DCL2, HSW, NSW , XSW , 
*01#wI»'*2»Rm0»wl) 

CALL  GLOUA,Be,CCX,OU,WL,NST,  THTO,KFO,FO) 

L 

L CALCULATE  SPRAYL 


1 0 Oc  RllsiU+l 

oo  14  1=1, NST 
fSAV=I-t 

tFlFNTRCt(r),Nc.o.O)  GO  TO  15 

14  CONTINUE 

15  sPRAYLslXSWlNln-XSWl  ISAV))*WL 

SX=0.0 
SY=0 .0 


fiKsO.O 


SN  = 0 . 0 


* 

* 

* 

* 


CALL  FINISX/SY  ,SK,SN,CR,CT,S, OMEGA) 

SFYV=DYV+FlNYV 
jFYRsOYR+FjNYR 
CFKV=OK VtFjNKV 
SFKRsOwR+F  j NKR 
SFNV=DNV+FjNNV 
SFNPsONRtFjNNR 

RC=0Y V*ONR_COYR-AN)*ONV 
<?Ch=5F  YV*SpNK-(3FYR-AM)*SFNV 
TFlIPOLM.Nrt.d)  RETURN 

rFUPOLR.Nt-.y)  WRIT F (0,209)  CCX,THTDPR 
wRIT  L (0, 22s ) 

rR  I TL  (0,229)  OYP,OYG,t)YR,OYV,OY4,DYOM,UYI)Q,UYUR,DYOv,l)YOW, 
PZP,DZG,nzR,czv,n/»,nzop,o/OG,DZOR,DZov,nznw, 
OK?,OKQ,UKR,UKV,n«R,OK!)P,OKD(l,PM)R,UM)V,r>iU)w, 
PMP,0Mi3,  OMR,  O.MV  , DMK,  OUOP,  0 MOW, DMDh , D«DV , ON  0»r 
OMM , UNO , DNR , ONV , ONR , UNDP , UNDO , ONOR , PNDV , GNQN 


rRI TC (0,233) 

rR1TL(  0,23,| ) F INYV,F  1'1Y»,F  IhKV  ,K  tNKH.Fl  NNV,F  INNR 
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! 


wRITL!t>,235) 

wRITLC 0,23t>)  SrYV,SFYR,SrKV,sFXR,SFNV,SFNR 
wRI TU  Co,237)  SC 
wRI TU (0, 233)  SCF 

wRITL  CUMULATIVF  INTCGRALS 

wR  I TL ( 0, 23 [ ) F0 
wRITL(0,232> 

wRI  TCI  0,233)  (I»0I(I)»0FX(I),DF2I(I),DF3I!l),UCI!I),DC2l!I), 
*DC3X(I)»t>CFllt)»DCF2ICI),0C2FlU).b3BI(I),!  = l,NST) 

pIlNVCRT  10  DtGRFCS 

0033  I = t ,5 

33  PHU(1)=PM0( I)/oTH 
wRITC  (6»239) 

wRlTll(o,2/iy)  (PHU(I),GX(I),SI(I),S1CI),TDHAF(I),I=1,5) 
cONVCRr  To  RADIANS 
o063  1=1,5 

63  pHU(i)=HHO(I)*DfR 
FORMATS 

loo  fORMATCSF ly.2) 
l 0 1 fORMAT(F1i>.0,1/|I5) 

I Ot  FOHMATC 1015) 

L 0 3 FORMAT  C20A/4) 

2(lO  FORMATC  lHt,  lOHINPOT  oata  //) 

If)  1 FOHMATllA,20Ai|) 

’02  FI'RMA  r C lx,  3 7 HOT  ,NSTEP,NPRNT  , IP,  IF  IN,  I PLOT  , 1 PT  , N JF.  t , INT  , I BUG  , I W , I 
*pR » ICO  /Fly. 2, 1215) 

2o3  fORMATUX,  i7HTn'f,PHI,PSI,Z  /8F10.2)  , 

cORMA  r ! I y . j 1H  A A . HH  . n: . DO  . AM  . nXIMI  . A I X . A T 7 . A I Y/ 


3 forma  r i u,  i7htmt,  phi  ,psi  ,z  /sfio,;.>)  , 

(I  fORmaTI  lX»3lHAA,rtU,CC,DI),AP,DX0U,AlX,AIZ,  AIY/ 

*/lF  10.3#  Flo.  1,F10.3,3F10.1) 

3 FORMAT!  Ix»23h.>L,SP,kho,Anu,CI)LL,CDNn/3F10.3,F  12.9.2F10.3) 

7 fOKmaTC  IX, (5mOmFGA,CR,CT,S  /MF10.2) 

B FORMAT!  lX,20HCCU,THT0,SPfl'RN,0FTH  /3H0.2) 

9 FORMAT ( in  I ,0M0RAFT  = ,F0.2, 1 0 X , 5H T K I Mr, F 6 .2) 

0 FORMAT!  lx,  l'IHXARM,ZARM,HACi:  /QF10.2) 

X FORMAT!  tX»UHYARM!I)  ) 

1 FORMAT!  I x , 9mdFl JLT ! I ) ) 

3 FORMAT!  1 X , 7HRMCP I I ) ) 

'I  FORMAT!  IX,  10HALPHA!!)  ) 

3 FORMAT! IX,“MNS1  / X 5) 

b FORMAT! IX,  20MB U8L,BUBb,w ALB, DtPfh  /HFlO,2) 

7 F'lRMA r ! IX, 1HMSLU0",SLSTRN, thcta  /ok  10,2) 

8 FORMAT!  IX, flMDRlSCd)  ) 

9 FORMAT!  IX.gHl  nvRCT.!  I)  ) 

0 FORMATl  1X,8HCHI,NC!  X)  ) 

1 FORMAT! tX,3MNSw) 

2 FORMAT! IX, -HXSW) 

3 FORMAT! IX, JMHSw) 

n FORMA T!/I5,2X,2ND1,2X,3F 1 X . 3/ U ) X, 8F 1 1 . 3) ) 

5 pORMA  T 17X,2M«1 ,2X,HF1 1.3/! l!X,8Fl 1 .3)) 

6 FI'RMA  r (7X,2Hn2,2X,8Ft  1.3/U  1X,8FU  .3) ) 

7 FORMA  I !///// IX, IbHCONVLRTEO  INPUT  / 1 X , 2'IHSP , AM, A I X , A I Y , A I Z , FRUUDt 
*/0G)2.'l) 
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cj  rj  rj  cj  rj  rj  rj  rj  rj  rj  cj  rj  rj 


23 

7V 


FORMAT  C/// i X * 27HNQN  DIMENSIONAL  DERIVATIVES  /) 


DYDV, 

DYOW/J 

OE 

12 

,4/ 

OZDV, 

DZDM/I 

OL 

12 

.4/ 

DKOV, 

DKOR/I 

Ofc 

12 

.4/ 

CMD  V , 

OMDW/I 

OE 

12 

,4/ 

DNUV, 

UNOh/1 

Ot 

12 

.4) 

AVITY 

= ,F  1 

0; 

2) 

31  fORMATC  )X,40HcEnTER  UF  PRESSURE  AT  LtMER  UF 

32  FORMAT  ( /ix,3fisr.c,  9X,2riD[,8X,3HOFI,7X,4HDP2i,7X,4MDPiI, 

*flX»3H0CI,7x»RHDC2I»7X»4HDC3I,7X»4HULFI»6X,3rtUCP2I»6X,SHDC2FIi 
*7X,4HB3HI) 

33  FORMAT  CIS,  1 IE)  J .3) 

3 'I  fORMATC/Ix,2jHSTaBILIZER  COEFFICIENTS 

* /7H  Finyv=,L1j.‘I/7h  FINYR=,E13.4/7H  FINKV=,E13.4/ 

*7H  FlNKR  = ,L'l  J.4/7M  FInNV=,E13.4/7h  F INNHs,e13.4) 

35  FORMAT  t /, 1X,33hSHIIj  PLUS  STABILIZER  COEFFICIENTS  ) 

30  fORMATC  bH  SrYV=,El 3.4/6H  SFYK=,E13.4/  UH  SFKV=,E13.4/ 

*6H  SFKR=,El3.4/  6H  SFNVS, L13.4/6H  SFnrs  , C 1 3 .4 ) 

37  fORMATC  /lx/  34HSTABILITY  CRITERION  FOR  SHIP  ONLY  = ,E12.4) 

38  fORMATC  l X » 36HST  ABILITY  CRITERION  FOR  SHIP  PLUS  FIN=  ,E1I 

39  FORMAT c///// IX, 21HGEUMETRICAL  VARIABLES  / , 1 X , 9HR0LL ( DEG)  , 
*inX,2hGI* !<)X,2hSI, 10X,2HS1,7X»SHTORAF) 

FORMAT C IX, f 10. 3,4012.4) 

fORMATC  IX,  3tthn»,FT#  RAMP,  «PER,  BETA,  ROEP,  XU,  iT0,CT~0,  T0/4F  10.2) 

FORMAT ( 1X,34MLDIS,RHG«A, A TM,PM|U, PHI l , T H TB , THT5/I  8 . 2, F 1 0 . 6 , 6 F8.2) 
RETURN 
|£NO 


40 

4 1 

4 2 


J 
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r r 


SUBROUTINE  DER(  AA,BB,CC,DD,THT,PMI,NST,Nl,N2,DEH,DtL2,HSx,NS«(# 
*xSW,DI  »wl  , w2#  KH'J,  *L) 

DIMENSION  HS»(  1 ) , NSN ( l ) f XSM  ( 1 ) 

0IMLNSION  0(25), F<25) 

DIMENSION  r(25),S(25),CSZ(2S),TCMPA(2S) ,TEMPR(25) 

DIMENSION  ol (25, 25), *1(25,25), *2(25,25) 

CChECK=CU-()B*PHl 

PHIU=CC/B8 

p(M )=XS*(n1)**L-AA 

p(Ni  )=CC-HSX(N1  ) -THf*F(NI) 

pO  1 Msl*NgT 

F(M)sXSwlB)*«L-AA 

0(«)=CC-»1SX(M)  -THT*F(M) 

fFlD(M) .LT.O.)  0<M)=0. 

1 CALL  SLC1 (m»D,D1  ,NSN,*1,W2,H,S,CSZ) 
hgi=oo-cc 

COMPUTE  DERIVATIVES  xIIH  PtSPUCT  10  F 

form  integrals 

INTEGRATE  aXIALLY 

CALL  IN7LG(d,n,F,S,CSZ,Nl,N2,DELl,DEL2,NST, IEMPA,TEMPR) 

CALL  NUNOIm(0B,HGT»UHO,WL,B,O,F ,CSZ, TFMRA,TEMPR) 

Rfc  ItJHN 
END 
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SUHKOUI  INE  SFCT(I.()»DUNSW»M»R2»B»S,CSZ) 

DIMENSION  B(2i,J  ,S  (20 ) .CSZ(25)»0I  (25,20),  *1  (25,25),  *2(25,25) 
DIMENSION  pll),NSrt(l) 

FLINLR(X,X2,XI,Y2,Y1)=YH-(X-XI)*CY2-Yl)/(X2-Xn 

R(1)=0.0 

S(I)=0.0 

eszu)=o.<> 

TFMPI=0.l< 

DHAFTsOU) 

jJ=NS*(I) 

Kllstl 

00  1 Js2fJJ 
B02sOKI»Jl 
RDl=Ul (I, J-l) 
rw12=w1 (I. J) 
p» i i =w i c i , j«  n 
rk22=w2(I,J) 
r*21=«2(1,J-1) 

IFIOHAKT.UE.0.0)  CO  TU  4 

IFCOMAFT.GE.OUI.J))  GO  TO  2 

rW12  = FLINf;rCO«aM.H02.R01.RW12,RHU) 

B*22=FLINtRCDf!AFT,R02»RDl/f<w22,HH2l) 

KL1  = 1 
r02=UHAFI 

CALCULATE  A«t A , GIRDER. AND  BEAM 

2  0ELD=HU2-W0l 
w ll)  = R;>12-Rwl  l 

w2D=*w22-Hw2l 

OELS=0.0*OEH)*(HW  (2  + RM  1 + RW22  + RW21  ) 

p ( I ) sRw I 2y  rw22 

S(I)=SCI)+oFLS 

rJM1=«W11+rw21 

CALCULATE  cFNTROIO  FOR  AREA  ABOUT  Y-AXIS 


rn2=lHI)-ho2 

SMOM=(  C02  + .j.'j*DFLUJ*UJMl»0ELD  + 
*(TD2M)ELU/3.J*o.5*DELI)*(HilUtR2U) 
TEMI’l=(tMPl+SMOM 
jFlKLl.EU.l)  GO  T 0 3 

1 continue 

3 rSZCi  )=TLPPl/S(I) 

4 RETURN 
end 


SUBROUTINE  INT[;G(8»0»F,S»C5Z»Nl,N2,0ELl»DUL2/N5T>TEMPAf  TEMPR) 

Ol MENS10N  TrMPA(ib),TUMPR(<i5)»fi2(,-:S)rB<.’K_'5),02F2(2b),82CSZC2S) 

*B2FCSZC2i>),B2IHlF(;,,j),rt2FDl>M2S)#U2(25),t,2F  ( 25 ) , D2F  2 (2b ) , 
*p2CSZ(2SJ  ,02CSZ2(P5),D2DD(-  125 ) , 02FCSZ C 2b) , OCSZDF ( 25 ) , P2FD0F (2S) 
*,D2CZUF 120) 

COMMON  /1NTF0L/  B2I,H2FI,fc)2F2I,B2CSZI,UFC3Zl»B2DDFI,BFDDFI, 
*o21#02FI,D2CSZj,DCSZ21,O2r,DFI»UFCSZI#0FDDH,U2l2I,DC/DFI 
DIMENSION  eUJ,0U),F(l),5(l),CSZ(i) 

COMPUTE  DERIVATIVES  OF  D AND  CSZ  HI1M  RESPECT  TO  F 

NlUhl-l 

OCSZOFl  t)  = <CSZ(2)-CSZ(l>)/0ELl 
0CSZUFlNl)=lCSZ(Nl)-C3Z(Nll))/DELl 
do  i I=2fNlI 

1 DCSzoFa)=.j.'j*ccszu+i)-csz(i-m/oEu 
N2I=Nl+l 

m22=nS I -I 

nCSZDF  CNST)  = USZ(NST)«CSZ(N22)  J/0EL2 
DO  2 I = N21 ,N22 

2 dCSZDFU)=u.0*(CSZ(I  + U-CSZU-1  ))/0EL2 

COMPUTC  ANo  STQHC  VARIABLES  FOR  AXIAL  I NTEGRA 1 ION, 

DO  3 I=t»*sT 

iFlBCI).LG.O.O.OR.D(I) .EO.O.O)  TCMPR(I)=I.O 

ifibu).eq.o.o.ur.d(I).eq.o.o)  go  to  o 

T£MPH(I)=ScI)/0(I)/D(I) 

'I  tFMPA(I)=2.4*Tl-MPRU)*0.4 
R2U)=B(i)*BlI)*TEMPRtI) 
a2F(l)=82(r)AF(n 

r2F2(I)=B2fC1)*F(I) 
r2CSZ C I) SB2 ( I)aCSZ(I) 

P2FCoZ(I)=h2F Ci)aCSZ(I) 

R2DDFU)=P2(I)*DCSZDFm 

B2F0DFCI)=H20Drt  I)*F(t) 

o2U)=0(I)*DU)*TCMPA(I) 

n2F ( l ) =02 ( I ) *F ( I ) 

n2F 2 ( 1 J =U2f  C l ) *F ( I ) 

n2CSZ(I)=D2(I)*tSZm 

n2CSZ2U)=o2CSZ(I)*CSZ(I) 

1)2001  C I)=DL’l  l)ADCSZDF(I) 
d2FCSZ(U=d2FCI)*CSZ(U 
n2FDDFm=n2Dnrm*FC  I) 

3 i)2CZOF(IJ=D2U)*CSZm*OCSZOF(I) 

PERFORM  AXIAL  INTEGRA  T ION 

Ol=0. 

r)2=0. 

Qj=0. 

0«  = 0. 

o5  = 0. 

156  = 0.0 

07=0.0 

oe=o.o 
(30  = 0.0 
olo=o.o 

Q l 1 =0  • 0 
0 12=0.0 


F 


t 

II 


fe'  - 


SUBROUTINE  NOUOlM(Bfl,HGTf  HHil , ML , B # D , F , CSZ  , TEMP  A , TCMPR) 
RtAL 


M.P  »K10  ,*1* 
MlUQrMlO*. 

»KlUP»KiOU,MOP,MlP 

» M I Q 

r M 1 W 

» M I UP  t 

K2P  #K2U  »K2« 
N2UQ,N2D'<. 

.K20P»K200,K2DW,N2P 

,N20 

,N2W 

,N2DP, 

KJP  » K3R  ,KjV 
N3UR*N3DV, 

,K3DP»K3DR,K3DV,N3P 

t N3R 

,N3V 

»N3DP» 

* 

K /IP  ,K/|R  ,K/|V 
M/|  DR  » M/IO  V , 

,K40P»KI|DH,K/IDV,M4P 

,M4K 

,M4V 

,M/|DP / 

i L3 


-L'l 


rL5 


DIMENSION  0a),O(l),FCn#CSZ(l)»  tE'iPA(l),TEMHRCl) 

COMMON  /INTEGL/  B,',X,a2Fr,B2F2I,fl2CSZI,0FCSZIfrt2DDFI,0FDDFI, 
*n2I,DC!FT»n2<.&ZI,OCSZ2I,r>znOF  T*OFCSZI#nFDOF  I , D2F2 1 . DCZDF  I 
common  /Non/  fiyP,or(},i)YR,orv,Dr^.oYDP,(jYuo,.TYOH/Oniv,orot*, 

* OZP,OZQ,OZK,DZVf  OZ'WDZDP,DZDU,DZOR,OZDV,DZD*«, 

* rKP,DKQ,D*P,  UKV,DKM,DKDP,DKD0,1»'D»,DK|;V,DK0M/ 

* nMPjDMQjOMR.OMV.DMWfOMOPjO^OGfOMOK/OMOV^MOw, 

* r>NP,DNQ,DNP,DNV,DM»»DNOP,DNCQ,DNDH,DNDV,DNDM 
Pl  = 3.1'lli><>27 

4KY=1. 

»KZ=1. 

H=  0.2S*PI*KHO 
Cl=  H 

cz=  C1*AKY 
C3=  M 

ca=  C3*AKZ 

1.2  = 0.5*kho*WL*«2 

L3=  L2**L 

L«=  L3*v»L 

L5=  L4*wL 

710*  = -82l*C l 

Z lDP=BB*ZtpW 

Z 10  = (U( 1 )*>2*F (i) )*C 1*TEMPR(1 ) 

ZlOQ=  B2F I*C l 

ZlH=  -rt(l)**2ACl*TEMPiUl) 

7 l P = UB*Zl« 

M lOMs  dZF 1*0 l 
M l UP*  OH*M 1 U« 

MlQ  = X-Bit  )**2*F(n**2*TEMPR(X)-B2FI)*Ct 
M 100s  -B2F2l*Ul 

mIh  = lBll)**2*F(l)*TtMPRU)tB2IJACl 
m l P = BP**!* 

K IDN=  BB*Z|0* 
kIOP=  BB*Z 1 DP 
kIQ  = BH*Zi« 

KlOQs  BB*Z 1 DU 
Kl«  = BH»Z|* 

KlP  = tiH*Z  1 P 
Y20w=  ” b 2 1 * C 2 
Y2DP=  BB*Y20w 

y20  = Cbl  I ) **2*F  t l) )*C2*TEMPR(1) 

y2DG=  B2F  I*f.2 

y2w  = “bit  )**2*C2*TEMPW(l) 

y2P  = bB* Y2« 

m2D*  = “<32F  I *C2 

w2DP=  BF*M2Dm 

m20  = (Bll)**2*F(n**2*TEMPR(l)  + U2i-n*C2 
M2D0  = b2F?j*C2 

m2*  = l-dl 1 )**2*F ( I) * K VPP( 1 )-U2I)*C2 
m2P  = BBAfjTi*. 
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r, 


o 


0 


( 


c 


( 


k2D*=  b2CSzI*C2-HGT*Y2DW 
k20P=  HH*K20w-HGT*Y?0P 

k2Q  = (-b(l>*A2*F(n*CSZU)*TEMPKU}-BFD0FI)*C2-HGT*Y2Q 
k2DU=  -bFCSZUc2-HGT*Y2DQ 

k2"  = l9U)**3*CSZCn*IEMPH(i)  + B2D0FI)*C2-HGT*Y2w 
K2P  = b8*K2w-HGf*Y2P 
Y3DV  = -02l*C3 

y3«  = (-0(/)**2*F(m*C3*TEMPA(n 

Y3DR=  -D2F j*C3 

Y3DP=  U2ESZl*t-3 

Y>V  = -0(1)**2*C3*TEMPA<1) 

Y3P  = (D(1)**2*CSZ(1)*TEMPA(1)+D20DFI)*C3  C 

m3DV=  -02Fj*C3 

m3R  = t-0(|  )**2*FCU**2*TEMPAU)-02FI)*C3 
m3DH=  -02F;>I*to 
w>DP=  0FCSzl*L3 

m3V  = (-U(i)**2*F(n*TEBPACn-02I)*C:S 

Ki'P  S (D(l)**2*F(l)*CSZ(l)*TEMPA(l)+D2CSZI+DFDUFI)*C3 

kZOV*  L)2CSZI*C3-Y3DV*HGT 

k3R  = (0(l)**2*F(l)*CaZ(n*TEHPA(l)+UFD0FI)*C3-Y3R*HGT 
k30R=  OFtSzI*C3”V30R*HGT 
KiOP=  -0Csz2t*C3-Y30P*HGr 

k3V  = (D(l)**2*CSZ(l)»TEMPA(l)t02D0FI)AC3-Y3V*HG,J 

K’P  = C-0( i )**2*CSZC1 ) **2*fEMPAU)-2.*DCZ0FI)*t3-Y3PAHGT 

74DV  = -02l*C4 

74R  = (-0(1  )**2*FCin*C4*TEMpA(l) 

740R  = -02F I *C4 

zaOP=  D2CSZI*L<I 

74V  = -0Ul**2*C4*TEMPA(l) 

34P  = (01  t )**2*CSZm*TEMPA(l  ) +D2DDF  X)  *CU 
m4DVs  02F I * C 4 

m4M  = (D(t)**2*F(l)**2*TEMPA(l)+D2FI)*C4 
m4UK  = 02F2i*C4 
m40P  = -0FcsZl*C'l 

M4  V = (0(t)**2*  F ( 1 ) * I t'MPA  ( l ) +021  ) *C4 

M4P  = (-OCl)**2ACSZ(l)*F(n*TChPAU)-l)2CSZ(-OK)Ori)*C4 

k40V=  bB*Z/|DV 

k4H  = bb*Z/  |H 

K 4 OR-  nb*Z/|DR 

K40P=  BH*Z/|OP 

K4V  = 0 H * Z 4 V 

K 4 P = llb*/4P 

OYP  - (Y2P+Y3PJ/L3 

0 YO  = Y20/L3 

0YR  = Y3K/L3 

OYV  = Y3V/L2 

HYW  = Y 2w/L2 

OYOP  = Cy2UP+Y30P)/L4 

0 YOG  = Y2CG/L4 

OYDR  = Y 3 0 R / L '• 

0 Yl) V r Y30Y/L3 

DYOw  = Y20R/L3 

OZP  = (Z|PtZ4P)/L3 

riZO  = zto/t3 

nZR  = Z4R/L3 

r»ZV  = z a v / L 2 

oZw  = Z|W/L2 

OZOP  = ( Z 1 0P+Z40P) /L4 

rZOO  = Z 1 DQ/L'I 

OZOR  = Z4DR/L4 


Y 


v >* 
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02DV  = Z/IUV/L3 

OZOW  = ZIDH/L3 

OKP  = UiP+K^P+K3P+K«P)/La 

OKQ  = U|3tK2'3)/'l4 

pKW  S lKjH*K(|W)/l./| 

OK  V = lK3VtK(|V)/L3 

OKW  = lKiwtKCw)/L3 

OKOP  = U|OP  + k2DP  + K3DP  + K4DP)/L5 

OKOQ  = U )0U+k2DQ)/L5 

OKOR  = (KiDR  + *«CiR)/L5 

OKDV  s (K30V  + M0V)/L4 

OKOK  = (k10*+K2DW)/L4 

0MP  = CMlP  + M(|P)/L« 

0*0  = MIQ/L4 

OMR  = M/IR/La 

pMV  = MrtV/13 

0MW  = «tW/L3 

pMDP  = (MlOP  + HilOP)/L5 

nMDU  = MlDQ/LS 

pMDR  = M4DK/L5 

0*OV  = M/|MV/L4 

rtMOW  s MlU»/U'l 

ONP  = (M2P*N3P)/L/| 

ONQ  = N2Q/L4 

OMR  = N3R/U4 

pNV  = N3V/L3 

ONW  = N2W/L3 

qNUP  = (N20PtN30P)/L5 

nuOO  = N2DU/L5 

l)NDR  = N3UH/L5 

nNOV  = NiOV/L'l 

oNOW  = N20H/L4 

OY«=OV(Jsl)Zv  = OZP  = DZR=OKK  = OK(J  = OMV  = OMPai)MR  = DNwrONO=0  , 

OYOR  = DYDU=DZL)VsOZDP  = OZPH  = UKr)W  = OKOT  = l)MDV  = i;f'DP=OMOK  = DNO^=ONOQ  = 0. 

OYV=2.*DYV 

0YP=2.*0YP 

l)YR=2.*0YR 

pZ*=2.*0Zrf 

0ZU=2.*0ZQ 

pKVs2.*0i\v 

OKP  = 2.*l)KP 
f)KR  = 2.  *f)KR 
0M*sC.*O*R 
0kqs2.*umq 
nNV=2.*0Nv 
0NP=2.*DNP 

ONR  = C.*l)NR 
OYDV  = 2.*l>YDV 
oY0P=2.*UY0P 
i)YDW  = 2.*l)YnR 
0Zf)Ks2.*l)Z()l» 

CiZD0=2.*UZ00 

OKDV  = 2.*IJKpV 
OKOP=2,*l)K0P 
OK  0W  = 2 . *I/K  OR 

nMOW=2.*IJy.n" 

OMO(j=2.*opoIJ 

0NDVs2.*l)fcoV 

ONUP=2.*liNOP 

0NDRs2.*0N0R 
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S'JUROUI  INE  GLOC*A,Ort,CC,DO,AL»NST,THff KFO,FO) 

OlMENSidH  Bl2b),CSZ(2b),r(25),G(25),S(J5),i)F(2b),DF2C25),0FJ(25), 
*r>CSZl2S) .DcSZ2C25).OCSZ3(25),UCSZF(2b) , 0CSZF2 ( 25) , UCSZ2F (2b) 
*,1)125) 

COMMON  /AHC/  DRAFT ( 25 ) 

COMMON  /LOe/  naiSF(?->),CNTRCi:(23),CHINLC23),HSPRAY(23) 

COMMON  /liiloMM/NSn  (25)  » *1  (25,2b) » "2  (2b,  2b)  ,01  C2b,25) 

COMMON  MS«(25)»OEL1»OFU2#MI,N2 

tOMMOn  /o/  GU25),5tC2‘v)*Slf2rj),PHO(25),TDRAF(25) 

COMMON/ X/  TSLLTC25),OI(25)>OFT(2b),i)F2I(25),OF3I(2,j).OtI(25), 
*OC2X(25)fOc5l(2b)»OCFl(25)»OCF2I(25)»DC2FiC2b)»03BIl25)»XS«(2S) 

STATCMUM  F0NC1I0N  FOR  THAPK/OIDAL  INTEGRATION 

TRAP(H,Y1,y2)=o.G*h*(Y 1+Y2) 


STATCMt'NI  FUNCTION  FOR  LINEAR  IN TCRPOLA T ION 

STATC(X,X2,Xl, Y2,YI)=Y1+(X-X1 ) * ( Y2- Y I ) / ( X2-X U 
L CONSTANTS 

c 

wL2=wL*wL 

wL3=*L2*RL 

wLR=RL3»hl 

WL5=WL'I*«L 

rS I ART  = 8B 

CSTAHT=CC 

nTR=3.i«i*3q27/iao. 
phIZ=2.*0T» 
r)ELTA=l.*OTM 
00  999  K=l,5 
rF3=1)STaR  ( 

CC=CST  ART 

PMl=PMi/-(K-l)*OCLTA 
PHI)  ( N ) =PnI 

CCHECK=CSTARr-BO*PHI 

PMlO=CSTAWT/Bb 

C CALCULATE  QRAFT  ANO  CC 

L 

P C N I )=XS*(Nl  ) *WL-AA 

rtCNl  )=CC-HS*<UI1)  -THT*F  (Nl)tBB*PHI 
r>0  R m=  i , mj  r 
cc=cs  r ar  r 
F(M)5XS«lM)*«L-AA 

0(M)=CC“MSwlM)  -THr*F(M)«hR*PHl 

XFIOID.LE.O.U)  0(N)=-(MT*  X s * ( M ) * *L  + PS  w ( 1 ) -MSN  ( M ) ♦ 0 0 * P H I 

TF10C1).LC.0,)CL  = -AA*  fMT  + MS*C 1 ) +H0*PHI 

TFIO(NI) .Lt .0. )0(M) =TMT*(XSP(N1)-XSW(M) ) a»l+MSW(N1 )-HSw(M) ♦bb*pmi 
!FlO(M).Le.U.O)  CC=(XSH(Nt)AWL-AA)*rHT  +MH*  Ph  I 
TFIOU).GI.O.O.Ani).D(NI)  .GT.0.0)  CC=CC*00*PHI 

(F  toil  ) .bT.0.0.  ANU.Ci(Nl)  ,G1  . 0 , 0 . ANO  . TH  T . L 0 . 0 . 0 . ANU  . CCHfcCK  .LE  . 0 . 0 ) 
*CC  = PHIU-*»jP  + CSTAW1  ♦(PHI-PHJU)m2.*BH 
if ID(K).LC.O.O)  D(M)su.O 
tFIk.eu.5)  dmaf r (M)=o(m) 

IFIA.EU.3)  CCUsCC 

L 

l.  CALCULATE  GlMPfR  AND  CROSS  SECTIONAL  ahf.a 
t 

I=M 


r 


I t 

I i 


\ 


i 


0RFT=0(M) 

rSZU)=o.n 

R(I)=0.0 

sm=o.a 

TEHPl=U. 

dFu)=u. 

OF2(1)=0. 
nF  3U)=0. 

OCSZCUsU. 

0CSZ2Ci)=i>. 

0CSZ3U)  = n. 
nCSZFcn=u, 
ocszrati)=u. 

nCSZ2FU)s.j. 
jJ=NSrtU) 
kL1  = 0 

C)P  b J=2,JJ 
b02  = L>1 CI» J) 

Roi=oi  (if  j-n 

rwi  i=«i  (i,  j-n 
r«22=w2 ( l , j) 
p*21  = v<2(  I , j-l) 
jFIDRFf  .LC.O.U)  CO  TO  9 
fFlDKFI  .C£.U1U,J))  Gil  TO  7 
p*t2=SfA IE(UHFT  / RD? , liV  1 , Rrf  1 2 » Hrt  1 1) 
rh22  = SI A lfc(0KFT  , R02 , PD  1 , Kw22» H *2 1) 
kL1  = 1 
pn2=UHF  T 
7 pELU  = Rl)2"Roi 
oEL02=OEL0aPEL0 
Wll)=Hrtl.?-WWl  l 
w20=Rw22-H*21 

pELS=0.5*0ELO*(Rwl2+Rwll t«w22+Rw2l)  . 

nELG  J =SUK1  (n  tP**i0+fVLL02) 

nELG2  = buRT  (/.2P*w2D  + CCLD2) 

oELG  = DLLGUOU_G2 

rCI)=R*12+R*22 

n J n J = W » 1 1 ♦ R w 2 1 

TL)2  = DU)-Rr>2 

Sf'UM=(TC2  + u.,j<DriD)*0JMl*DELD+(  1 02  cOEU  P/3  , ) * 0 .5*DEL0*  ( w 1 0 + W2D ) 
rEMPl  = (CMHl+SM(lrt 

stnssuj+prLs 

r,ci)=GiiJ+oruG 

cszti)=TLMpi/scn 

tFlKLl.EU.D  GO  TO  12 
3 cnNTXNUE 

l2  TFIK.NE.3)  GU  TO  9 

r>FUJ=0(I)»F  U) 

0F2(  11=01-  ( J ) *F ( 1 ) 
pF5(X)=0F2( [)»r (I) 
pAHMsCSZl I ) +00-0(1 ) 
nCS/(I)sU(T JAFARM 
0CSZ2(I)=DtSZU)APAWM 
i)CSZ3U)=D(.SZ2(l)APAHii 
0CSiZFU)=O(I  J*r(I)*PAHM 

orszr2(i)=or.szr(i)AFtn 

o CSZ2F ( I J=pC 3ZF  Cl )*PANM 
9 CONTINUE 
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10 

n 


9<*9 


(N2,DLL2#S(N1) ,ql) 
(N2,DEL2»G(Nl),Q2) 
(Nl,DELl»S(l),Wt) 
(Nl » DEL 1 » G ( 1 ) » R2) 


INTEGRATES  FUR  WCTTED  SURFACE  AHCA  AND  DISPLACEMENT 

0 1=0.0 
r2  = 0 . 0 
nS  ri=NST-l 

IFIN1.NC.NSH)  GO  TO  10 
N l =N  1 + 1 
G<1  TU  11 
CALL  SIMPSD 
CALL  SIMPSQ 
CALL  S1MPSO 
CALL  SIMPSO 
Si  (K)  = ((U+pl)/WL3 
C.IU)  = (02  + r2)/WL2 
SIU)=S(1)/nL2 
TDRAFIK)=D(1)/WL 
jFU.NL. 3)  GU  TU  999 
pIU)=U.‘ 
pFIf 1 )=0. 
pF2nn=G. 

OF3IC1)=0. 

0CIU)=0. 

0C2i(l)=(i. 

0C3i ( 1 )=C. 
ijCh  1(11=0. 
pCF  21 ( 1 ) =0. 
pC2F I ( 1 ) =0  , 

p’B I ( 1 ) =G . 

pO=LSZ (kFO) fUO-O(KFO) 
pfl  I I=2»nsT 
H = XS«(I)-XS^U-1) 

H=H*WL 

a i = thap(m,0(  n ,Du-n ) 

42=TKAP(M,DFiI),0F(I-l)) 

4:=TNAP(Nf nf  2( II , OF? (1-1 ) ) 

4UsTKAP(M,oF3( I) ,DF3(I-1 ) ) 

45  = THAP(H/B(l),H(J-m 
AG=TRA('(h,nLSZ2(I).DCSZ2(I-l)) 

47=TRAP(M,DCSZ(I) .DCSZ(I-l)) 

ABsTWAPthjDCSZr (I) .DCSZF (1-1 ) ) 

49=TWAP(M,0LSZr2( I) ,DCSZF2(I-1 )) 
410=(RAI>(H,DCSZ3(1),OCSZ3(I-1)) 

A 11= TRAPlH,UCSZ2F( i) ,0CSZ2F<  J-l) ) 

OlU)=DT(I-l  J+Al/Rl  2 
ilFI  (1  )=DH  ( l-l)+A2/HL3 
pF21  ( I ) =NF 2 1 1 1-  1 ) ► A3/'-«L‘t 
nF5 i ( 1 l=OF  jI  U-|  ) f a<I/«L5 
pC  1 ( 1 ) =0L  I(1-11>A7/kiL3 
nC2l  (I)=DC2  IU-1)  FA6/-II.9 
0C31  (1 1SUC3UI-1 ) FA10/«L5 
OCF  I ( I J=UCF 1(1-1) 1Att/-Ld 
0CF2i(i)=0cF2£(I-U  + A9/'»L5 
0C2l:I(l)=Dc2t  I ( 1-1  )+A  U/WL5 
P3H1  (I)=t»3Bt  U-l  ) tHB**  3*A5/wL5 
CONTINUE 
CONTINUE 
CC=CCU 
PETUHN 
pN  D 
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SUUROUI IN£  F IN(SX,SY,SK,3N,CR,CT,S, OMEGA) 

COMMON  /All/  *R » Ctt AR » C0$0»  N£ , S l NO, U2 

COMMON  /b/  P»<.,R,X,YY,Z,U»V,M,PMI,THf,PSI 

COMMON  /F  C(jC  F / FYNCL,F INYV,FfNYR,F  I NK V , F I NKR , F X NN V , F I NNR 

COMMON  /KImVUR/  AjBBP  *UtLl»TcBAR,XFN,DOP 

COMMON  /IN/  AA, AIX,A U, AM , HO , Cb , DUMMY , D t H , QXDU , F 0, G , NST , NV AL, 
*Pl,RHU,SH,UU,i-.L»XLG,XFG,C0LL,CDNN,|--KOUDE,CC»0U,ANU,ALU0,CLD 
COMMON  /LIfT/  L TA<3fl),CLlF  f , GAMMA , XL AM 

COMMON  / I YLC/  X ARM, 2 ARM, B ACE, Y ARM ( U) » DEL JE T ( 4 ) , UUCP C 4 ) , NJE T 
1|2=U*  A u 

iKlSX.Nt.rt.O)  GO  TO  5 
nE=1X 

CLIFT=U.2*pl 
SlN(isSXNlOMrCA) 
c0SU=CUS IOmI GA) 
tc«ah=o. 1 
CRAR=(CR+CtJ/2 
4=char*s 

aR=S**2/A 

XLAM=CT/tR 

GAMMA=ATAN(0,75«CR*(1 .-XLAMJ/S) 

BHP=aB-SASlNU/2. 

r)OP  = DO  + S*CnSD/C,+OACC 

xFN=-(XLG-LHAR/(2.*rtL)  ) 

0tL=S/(NL-[) 

OELI=l./lNe-l) 
ftau)=o.p 
no  a 1=2, ne 
A pTAd  )=£IA(l-l)+OEL 
VORX=0.0 
yORYsO.O 
vnwK=o.o 
VOWNSO.O 

CALL  FlNCi>F(CR,CT,S, OMEGA) 

5 TFUHl.nt.y.b)  GO  TO  10 
rE  f A = « ( V t XpG*R ) /U 

CALL  VUR FEx(VURX,Vr)RY,VORK,VORN,rtETA  , C R , C T , S , OMt GA ) 

[0  fBE TA=-( V+xFnar) *COSU/U 
pNZIJ*CbAK/4NUASP 

CF=0.0A«/(rNA*U.1666) 

CO  = 0.1..’!3*Pr*TLHAM**2 

l-|RAf;=(CD  + 2.*CF.|.(FYNCL*FBETA)**lV(PI*AR) ) AA/wL**2*Ui> 

FlNX=-2.*0«*C 

.<?X=FXNX  + VOrX 

rY=(FINYV*v+FINYH*R) au+vory 
«5K  = < F INK  V*v  + F INKRAR  ) *U  + VOFtK 
5 N s ( F I NN  V * y AF INNRAR)  UI  + VORN 
RETURN 
END 
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SUBROUTINE  F INCUF(CR»CT .SfUMEGA) 

COMMON  /ALL/  AH,CRAR,CUSO»Ne,SIMO,ua 
CUMMDN  /b/  P,Q,R,X,YY,Z,U,V,*,(JHI,THT,FSI 

ftlMMUN  /f-CnrF/  FYNCL.F  INYV,F  J N Y R , F I NK  V , F I NK  R » F i NNV  , F I NNR 
COMMON  /FIMVUR/  AfBflP.DCLl# I CHAU , XKN, OOP 

COMMON  /In/  A A, XIX, AW, AM,8B,CB,CF,0rR,0X0U,FU,G,NST,NVAL* 
*PI,RHU,SP,ijO,wL,XLG,XFG,CULL,CONN,FRUUOt,CC,UO,AMU,ALOD,CLD 
j VOR=0 

CALL  LIFTC(0.,CL,CL,CM,CT,S,QMCGA, IVUh) 

fyncl=ll 

fbc=coso/uo 

FlNLV=CL*A/*,L**2*UO**2*F0C 

FINLRsRNLVAXFn 

SFV=FlNLV*tOSO 

SFW  = T'INLR*COSO 

VFV=F1NLV*SIN0 

vFH=F 1NL«*SINU 

fINYV=-2.*$FV 

fINYR=-2.*sFR 

FINKV  = -2.*vF  V*R0P/RL  + 2.*SFV*DDP/wL 

FlNKhi  = -2.*vl;M*BUI,/ML  + 2,*SFM*00P/V’iL 

FlNNV=F INYV*XFn 

fINNK=FINYR*XFN 

RETURN 


SUBROUTINE  LIF1C(BETA,CL,CLR,CK,CT,S,0MEGA,1V0R) 

DIMENSION  cLC(30)  ,CLCI<(30) 

CHMMUN  /ALL/  AR, CHAW, COSO, NE,SINO,U2 
COMMON  /d/  P,U,K,X,YY,Z,U,V,[)UM,PHI, rhI,PSI 
COMMON  /(•  InVUI.V  A,UBP,DCLI , 7CBAR,  XFN,  Dl)P 
COMMON  /LIFT/  LIAC30), CLIFT, GAMMA, XLAM 

COMMON  /IN/  AA,AIX,AIZ»AM,BB,CB,CF,I>TK,DXDU»FU»G,NST,NVA(.» 
*pI,RHU,SR,UO,aL,XLG,XFG,CDLL,CONN,FRiJUOE,CC,OD,ANU,ALUO,CLO 
w=0.o 

wP  = 0.0 
AlPHA=l.li 
ALPMAR=1 
oO  40  L=1,NC 

tfuvok.lo.oi  go  to  20 

XFlBCTA.EG.0,0)  GO  TU  19 
CALCULATE  s IUEW ALL  PARAMETERS 
SlNTsSIN(THT) 

OT=CCtAA*SlN( 

nF=CC-CWL-AAjASINT 

0=0F 

XFlOf.GC.O.O)  GO  TO  10 
0=-nl*SINT 
!>T  = 0.0 
10  o2=0**2 

CALCULATL  LIFT  ON  SIDEWALL 
XLIFT=CLIFT*U2*02*BETA 

CALCULATE  VORTCX  STRENGTH  and  POSITION 

5JNH  = SINUTAN(BfTA)) 

H=0.25*P1*0 

GHK  = XLIF  T/(U*M) 

Yl=SiNM*«L 

Y2=E TAlL)*StNO 

y=yi+y2 

YP=Y1-Y3 

Cl=ETA(L>*tOSO 

c2=H-Cl-'H 

c;=h+cuut 
ri=sqrt (ca*>a+YA*a> 

plP=SuRTlC2**2fYP**2) 

ni=GKA/(a.*Pi*Ri) 

OlP=UHK/l2.*PI*RlP) 

TFCY.EU.O.yJ  XMU1=PI*0.5 
TFIY.EU.G.o)  GO  to  22 
yMUJ=ATAN(AHS(C2/Y)) 

2c  wI=uj*sin(xmui-omi;ga) 

tFUP.EO.O.OJ  XMl'l  = Pl*0.5 
tFIYP.LG.0.0 J GU  TO  23 
xmuj=aTanuHS(c2/YP)  J 

23  w IP  = 01P*9  IN  ( XMIII  -UMFGA  ) 

sIuerash  calculation  for  image  vortlx 

p2»SGRT (C3a»2ayaaP; 
b2P  = S'JM  lC  j*A2tYP»*?) 


' 


! 


I 

* 


02  = 

o2P  = GRa/C2.*PI*h;:P) 

IFIY.EU.O.-J)  XMU2=PI*0,5 
tFIY.eo.g.*j)  go  to  n 
yPU2  = ATANUbSCC3/Y)) 

11  w2=N2*S1N(xmU2+UMEGA) 
rFlYP.LQ.O.OJ  XMUPSPTXO.S 
IFIYI'.LU.O.O)  GO  TO  12 
y-U2  = AlAN(/,t)b(c3/YP)) 

12  w2P=Q2P*SInUmu2*UMEGA) 

w=w) +w2 

wP=w Jp+w2p 
aLPHA=-wP/u 
*LPMAR=-«/U 
t 

L CALCULATE  FORCE  ON  FINS 

L 

20  cETA=CH-CH*F  TA(L)*<l.-XLAM)/$ 

CLOR  = AR/ (2.+AH) *r.*PI*ALPHAH 
CLU  = AH/C2.tAR)*n.*PT*ALPHA 

CLCR  (LJ=y.5*(CCrA/CBAR-f((«./PIJ*30RT(l-(ETA(L)/SJ**2)) 

* -U.-ETAIL)/S)*(‘I.*(1.-C0S(GAMHA))  ))*CLOH 

CLC  CL  J=<j.'j*CLCTA/CbAH  + ((«./PI  J aSUHTU-(ETACL)/S)**2)) 

* -C1.-CTA(L)/S)*(0.*C1 .-COS(GAMMA) J ) )*CLO 

40  CONTINUE 

CALL  SIMPSOCNF, DELI, CLC  ,CL) 

CALL  SINPSO(NC,DCLI,CLCR  ,CLR) 

rF  TURN 

fnd 

vy 


c 


■ 

l 

L* 
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subroutine  vurtfxtsx.sy.sk.sn.rcta.cr.ct.s.omega) 
common  /all/  ar,cbar,coso»ne, siNO.ua 

COMMON  /HHJVUR/  A.hMR.UELI.TCUAK.XFN.DDR 

COMMON  /IN/  AA.A1X.AU, AM, BB.CU.CF, DIR, DXDU.KO.G.NST.NVAL# 
ApI.RhO.SR.uO.WL, XLG.XFG.COLL.CUNN.KRUUUF.LC.UU, ANU.ALUO.CLL) 
1 vOR=l 

CALL  L IF  IC(BF:Ta,CL,CLK,CH»CT,S,UMEGA,IVOH) 

FlNLR  = CLK*A/',L*Aa*U? 

pINL=CL*A/wL**2*ua 

0RAGR=LLRA*2/CpI*AW)*A/WL»*aAUa 

D«AG=CL**a/ (RI*AR)*A/wL**2*ua 

SFK=r INLK*cOSO 

SF=FINL*C(ISU 

SX=-CUHAG+dWAGR) 

VFH=riNLH*SlNO 
VF  = ”FI NL  *S I NO 
SY=ST +5FR 

SN  = SY*XFN+(DKAr,R-ORAG)  *BBR/WL 
SK=CVFH-VF)*BbR/WL-SY»DDP/RL 

RETURN 

end 
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SUBROUTINE  DLRIVC(T,N,Y,YP) 

REAL  kc,mc,nc 
nIMENSIUN  A(0,5),B(5) 

0 I M£NS I UN  AlU, j),f)l(i) 
pIMENSlON  Y(l2),YP(t2) 

(COMMUN  /A/  PUOl ,QDOT» ROOT, PHJ DOT, TH TOOT , PS IDUT , UDOT , VDOT , ROOT, 

*yl)UT  » Y DO 1 > ZUU T 

C(1MM0N  /Al>C/  0RAFTC25) , WEIGHT, DUMMIE (5) , THETA 
COMMON  /O/  P,Q,R,X,YY,Z,U,V,W,PHl, THI ,PSI 

COMMON  /NOD/  PYP,OYQ,l)YH,DYV,l)YW,DYOP,DYDQ,DYl)R,DYOV,DYDW, 

* PZP,DZO*l>ZR,DZV,OZW#DZDP,DZUG,OZDI<,OZDV,DZDh, 

* nKP,DKQ,DKR,l>KV,OKW,DKOP,OKDU,OKDR,CWI)V,DKDW, 

* PMP,DMr),OMR,Of  v,umm,dmop,omou#omdu#omov,umdw, 

* PNP,DNQ,DNR,Or.V,l)Nw,ONOP,ONDG,ONDK,UNOV,ONOw 

< CC1MM0N  /OERV/  nk, DELTA, FX,FY,FK,FN,XUUELU,0SAGY,0RAGN, 

*oELTAY,OELTAN,oHACK,dETAR,OFTH,UCLTAX,OELTAK.OFTHl, T HR ATE 
*,DELP/OELS,KPM, IFOIL, IF  IN 
COMMON  /FOyL/  Cf ALFA, GAMA, XF 

COMMON  /in/  AA,ATX,AIZ,AM,HEt,CB,CF,OTR,OXDU,FO,G,NST,NVAL, 

*pl , RHO, SP , UU, WL , YLG , XFG , COLL  » CONN , KROUOt , CC , OU, ANU, ALUD* CLD 
*,NO,NG,SPTtiRN,  IPLOT ,IPT,AIY 
COMMON  /PReS/  CD I S , RNOW A , PH In , PHI 1 , A TM , PM AX , AC , DEM 
CflMMUN  / I EmP/  SX,SY,SK,SN 

common  /tempi/  thigh, tlow,shtpdg, totloc, tx 

COMMON/ Y/  DNPHI,t)KPHl,OYPHI,oCLN,OELK,OELY 
COMMON  /INOFR/  CR,CT,S, OMEGA 

COMMUN  / 1 V c C / XAHrt,Z ARM, RACE,  YARM  ( «)  , DEL  JC  T ( <l)  , RMCP  ( «)  ,NJE  T 
COMMON  / Vs  A v/  pwF.T, RAMP, vsPEH,CL'L, CAY,  IBUG,F(25), BETA 
C0MMUN  /Vs GT/  BU0YAH,TNWGT,WM0,WX0 
COMMON  /VLDUT/  voloot 
COMMON  /VOl.M/  VOLP 

COMMON  / 1 UVR/  PRM,PwZ,PMC,PZC,PSLZ,PSLM 

common  /pseal/  thtb,thts 

COMMON  /Hnw  PC,  OF,  QO,  VOOTP,  AOP  , A I 
II  = YC1) 

V=Y(2) 

W=Y (3) 

P=YCR) 
q=YC5) 
r = Y CO) 
y=YC7) 
yY=Y(B) 

Z=Y(P) 

. •.  pHI  = Y (10) 

THTsY  ( n ) 
pSI=Y (12) 

BFTARs-A  Tan( IV-XLC*R)/U) 

CALL  OHAO(ORAGY,I)RAGK,ORAGN,P,R,  V) 
tFCT.EO.li.) 

*C ALL  AUXILY(Pht »U»XUDELU,ONPhIF,DKPHIF) 

TFCIFIN.NE.O)  CALL  FTn(SX,SY,3K,SN,CK,CT,S, OMEGA) 

CALL  TMI U$T (U, VMT,DFTH,TX, TY, T K , Tn , SH I PUG , TOTLPG ) 
4R«ALL=2.*crZPL 

ECUET  =0.M 

C0I=2.*U>.5*UYV*V*U)»*2/(P1*aPrALL*ELUEF) 

L 

c CALCULATE  ijPUT  , V OUT , WPtJT,  HOOT , OOUT  , WUUT 


a(I,d=am-oyuv 
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A(1*2)=-UYDK 
Af 1*3)=-0Y0H 
A(l,«)=-i>Voa 
A(l,S)=-IJYDK 

a (2* l ) =— uznv 
a(2,2)=am-0zow 
a(2,3)=-OZoR 
A(2»4>=-UZdG 

A(2/S)=-UZDH 

a(3,1)*-okdv 

A ( 3 r 2 ) S"0K  Drt 

a(3,3)=AIX-0KDP 

A(3»«)=-I)KDQ 

A(3,5)=-iJK0K 

AC4,l)=-0MnV 

AC«,2)=-0Mn* 

A (4»3)=-0MpP 
AC4»4)=AXY-OMno 
A(4*S)=-»M()R 
ACS, 1 )=-Oh0V 

a(5»2)=-&Ndw 

a(5,3)=-0n0P 

A(b»4)=-I)ND0 
a (5*  51 =A  X Z-ONPr 

CALL  LlNYt.LCFXLV,FYLV,FZLV,FKLV#rMLV,FNLV) 
CALL  INERTxA(FXlC»FYIC,FZXC,FKIC,FMIi;,FNIC) 
sz=o. 

SM  = 0. 

SET  THRUST  FQUAL  TO  SHIP  DRAG  AT  T=0. 
IFIT.EU.L)  TXO=3HIPDG-SX 
jX=TXO 
TZ  = 0. 

tm=0. 

n«AGZ=0. 

0PAGM=0. 

,->XPHiF=0. 

nYPHIF=0. 

OZPH1F=0. 

r>MPHIF=0. 

0PAGX  = XUD(:t  U-SHIPOG 

CALL  SLAhAVlWX,WY,RZ»AK,HM,WN,VOL,AO,Y,T) 
CAUL  DRAUV ( 0Z2  , IJMP  , 0K2,  F,  W,  G,  2 ) 

CALL  DRAGV(UZ3,DM3,DK3,F,H,G,3) 
nZ  = DZ2»l'Z3 
CiM=OH2  + DM3 

llK=0K2  + DK3 
VHLPSVUL 

CALL  l>HFSS(T, Y,VOL,AU,XC»YC, ZC,KC,mC,NC) 
p*M=wm*0EM*»L/2240. 

PWZ  = HZ*1)Um/22'4o. 

PMC=mC*I;Lm*wL/2240, 
pZC  = ZC*DUm/22'M). 

sow  seal  forces  and  moments  with  roll  moment 

HEV=Y (9)*Wl 

mt  = i i 

tfsth=amxm[ tun  A/DTR+.Dooi»  lao.) 

TFttfL  rA.£Q.0..CJK.  TLSTtJ.EO.180.)  NI  = t 
0ELSL=2.*HH/NI 
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X$L=XFG**L 

YSL=-(BBtO,5*ncLSL) 

ZSL=DU-CC 

ROwzsauwKsoOrtMso. 

STNZ=S I NM=sTNK=0. 
pCGAGE=PC-ATM 
jFIPCGAUL.lt. 0.)  CO  TO  3 
; 00  1 Ii:l>Nl 

YSL=YSL+OELSL 

pALL  S«AVE(XSC,TSL,ZSL,Y,T,(:tASL) 
oBOR=ETASL 

IFIOBUW.LE.O.)  GO  TO  1 

wBU-sDBfJw/Sl'NCrHIB) 

r>FLZ  = -OLLSLA'<eoiN*PCGAGL*COS(THTB) 

R0KZ=rtU«4forLZ 

D *RMX=XSL 

ARMYsYSL 

rOHM=BO«M-01 LZaARMX 
rOWK=BUWK+oCLZ*ARMY 

1 c<JNTiNUE 
L 

C STERN  SEAL  mncFS  ANO  MOMENTS  RITrt  ROLL  MOMENT 

ySL*-XLG*ml 
ySL=-(BB+0.5*0cLSL) 
oo  2 ix=i,ni 

YSL=YSL+OELSL 

CALL  S«AVE(XSL,YSL,ZSL,Y,T,ETASL} 
nSTN=ElASL 

iFlDSTN.Lt.O.)  go  to  2 
MSTN=OSTN/SIN(TMTS) 

HYORCJsO.SiARHOAGAOSTN 
i)ELZ=-OELSL*MSTN*PCGAGE*COSfTHTS} 
sINZ=SlNZ+nELZ 

4rmx=xsl  , 

jRMYsYSL 

STNmsSInm-oELZaAHMX 
<;TNK  = S INK  tot  LZ*ARMY 

2 rYT'M  T x •-AOE 

3 cHNTiNUE 
pSLZ=BO*ZfSTNZ 
pSLZ=PSLZ/2240. 
pSLM=bO«M*s I NM 
pSLM=PSLM/22«0. 

SLZ=(ljO*ZtStNZ)/OCM 
SLK=(oOUA+sTNk)/OEM/wL 
$LMa ( BURM+sTHM ) /OCM/hL 

wZ=RZ+SLZ 
wKSRKtSLA 
wMswm+SLM 

XKGT=rLIUMT*SXN(  THT)/UL'M 
r =*L ION r*f OS (IHT) /DEM 

pX=-AM*  (U*w-H*V)  ♦ XLV  + SX+t-  X I C *TX +0R AGX *0 X PH 1 F -CO I + R X-X rtGl  + XC 
fY  = -AM-»(K*|i-m*wT  YFYLY»SYtf-YIC»TY  + DRAGY  + OYPHlF+*>YYYC 
p7s-AMA(P*V«o*lj)»FZLV  + S/  + r-/Ic»tZH)>vAUZYDZPHlF+AZ  + ZRGT  + ZCtOZ 
FX=-tAXZ-AlYJAQ*R»FKLvtSXfFMCt  I A ♦OR  AUK  + 0*  Ph  I F + «K  + KC  + DK 
pMS"C AJX-AlZT  »RAP»FMLY*SM*FMICt|M  + ()RAGMtOMPHiF  + RM  + MCtDM 
FN  = -lAiY-ATX)*P*Q  + FNLVFSNtFNJC  + TN>DRAr,N  + ONPMIF+i«N*NC 
q ( 1 ) sFY 
q(2)=FZ 


) 


>1 
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p(i)=FK 

B(4)=Kh 

i RCS)=FN 

n£U  = 5 

CALL  CUM0(a»NL-Q,o,B,1»MEK.DET) 
uDOT=FX/Um-OXDU) 

> vOOTaBtU 

W0UT=B(2) 
pOOT=8(3) 
nOOT=HU) 
rOOT  =B  (5) 
yPU  )=UDUT 
YPl2)=VnUT 

, . yP<.3)=wDUT 

yP14)=PDUT 

yPV5)=Q0UT 

ypu»j=«dut 

c 

L CALCULATE  xdot,ydot,zoot 
c 

cnSTH=cos(rHT) 

SINTH=SIN(THT) 

cnsPHi=cus(PMU 

SlNPHI=STM(PHI) 

c0SPSI=CLIS(PSI) 

SlNPSI=SXN(PSI) 

XOUTS  U*CnsTHACUSPSI+V*(SINTH*SINPMI*COSPSX-CUSPHl'*SINPSI)+ 
*W*lSlNTH*COSPHHCU8PSItSINPHI*SINPSI) 

YOUTs  U*C(igTH*slMPSI+V*CSlNTH*SINPHl*SINPS t +COSPM I *COSPS I ) + 
*W*lSiNlH*COSPHI*SINPSl-SiiNPHI*CiJSPSI  ) 
2DOT=-U*5>InTH+v*CGSTH*aiMPHItw*CUSTH*CUSPHI 
YP17)=XDUT 
yP18)=Y0UT 
YPl9)=Z0Uf 

l>  « 

c calculate  phiodt, thtdot.psioot 
c 

Aiu,n  = i, 
aU1»2)=«. 

A l U / 3)=-ST*TH 
a 1 n=c. 

A U2/2J=COSPHI 
A 1 12, jJ=CUsTh*3lNPHI 

• ••  AU3/i)  = t/. 

A 113/2) =“S  r^PBI 

A 1 13/3)=C0STH*C0SPHI 

RlU)=P 

R l 12) =U 
rU3)=H 
nEQ=5 

CALL  CUMB(aI/NlO,0,BT,  1 , t<ER  / (JET  ) 
yPUO)=U1  ( l ) 

YPU1)=U1(2) 
yP  U 2 ) =U  1 Ci) 
iFlKP.Ltl.ft)  *P  = 0 
RETURN 
£ND 
vv 


O 
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SUBROUTINE  DKAG(DY,[)K,DN,P,R,  V) 

C(!MMUN  /In/  AA,  MX,AU,  AM,t(H,CB,CI-»0  TH , DXDU, FU, G , NST, NVAL, 
*PI,RHO,SF,UU,  *L,  XLG,XF  G , C DLL , CONN 
cOMmON/x/  iSL(.H20),DI(2*5),0FIl25),0F2J(25),DF3I(25),DCIi25), 
*DC2iC25)»DLM<2!>)»DrFU25),DCF2I(25),Cu;2FI<2S),f)3L'K2S),XSK(2S) 
COMMON  III  AK, aKL,ARL2,AHL3,AMF,AR|-2, AKF 3 , APF L , AHFL2 , AKF2L , B3B 
p2=P* P 
r2=R*W 
v2=v*v 
rP2=R*P*2. 

Vp2=V*P*2. 

V»2=V*H*2 , 

VtlsV-FU*P/wL 
nNt=SIUNU.0,VO) 
tFlR.EU.U.)UU  TO  7 
X0=-VU/S ' 

7 CONTINUE 

4RCA=t)I  (NST) 

ARE.AL=OFI(mSI) 

4REAL2=OF2j(NST) 
aRL AL3S0F  3l INSl  ) 
aREAFsOCICnST) 
aRLAF2=uC2i(NST) 
aRF.AF3=OC3i  INST) 
aREAFL=OCFI (NST) 
aFL2=DCF2i (NST) 

AF2L=0C2F I (NST) 

0Y  = -CI>LL*(V2*ARFAFR2*AWEAL2fp2*AREAF2  + VH2*ARL:AL-RP2*AREAFL- 

*VP2*ARLAF ) 

r)N=-CULL*(v2*AREAL+R2*ARf AL3+P2*AF2Lf VW2*ARfc■AL2-KP^*AFL2- 

*VP2*ARLA^L) 

nK  = C0LL*(V2*AI?L  AFfK2‘*aFL2FP2aARLAF3FVR2*AREAFL~RP2aaF2L“ 
*VP2*AHLAF2) 

r>KVs-CUNN*u3llI(NST)*P*ABSCP) 
lFlK.EU.O.)Gu  TO  2 
TFIX()  + XLG)  2,2,1 
t tFUD-XFG)  3,2,2 
3 C*LL  GLUM ( x 0 ) 

AY=-CnLL*(v2*AR+R2*ARL2+P2*ARF2+VR2*ARL-HP2*ARFI.-VP2*ARF ) 
AN=-CULL*IV2XAPL+U2*ARL3+P2AAUF2LFVR2XAHL2-RP2*AKFL2-VP2*ARFL) 
AK=C0LL*lV2*Aj<r  + H2*ARFL2  + P2*ARF  3 + VR2aauFL-RF2*AR|-2L-VP2*AKF  2) 
nNtlJ=SIGN(  | . 0 , - X 0 ) 
llYs(!lY-A»*2.J  AlINCP 
• oK  = TDK-ak*2.)*oncP 

pN=(0N-AN*2.)  *l)NCP 
2 0YsnY*uNE*2. 

OK=2.*UK*OnF  +OKV 
ijN  = ON*OF'L*2. 
return 
end 

vy 

f 


i 
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SUBROUTINE  GEL'mCXO) 

COMMON  /IN/  AA,  AIX,AIZ,  AM,Ub,Cn,CI-/OfH,DXDU,Fll,G,NST,NVAL, 
ApI,RHiJ*SH,i|U*wLf  X LG, XFG, COLL, CONN 
rOi-iMUN/X/  ISIC.H25)  ,01 125)  ,DF  I (25)  ,OF2K2rj)  ,0F3I(2S>  ,DCI  (25), 
*OC2I(25)«DciU25),DrrH25),OCFCnCS).0i:2ri(25)*liiHI(^S),XSW(JSI 
COMMON  ///  AN,AWl,ARL2, ARL3,ARF,ARF2, ARf  3, AHF L , ARFL2, ARF2L/H3B 
YaU«t»Xt#X3,Yl,Y2)=Ytt(X0-Xl)*(Y2-Yt)/(X2-Xl) 

XO  = XI)  + AA/«l 

do  1 I=2»NS> 

K = I 

iFlXO.Ge.XswU-D.ANO.XO.LT.XSHCI))  GO  TO  2 

1 CONTINUE 

2 kI=K-1 

Xl=XS*lKl ) 

X2=XS*(K) 

AK  = YOCXO  ,Xl»X2,DI(K1)  ,DHK)) 
aRL=YOCXG,x1»X2,OFI(KI),DF1(K)) 

4WL2=YU(X|1,  x l,X-’/DF?I  IM  ) ,DF2l  (KJ ) 

ARL3  = YU(Xd,x  l , x2,OF3I (KJ)  »UF3l («)  ) 
aHF=Y0CXG,xI*X2,DCI(K1),UCI(k)I 
4KF2=YU(X0,X  t , X2,0C2IIK1 ) ,0C2l  (K)  ) 
4RF3=YUlXn,xl,x2,DC3KKn,t>C3l(KJ) 
aRFL  = YU(XO,X  1,X2,0CFIIK1  J.OCFKKJ) 
AKFL2=YlKx.J,Al,X2,DCF.M(l'.n,DCF2IOO) 

ARF2L  = YOlX.j,Xl  ,X2,0C2F  1(  K I)  , OC2F I (K  ) ) 
e;B=YOlXt,xI/X2,B38I(Kl),B3B/(K)) 

RETURN 

end 


' ' 


I 


» 
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SUBHUUTINt  AUXILY(PHl,U,X>JDlzLU'ONPHlF,OKPhlF) 

COMMON  /»/  P,Q,H,X,YY,Z,UUM1,V,-*,0UH,THT,HS1 

COMMON  /if./  AA,  M X,  A U,  AM,  Hb, tb, GUMMY, GTR,OXOU,f(j, G» NS T.NVAl, 
*Pl,Whl),SH,uO,WL»XLG,XFC,,CULL,CONN,l:KUUDU,CC,DO,ANU 
COMMON  /U/  GU25),  SK25),  31(25),  PHO  (25  ) , T0RAK25) 

COMMON  /*ALU/  VOLO » OR  AGO  , D£L  l)RG 
Y0lXO,Xl,X2,t l,Y2)=Yt+(X0-Xl)MY2-Yl)/(X2-Xl) 
COlDTRA)=2./lSP/SQRT(G*DJRA) )**2 
ny  al=5 

jFlU.NL.l)  GO  TO  10 

kO=NVAL/2+i 

ut’=l. 

rN=UO*SP*HL/ANU 
ARG=(AL0Glii(HN)-2.)**J 
CFU=iI.07‘j/aMG  + 0.000<I 
rN=U*SP*wL/AMJ 
4RG=CALOGlg(KN)-2.)**2 
CF  = 0.075/ARG  + .OGO'I 

S*A0=G1 1*0) 

SPAO=Si (KO) 

\/OLO=S  1 ( AO) 

TORAFU=TI)RaF  IK0)*RL 
C8U=I).0 

jFlSBAU.Le.O.O)  GO  TO  9 
CFB=GFu+G» AO/SB AU 
CPU=D.02R/hURT(CF3) 
cFR  = C'UTI)RaFU) 
jFICFR.LI . t BO  1 CBU=CFW 
9 0»AGU=(CF|)*S»A0+CB0*SBAUJ*U0a*2 
10  oO  1 1=2 , NvAL 

K = 1 

iFlPHI.Gt.PHU(I).AND.PMl.Lf.PHOCI-l))  GO  TO  2 

1 cONTXA-Ot 

2 SKAR  = Y(jCPHI,PHf)U)  ,PMU1K-1)  ,GI  CKO  ,GI(K-1)  ) 

SflAR=YU(PHf ,PHu(K) ,PH0(K-1) ,S1 CKJ,3l{K-t}> 

VOLW  = YOIPht«PHOCK)  ,PmUik-1)  ,SIUJ,SI(X-1)) 
TnKAFMsYU(phl,PHfi(K)  ,MH0(K-1  ) , TOKAK(K)  , TORAF  ( K- 1 ) ) *Wi. 

CHRsD.O 

IF ISBAK.LC.O.O)  GO  TO  11 
rFB=CF  *GwaH/SRAH 
CRR  = I).U29/SGK  f (CFB) 

CFH=CQ( TOWAFH) 
fFICFH.L! .tBR)  CBR=CFR 
pHlMa-PHl 
|1  no  3 I =2  f N V AL 
K = I 

TFlPMlM.Ge.PMU(I) .ANO.PHIM.LT.PHU(I-l) ) GO  TO  A 

3 cO NUNUL 

<1  sWAL  = YO(PHIM,PhO(K),PM(UK-1),GI  (K)  ,GI(K-1)) 

$6At.  = vU(M|trft,PHi)(K),PM0CK-n,S!  <K),S1  (K-D) 

VOLL  = YU(PhtM,PmO(K),PH0(K-1) ,SI (K) ,SI(K-1  )) 

TORAFL  = YU(PfiIM,PH(J(K  ) ,PHQ(K-1  ) , TURAF  (K  ) , TDR  A F ( K - 1 ) ) * WL 

C B U = 0 . 0 

TF ISBAL.LG.O.O)  GO  TO  12 
CFH  = CF  *t)».AL/SHAL 
CBL=0.029/SGRT(CFH) 

cfm=c(uin»afl) 

TFICFR.U  .cHU  CBL=rFR 
l2  cONTlNUL 

ll2  = l)A0 
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pRAGH=(CF  *SwAR+CbR*SHAR)  *U2 

0RAGl=(CF  *SrtAL  + CHL*S  1AU  *U2 

r)WAGV  = URAGR  + l)RAGL 

x(iR  = ORAGU-oRAGR 

XUL=0« AGU-oHAGL 

XUU£LU=XUR+XUL 

0FU0KG=-XU0FLU 

nPN  = db/'«L 

l)MPHiF  = -(XllR-XUL)*BrtN 

nKPhlF=-  IVrjLK-vULL)  *UBN*2./FR0UDE**2 

rEH/RN 

ENO 

w 


SUBROUTINE  THrtuSl (U, TUT, DPI H, TX, TY, Tk, r N , SH 1 HDG , TUTL DO ) 
OlKpNSlUf*  DFLJ  C'D  ,OP(<l)  , I JLT  (4) 

COMMON  /IN/  AA,AIX,AU,AM,8f',CH,CF.DrR,l>X0IJ,FU,G,NST,NVAL, 
*Pl , RHU,SH,uU»  HL>  XLG,XFG,CDLL,CUNN,FRUUUt ,CC,  DD,  AMJ  , ALOD , CLD 
*,wC,nG#SHTu«N 

commun  /iHpsr/  i coni , i con^, rci)N3, id»ag,ccu, tmto 

fflMMUN  / f V (.  L / XAHM,ZA«M,0ACE,  YARM(a)  , DLL  JO  ( U ) , HMC P ( « ) , N JE T 
*,ALHHA(d) 

OFL T (XX J=XX 

0FlH(  YT  )=i».o  l 33'UYY**2t0.2u67*ABS(YY) 

FMIlMSb)  = lo.<J  ' (SS/1 . 669 )**2- (9U  .*  (SS/  l . 009 )+ 5280 0 0 , 

OGHOM («wj=ww/ 1.089* ( 5900,-350. *( U.- TtON<l) ) 

»hlP(SS)=  2.8a(SS/1 .6O9)**2-10,*(SS/1.6H9)*82000. 

TOONS'. ' 

TCtlN«  = RMCH( 1 ) +RMCPC?) + RMCP ( 3} tRMCPCH) 
jFlU.NL.l.)  UU  TO  1 

CALL  KLSULnCSP.CCO, THIU,SM1PDG, TOTLDG) 

THMEAN=TUTLDG 

IFICC.NC.GCO)  CALL  RLSI)LD(5P,CC, IMT, Dummy, OUMMT) 

V = SP 

thmip=fmxp(v) 

C(IPP=.5*Kh(v*«L**2*SP**2 
cONSTl  = tOO<jOU  ,/COFF 
gUNS  f 2=0 liOyO  . / CUFF 

thmipu=tmmip/coff 

THMAWG=TtONj/LnFF 

THCON  T = T HM  I HO- VMM ARC 

THMCP=  f HClJNT 

THNEVS  = Ki1IP(Y)/C0FF 

THTUHN  = TMMi;AN 

jF  ( SPTUKN  . SP)  CO  TU  5 

CALL  i<LSULpCSPTURN,CC/THT,  DUMMY,  THTUNN) 

thmii’o=fmip(spturn)/coff 

THCONI=IMmxHu-1MmARG 

3 YFlTHTURN.gT.THCONT)  TH fURN=THCUNT  • 

0 I FF  = T HCUN T-  I HI  IJ«N 
THK9D=Thmcp  -diff 
]■> o o x=t,N,jn 

nELJCI  J=DEL T l DLL JET ( I ) ) *DTR 
nPU)=L>CLH(OI:LJCr(I)) 

|F  (AbS(DLL.jr  I (I)).  CD.  90.)  OP  ( II  =0  . 

TFIDLLJCI(I).CQ.i80.)  DP(I)=0. 
j c()N1iNUF 
r.O  TU  5 
1 Cf)N  T l NUC 

v=u*sp 

CALL  WLSUI n C V ,cc, THT, SHI  HOG, TOTLDG) 
fHMIP=FMlP(V) 

T M M 1 P U - T H M } H/LOFF 
THCUN  T = TMMjHO- 1 HMARC 
TM.XIJl)=  ThLOmT-OIFF 
THKfcY5=RMlP(V)/C0FF 
5 

rF(oi ■tm.«e.oj  go  in  20 

oo  25  1=1 ,nJLT 
»NJE  r=NJLT 

T JF.  T ( 1 ) = I Hf  liNT  / AnJCT 

TJET  (I  ) = l Jc  Tl  X)-(l  ,-HMCPd  n x CONST  1 

|F  IDLLJL  I ( J ) . fJ . 1 8 0 , ) T JF  i ( I)  = THHtVS-(  1 ,-PMCPU)  )*cOnSI2 
fFlAbSCDLL.jr I ( I ))  .L  ).90  .)  TJU  r ( I)  = THHLVb-(  1 .-WMCP(I) ) *CUNST2 


TJf  TCI)=TJETU)*U.-DFCI)/100.) 

2$  CONTINUE 
r.o  TO  10 
i6  cnf,TiNUE 
pAlR=NJEf/2. 

TMIGH=(1mkjPU"»0GMUM(V)/C0EF)/NJC  f 

pGO  = fhHQU/RA  IR 

|Fl  THIGH. GLr.hRo)  THIGH  = »00 

tlu»=wud-thu;m 

nJTsNJLT”! 

00  UO  X = 1 , NJ T , C 
TJET(I)=IhiGH 
(|0  TJETU*t  JsTLUw 

TFIDFTH.GT.O)  GO  to  10 
oO  50  ls>,N Jl,2 
TJET  CU  = lLort 
GO  tJET(I  + U=THIGH 
lO  TX=0. 

tt=o. 

tk=o.o 

JN  = 0 . 

r>u  jo  i = i,nJft 

Tl=TJET(l) 

OtL 1 =OEL J ( I ) 

ALFJ=ALFHA(I)*0TR 

TX=TX+I  l*CDSlALFlJ*CUS(DELiJ 

T Y*T Y ♦ T l *CnS  l Ai_F  1 ) *S  t N ( DEL  i ) 

TK  = TK  + Tl*SjMALFl)*YAHM(I) 
tn  = TN-I  l-*L0SlALFl)*C0!»iDELl)*YAHM(I) 
3 0 c'JNTINUC 

iFlDrTM.tO.O.)  TX=TX-OGMO«CV)/CUFF 

XKSTK- rY*lAkM 

TN=Tn-TY*<akm 

OF  turn 


r 


SUBROUTINE  Rt  SOLDI V,DHF  T, TRIM, SHI PDG.TOTLDG) 

COMMON  VQL'J,0RAG0,DEI  OHO 

COMMON  /IK/  AA,  4 IX,  A t / , AH, OH, Lb, CF  , 0 fH , DXOU, FO, G , NSI  ,NVAL, 
*PI,RH0,SH,||I  ,WL,XLG,XFG,CPLL,LOUl,FKUUOE,CC*DD,ANU,ALUD,CLI> 
COMMON  / ICmR/SX/ST,  SK  , SN,  HAVEni,,  AL  RIIOG,  HYUHUf  ,SPR  YOG,  SEALDG, 
ASRINGG'F iNpG 

COMMON  /Abe/  DRAFT  I2S)  , *E  I GHT  , BUBH , 6U6L  , H ALB , SL  BGh  , SL  S T RN , THETA, 
*OEMTH,SPHayL 
rO=RHO 

tfiv.nl. sp)  co  ro  10 

HOL=dUbB/BuUL 

HOL=DRFT/euBL 

(0  c=  VULU*Rll*G/WtIGHT*2.**L**3 
F = V/SURTir,*HUHL) 

CVT=wEIGHT/(R0*G)/BUBL*A3 

CALL  "AVEIh* L,hOL,C,F,  WAVED 

wAVE0G  = 0.G*ROag*BUHI.  »*jM:VT*«2a«yAVET 

call  ALRUCxL,DLPTH,bUbH,w>'.Lb,DRF  T,  v , AERODG) 

CALL  SRRAY(V,SPRAYL,SPRYDG) 

CALL  SLAL(pUbb,V,SLBO*,SLSTRN,THCrA,SCALDG) 
pREO=O.S*  RU*RL**0*SP**2 
SKINOG=2.*oRAG(|APREO 
FINOG=-Sa*pNLO 

SHlPUG=l»AvrUG4.ACR(10G  +SPRYDG+SEALDG+SK  INDG) /PHLO 

TOTLOG=SHlpOG+Df LURG+FINOG/PREO 
rETUHN 
END 
vv 
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»v 


SUbROUUne  BAVElBUL/HUL/C/F, TOTAL) 

Pl=s.l'll‘j<» 

M = 10. 

w 1=2.*" 

T(1TAL  = 0. 

wAVEDG=0. 

nIFF=l. 

r>0  1»  M = l,20 

FRSL=i. 

AM=M-1 

jFlAM.LO.O,)  EPSL  = 1. 
r.AMAsl  .-c 
pliJL=C/C9./:,.*H0L) 

ALFA  = ‘J.*tJI*AMAr**;VW 
flE  TA=2.*PI*BUL*AM/h1 

FACsU  . + 5grU  l.+ALFA**2)  )/SQrT(  1 ,+ ALFA A *2) 

SB=SUK I (y.b+.SASURI l 1 . EALFa**2 ) ) 

A K 1=0.5/F**2 

SIGBAsCOSCaK  l*St3)/SB-SIM(AK  USHJ/(AKl*Sa«*2) 

OELTAsAKI  *S0**2*HUL*2. 

A=a.*blOL/ (AM*SQRT(  W /F*»2))*C0S(BETA)*(l.-EXP(-DtLTA))*SIGMA 
» /SB**2 

fF l AM ) 5/6,5 

5 pSI=SIiS(Ul-TA)/(Pl*AM/Hl) 
p,Q  TO  7 

6 PSI=2.*0UL 

7 p=2./6UL*SQHT(AKl/  Wt)*GAMA*SIN(AKl*SB)*PSI 
wAVEDG=( A-h)**2*FAC*F»*2*EPSL+»AVE0G 
TFUUTAL.CO.U.)  go  to  b 

r>IhF=AOSl  (uAUOG-TUTAD/TVrAL) 

B tot  AL  = WA  VErjG 

fFlOIFt  .LE. 0.001)  GO  TO  99 
tO  CONTINUE 
<j9  p E T U R N 
pM) 


I 
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iJ 


- > r 

- its.i 


jUHHOU  T l t»E  AtKO(SLFT»Ol-PTH,Hf  Bl.UHK  I »V,  ACUOOli) 

ANU  = 1 .!>(>£— j4 
RUsO.OOCiB 

pENOLO=V*SLFr/ANU 

CF  = 0.45b/4|_UCluC»CNOLO)*<r2.5e 

4REA  = SLFf  Afb  + HJt  (Oi;PTh-nHF  n*2,) 

FPJM  = CDUPT^”('RF  1 )-*(b  + Hl*2.) 
pR£  = D.i)*i<0*VAA2 

FRMT06=Pt<E*0.fa*FRUNT 

SKIhDG=Pi<E*CF*ARCA 

aF«OOG=F  KNTUG  + SKINOG 

RETURN 

pNO 

vv 


SUBROUTINE  SPRAY tV,SPNAYL,SPRYDG) 

C OMMJN  /A6c/  f>R*F  TT251 , WEIGHT,  huba,8Ut>L,w  ALB,SLfiOw,SL  ST  KN,  THt  T A, 
*pEPTH 

COMMON  /CDE/  nRISC(2iJ,ENTRCli(J3),CHINU'(l*'4),M5PHAY(23) 

COMMON  /IN/  AA, AIX,A  U , AM , «(< , CP , CF  , D f R , DXOU,  HJ,  G,  NS  T , N V AL , 
*pI,RHU,SP,UO,w|_,XLG,XFG,CULI. , GUNN , FROUUE , CC » GO , ANU 
CdMMOIM  / S F. s / hS«C25) ,0LL1»&EL2,NI,n2 
RU  = RHiJ 

fAC=5. 10150/180. 
pO  10  1=1, NS f 
to  mSPRayu;=0. 
pO  30  1=1, nST 

ANG  = SINCI)RlSL(I)*FAC)*SIN(ENTRCCCn*FAC) 
mSPRAY ( U=v*<2/(2.*G)*ANG**2 
CHA=CHlNt( j)-PRAFT(I) 

TFICHK.U .y.O)  CHK=0.0 
IF IHSFKAY ( j) .GT.CHK)  HSPRAY ( 1 ) =CHK 
3O  cONTiNOC 

CALL  SlMPSdCNl ,PCL1 ,H8HPAY( I ) , ARE A 1 ) 

CALL  SiMPSMCN2,DCL2,HSPRAY(Ni ),AHEA2) 

area=ah£awakfa2 

m=N1+1 

U=V*C0SlLNTKCrcM)*FAC/2.) 
r RF."iOLL)=V*SPWAYL/ANU 
CF =0.0/5/ (aLOG10(RLNOLD)-2,)**2  t 0,0000 
pHt=0.5*H0*U**2 
SPNY0G=PHE*AWfA*2.*CF 

RETURN 

fNO 


SUBROUTINE  SLALCB»V,SLhOR»SLSTRN,  lME  f A,SC.ALOG) 

COMMON  /AHc/  PRAFTC2S) 

COMMON  /IN/  AA, AIX,A U , AM  , OH  , C 0 , / / , D I W , f>XDU  , F 0,  G,  NST  , N V Al.  , 
*Pl,WHU,St',UL,»*L»XLr.,XFG,COLL,CI)rJN,l--R0UOE,CC/OO,  ANU 
COMMON  /Sts/  HS«(r>5)  ,L)IL1*DF.L2»NI,N2 
►I  j = N 1 + I 
PO=KhO 

pRE  = 0.ti*Kij*V*»2 
4NG=f«ETA*".  I 4 159/180. 
rHO*=0. 

Stl=URAF  Un3)/SIN(ANG) 

TFlSLl -Gt.sl.UOw)  SLl=SLBOW 
rOwSL=SL1*(,USIanG) 
tFlKU«SL-tt.U.)  GO  TO  10 
aENOlO=»U«Sl.  *V/ANU 
CF  = 0.0'l‘l/ CrTNOEO**!  1 ./6.)  ) 
rBUKsPRCA’UBORSL  *CF 

to  CONTINUE 
rSTWN=U. 

SL2=0RAFI(|  J/SINCANG) 

TFISL2-GL.SLSTRN)  SL2=SLSTKN 
sTKNSL=SL2aLUS(ANG) 

TFISIMNSL.tC.O.J  CO  TU  20 
r9£N0LD=S I RmSL  AV/ANU 
-cF=0.0<WI/(RrN0L0x*Cl./6.)) 
rSTKN  =PKE*a*STRNSU  *cf 
20  sEALOG=R«OwtHSTRN 
rETUKN 

end 
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gllbKOUf  1NE  UNVFL(FXIV,FYLV,FZLV,FKLV,FPLV,FNLV) 

COMMUN  /d/  P,Q,R,X,Y»Z,U,V,'/(,PHI,ThI,PSI 

COHMUN  /NOl)/  nYP,DYQ,UYR,OYV,OY«,OYbP,OYUQ,DYDr<,OYt)V,DYDW, 

* OZP,DZQ»UZW,PZVf  O/H,HZ')P,L)ZL)U,nZDH,0ZOY,0ZDW, 

* ''KP,DXQ,OKH,DKV,DKW,m<OP,DKDQ,DKOI<,DKDV,DKOW, 

* nMP,OM'J,lJPR,CMV  ,DMrf,  OMOP,OMOlJ,DMDHf  DbDV.nMDh, 

* nNP,ONGl,ONR,OHV,DN«,DNPP,L)NDQ,UNUK,DNOV,ONOW 
FXCVrO. 

FYLV=(OYV*VHJYwAhtOYP*P+DYfi*Q+OYR*R)*U 

FZLV=(OZV*vtUZw*w+OZrj»P+OZHAQtUZH*H)*U 

FKLV  = (UKV*V  + l>KW*H  + OKP*P  + OKQ*Q  + l)KR*K)*U 

EMLV  = (0*v*v  + UM/i*W  + DMP*P  + 0NIJ*OYOMRaR)iMJ 

FNLV  = CDMV*vtPf<W*w  + ONP»P  + l)MIJ*(l  + ONR*R)*U 

bETUWN 

£NO 


1 


; 


SUBROUTINE  IMRTlA(FXXC,FTIC,FZIC,FHIC,FNIC»FKiC) 

COMMON  /o/  P»G,R,X,YY,Z,U»V,w,PMI,THT#PSI 

COMMON  /NOn/  nyP,OYf),OYR,DYV,OY.'l,DYi)H,OYDO,t)YOI<,OYOV,DYO«/ 

* PZP,DZG,OZH,OZV,PZw,nZUP#OZO(,,OZURf nZDV,DZ0R, 

* PkP,OKQ.i)KR,OKV,(iK*,OkOP,OkDU,DKOH,DKOV,OKDh, 

* PMP,OMU,UMH,OMV,DMW,OMi>P,l'-MDQ,OMUR,OMOV,OMDH, 

* nNP,UN(J,ONR,DNV,r.NW,nNUP,UNUr,,ONUI<,|JNUV,(JNO« 

frX IC=  -OYOV*«*V-OYUP*R*P-UYOR*H*KtOYDw*n*a+DZUO*U*U+DZOP*P*Q 

fyic=  -r>znw*«*p-ozoo*p*o-oznp^p*p 
fZIC=  oyov* v*p+uYDR*p*R+uYOP*r*p 

pKIC=  -DTOR*P»VtOYOP*R*V+(DKOP-DNOH)*P*fi+ONDP*(H*R-P*P) 

* tOZDP*M*R+UMPP*Q^R 

FNIc=-UZUP*w*Q+UZOQ*«*Pf(DMOO-l)KOP) *P*Q+DKOQ* (P*P-Q*Q) 

* -OYDP-»V*(i-ONnP*Q*H 

FKIC=CUZUw-OYnv)*V*w-(DYDRfr)MDR)*R*«+(OZOQ+ONOV)*V*Q 
*+(0NUH-OMOQ) *R*U-0YDP*P*rt+UZ0P*Y*P-l)MDP*R*P+0N0P*0*P 
RETURN 

end 

y* 


i 


: 

; 


! 1 
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SUBROUTINE  St AwAV(wX,*Y,WZ,WK,HM,WN,VOL,AO,Y,  n 
SEAL  MwK,My(M,MwN,MTM,MTN 
p I BENS  I UN  Y ( 12) 

COMMON  /IN/  AA,AIX,AIZ,AM,BH,CU,CF,OtR,DXDU#FO,G,NST,NVAL, 

*P I » WHO  »SR»uU#/>L»YLG,X('G»CDLL,CUNN,FKUUUE,CC,IjD»ANU*ALI)0,CLO 
a.nc.ng.sriuhn, I PLOT , IRT , AIY 
COMMON  /SCS/  HSW(?‘j)rUiai.0EL2»Nl,N2 
COMMON  /(Ran/  F rx#MTM»MTN 

COMMON  /"FnR/  r/»X(25),FwY(25),F*Z(25),MWK(25),MWM(25),MwN(25), 
*AREAC25)>F(.FAK(2'j) 

COMMON  / * G T / BUOYAN, INWGT, RMO, WXJ 
COMMON  /USlCAK/  BLEAK, SLEAK 
po  I I=t»NsT 
I CALL  BUOY(j,Y,T) 

jnTEGRATL  for  WAVE  FORCES  and  MOMENTS 
0 l =32  = 05  = 0/1  = 05=06  = 0 . 
mSTi=nST-1 
nT  = n1 

iFINl  .NE.NS (I)  GO  TO  5 

nT=NST 

C,0  TO  6 

5 CALL  SIMRSO(N2,OCL2,FWX(NT),OI) 

CALL  SlMHSn(N2,DCL2,FwY(NT) ,02) 

, CALL  SIMRSnCN2,i)CL2,Fw/(NT),03) 

CALL  SIMMS0(N2,DCL2»M^K(NT) ,06) 

CALL  SImwsoIN2,DCL2,m*m(nT),q5) 

CALL  SIHMSnCN2,0rL2,M/.N(NT)  ,06) 

6 CALL  SlMKSn(NT,DCLl,F*X(l),Rl) 

CALL  SIMMSOINT ,DCL1 ,FwY Cl ) ,R2) 

C ALL  SlMR6(HNT,ULLl,FWZ(l),Hi) 

CALL  SImMS(1(NT,0CL1  ,MWK  ( 1 ) ,R4) 

CALL  SIHRSnCNT,DCLl ,Mwh( 1 ) ,R5) 

CALL  SIMRSlUNT,OCLl ,Mim( 1),R6) 

CALL  SIMMSnCNJ  ,ULL1  , ARl:  A ( 1 ) ,R7) 

CALL  SIMMSijC  n(,!)CL1  ,Ft.EAK(  1)  ,38) 
pfMENso  .‘j*rMu»wL  **2*SR**2 

WX=(01+R) »FtX)/OCMCN 

wX=(OltRl  J/OCHCN 

wY=(02tR2)/DLMEN 

wZ  = (0itR5)/DLML'N 

W<  = (Utl*K**)/DEM|;N/wL 

wM= (U5+R5»mTM)/0CMCN/rL 

WM=  ( U5  + R6  ) /OCmCN/ziL 

wN=(OotRo+MTN)/OCNEN/«L 

wN=(u6tRo  )/i)CMLN/™L 

V 0L=H7 

aO=R3 

wLFAK  = KtAK(Nl)-*O.S 
SLEAK=)LtAK( l ) *G  .5 
iFUNWGT.Eg.U)  GO  TO  2 
wx=«x-«xo 
^ms«m-«MO 
i CONTINUE 
bE (URN 
pNO 
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SUBROUTINE  UUUY  ( T , Y, T ) 

i*  definition  or-  parameters 

BLAL  MnSK,M,.bM,MV<SN,MwPK,M».pM»MWPN,MWK,M*M,MWN,MBK,MBM,MTM,MrN 
DIMENSION  .!  IKAfUa)  , SGN(4 ) »OINT (4J , Y( 121 
DIMENSION  awl  AD(25),r>FT(25,4),bCAM(2'i,4) 

COUPON  /(ifl)M'V  nSm(25),«1'.25,2S), 42(25,25), 01(25,25) 

COMMON  /I\/  AA,  A IX,  AIZ,  AM,bH,Ch,CF,0  I R . DXDU » M.1,  G,  NbT  , N V AL , 
*Pl»RHU»SF\l.'W,"L»XLG,XK-,i:DLL,CDN.N,|--RUUOe,LC,UO,ANU,ALUO,CLO 
A,NC,NG,Sr'rijHIM,  IPLOT,  [PT,A  [Y 
COMMON  /3LS/  HS*(25) 

COMMON  / I w R b I / fCONl,  IC0N2, rcONi,  IOHAG,CCU,THTU 

COMMON  / I WAN/  FTX»MTM,MTN 

common  /wav/  r>WFT,WAMP,«PEH,CCL*CAY,lBUG,F(25) 
common  / WFi)H / F*X(25) ,FmY(2S) ,F«Z(2S) ,MWK<25) ,MRM(25) ,MWN(25), 
*aRLA(25)  ,F|_FAKC2a) 

COMMON  /X/  DUMMY (300) ,X$W(25) 

COMMON  /bCCM/  hEME (25) , BEAM  I (25) ,8tAMf I (25) » BE  AM|'2I (25) 

*,BLAMF 51 (2S) 

COMMON  /ciFEM-/  UCM2(25)  ,t)EM3(25),  ARMS,  ARM  P 
COMMON  /OSCAL/  ZARMSL , AREASL 

D A T A SGN/-J  .,l.»«l«»  I ./ 

TFU.EU.1  .AND.IBUG.NC.O)  WH  I TE  (6, 202)  A A , BB  , CC  , DO,  Y ( 1 ) , Y (7  ) , Y (<>) 
«,YU0),y(t|),Y(l2) 

2()2'  FORMAT  ( Ml,  aAa,BB,CC,UD,U,SURGC,HFAVC,PH1,ThT,PSI*/10GU’.4) 
FWX(X)=FNY(t)=FwZ(i)=M«K(I)=MwM(I)=MWN(I)=0. 

aHLA(I)=0. 

aREAD(X)=». 

BfMgd  )Sli, 

REM2(I)=Bf-.M5in=0. 

JJSHSWCXJ 

tFUJ.LU.l)  R|.‘l  urn 
rHOG=RMO*G 
1)0  I K = 2, 3 
jIRAN(K)=i> 

nFT  ( 1 , K)  =rlEA.-'(  I , K)=0. 
pO  2 J = 2,J.J 
0 l TOP-O t ( I , J ) 
nltfUfsOl II ,J-1 ) 

Ml  TOPswlU.JJ 
M 1 BO  fswld,  J-l) 
w2TOIJs«2  ( I , J ) 
w2BOfs»2(I,J-l) 

CALL  SLGAL(K,W1 fOP, W t DOT , R2  TOP » M2BUT  » D 1 TOP , D 1 DOT , Bt T A , H YPOT ) 
w T s t*  b 

fFlK.Gf.2)  wIs.MT 
Z = OD-MSw(r)-(JllOP 

wGT  = CC-nl  TriF-wS-Kn-f  ( 1 ) * T AN  ( Th  TU ) 

CALL  S*A/E(FIX),WT,2,Y,T,ETA) 
hCHKshG  T tfc  t A 
J F (MCHK.GF ,0.)  GO  TO  2 
nFT  ( i,B)=n  j 1 *jP MCHK 
«t  AM(I,k)=wX  IC'P 
jTRAN(K)=J 

pI-IT(A)s.JTT(t*X)-01HOT 
J F ( 0 X N 1 ( A ) # G T • 0 • ) CO  TU  t 
,jTRAN(iv)  = t 
CO  TO  1 

2 c""TlMIJt 

jTH  AN  ( K ) a J,J 

pFT(i,rO=0|  I UP 
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BEAH(I,KJ=w1 TUP 
OlNT(iO=iilTl)P-01BUT 
! C0NTI'JUC 

RE*2CI)=BEAfUI,2) 

rF«3CI)sbEaHU,3) 

1 F l OF  I ( I » <2 ) . LT . G . ) OF  T ( I ,23=0. 
jFinrUI,3).LT.&.)  OF  T ( I ,33=0. 

A«EADU)  = rt.b*Or  f (1,2)  *(Bf  AMC  / ,2)+wl  (X,  1)  )+0.5*OFT(I,  3)*(9EAP(  1,  3)  t 

I "111,1)3 

4«MS=brt  + 0.25>(PFAM(I,a)  + wi(l,  I)) 

4BMl3  = -Cbtlt.j.2^*(BCAM(  1,3)  tWl(I,  1)  )) 

nINTU)=DlNU2) 

rtiNr(<i)soiNTi3j 

pAl.L  vULUM,r(l  , ao,JTHAN,OXNT,DwCI,FLK,AK,ASEAUfHSCAU 

aPIlAII  JsAW 

fLEAK(1)=Fl»< 
pH=Y(l )*SP*"L/ANU 
tHU  = (ALUGl.j(HN)-2.)**2 
CF  = 0.075/AhG+.0G0'I 
nPU=CC-F  U)*  f AN(  THTQ) 

FwXCl)=  -C.F*(OFT(I,2)tOFTCr,3)-2.*OPO)*UHO*(SP*Y(l))**? 
pWY(l)30« 

, F"Z  ( X)=-AK|5  AO(  I ) aHHOG 

s4  M K ( X ) 3-KHpG*  (APMS*0.lj*OFT(I,2)A(HEAMCI,2)+Wl(I,  1 ) ) ♦ 

* AHMP*0.‘j*nF  H I«3)*(UEA*(l,3)ti»l  (I,  l)  ) ) 

MHMQ)=  aW|£AL)(i)*F(I3»HhOG 
mWN(I)=0, 
rRTUKN 
p NO 
vw 


SUBRUOflNE  S*AVF.(XC,  YC,ZC,Y,T,b'TA) 

DIMENSION  yU2) 

COMMON  /IN/  AA, A ix, AIZ, AM,B8,CH,CF,0rR,0XDU,F(l»G,NST,NVAL, 
*Pl,RHU,SP,ljU,»L,XLG,XFG,Ci)LL,CU''lN,FRl.RJOL,CC,ll!T,ANU,ALUO,CLll 
* , NC  » NG, SP  HlNN , I PLOT , I P T , A I Y 

COMMON  /aav/  P"FT,«amP,hPER,CLL,CAY. IBUG,F (2S1 /BETA,  IN, WOE  P , OF  E SE  T 
*, WLG, ICU,XO, HO, E I AO 
COMMON  /«GT/  BLNlYAN, [N*GT,WMO,«XU 

SINH(U)  = lCxPlU)-CXP(-iJ))/2, 

SFCHCAKG;=n./(CXP(ARG)*EXP(-ARG) ) 

nATA  CU'.l,C01»U)2,C0j,C(la/ll  .5i‘>206‘.i<>,-52.7(.71<>255,  10  7.1876202, 
A.iQO.OOBBHid, 3G.2307 187a/ 

HEAYL=Y (V)**L 
pHI=Y (101 
jHisYinj 
PSI=YC12J 
pSO  = BE  T A“P$  l 
COSI=CiJSIThT) 

Cw  = T*  aL/8P 

l|T=(Y(7)AC0SlBCTA)tY(B)*3lNC0ETA))A"L 

apgj  =xc>  i ani  imi  ) -me a ve/cos r-YC  a tan  (ph i )/cosr 

4Ri,2=(XC*CnS  r + ZCASINlThD  ) *COS  ( PSill  t YC  * S IN  ( PSO) 

' qO  TU  (1/2,31,1* 

c 

c sINUSUIOAl  wave 

c 

1 pTA  = -"AMP*«5lN(CAYAXC) 
tf  uNNi;r.Ea.o)  RcruRN 
ct=cll-*t**l/sp 

pi A=-"AMP*S IN(CAY*( AHG2+UT-CT1 l/COST-AHGl 
pETURN 
L 

l.  SOLITARY  *aYL 

L 

2 nFFSLT  = O.G'j*CrLA*PER 
A l=CAY*(AC-0FF5l'TJ 
PTA  = «AMH-*SeCH(A1  )*■»? 

IFUnaGT.EQ.O)  RCTURN 
CT=ClU*T" 

T = AdSC*Jr-CT  ) /*lg 

A t=LAY  * ( APfj’  + Ur-CT  ti.)F FSCT  + 1 * wLG) 
p T A = » AMP*SgLH ( A 1 )**?/COST“ARGI 
pETUHN 
L 

L pXPLUSIUN  wA  YE 

L 

3 fTASO. 

lFUHnGT.FQ.01  RETURN 
H=*OCP 

TO=XO/3QHT(G*H) 

|i)  = BO. 

Tw=Tn»TO 

psOTtXU 

pf-sR/  I */bQpl  U*H) 

TF(Rr .GE..31  GO  TO  0 
y=l  ./('l.*«F**21 
C.O  TO  B 
/|  qf2s«F*KF 
pFj  = HF.:*HF 
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rF«=«F3*KF 

X =CU<l*KF/|tcUJ*rtF3*C02*fiF2+COl*RF>cau 

CAYsX/H 

nMliGA=CAY*SUNT(G*TANH(X)/CAY) 

f fc  L = UMLG A/c  A Y 

CT=CLL*r>i 

4*2=2. »X 
SHk2=S1NM(hK2) 

„RG  = h*..’/Shk2 
aRG3=1  . + ARq 

j«GU  = -ARG3/CARGA  ( t ,-HK2/TANH(HK2)  )+0.5AAf(G3**2-ARr,i) 

R 0 K = C A Y * H C) 

CAUL  HLSiiCu(i»RtJK,BJ3i 

aR05>=(cay*  ( akgc+iit-ct)  j/cos  r 

pT  A = (£  f All*RU/l.’)  *SURT  ( Afi'G'D  *f)j3*COS{  AKG5J-ARGI 

rF  TURN 

fNO 


/ 


S'JBHOUI  INC  VULUMCdf  HH,  JTRAN,OINT,D*ET,FLCAk,  AREA) 
t sUbHUUriNE  TU  CALCULATE  AUtA  BC!"t.tN  CALM  WATER  SURFACE 

l.  a Ml)  taE'f  DEfjK  and  LEAKAGE  FUR  CROSS  SE  T ION  I 

DIMENSION  ,J  TRANCE)  fDlNT(<l) 

rUMMUN  /UFOMM/  NSw(?S),wi(25,25),R2(25,25),Dl(25,25) 

FLEAK=U. 

C sTAWDUARi)  s t DCW  ALL 

dS=01MT(2) 

JSSJ TRAN(2) 
jFUS.LQ.I)  GO  TO  1 
jSl=JS-l 

hGTS=0*ET-01(I,JS1)-0S 
r,0  TO  2 

1 hGTS=D«LI-dS  C 

fLEAK=-03 

l POUT  SIDERaLL 

2 op=GlNI(i) 

JP  = J IRANll  3) 
jFIJP.CG.t)  GO  TO  3 
jPl=JP-l 

hGTP  = D»E  f-oUI,  JP1  )-DP 
r,o  to  'i 

3 hGTP  = u«L  I -Dp  C 

' fleak=flEak-dp 

R »Rt A = Mt>* (HGTP+HGTS) 
return 
f nd 
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I 


SUBHOUTINC  SI  GAL (K,  rOP,wlBOT/H2TOPr*J0l)T,OTOP,D8Ul,0tTA,HYPUT) 
SUBROUTINE  TO  CALCULATE  SEGMENT  LENGTH  AND  ANGLE, 

OOIFsDIuP-dHlIT 
hDIH'  = «1  TUP_»  IBOf 

iFiK.Ld.c.nP.K.FQ.3)  Ht)ir=H2Tnp-H2i3or 

rETAsA  fAN2f rt DIF /DO IF) 

MYPLl(=S(j«T(HUlr**2  + 00IP**2) 

PC  T URN 
pND 

4 


* 


SUBWUUI  INE  OitAr,V(DZ,UM,OK,F,K,lj,K) 

fOMMON  /BEE*/  BlAM(25),6CAMI(2B),ijEAMPl(25),BEAM|'2I(25) 

* , FEAHF31 l?5) 

COMMJN  /BEct'2/  Bi;m2(25),0EI"5(2B),ARmS,ARmP 

COMMON  /IN/  AAjMX,  AIZ,AH,BB,CH,r.F,0  I H , 0M)IJ , F(j,  G , NST , N V AL  , 
*PI»Rri()»SP,u,1»'',l.»XLG»XFli, COLL, CONN 
COMMON  /A/  DUMMY (300) , X3w (25) 

CPMPIJN  ft  1/  i)li,BRL,HKL2,BRL3 
DIMENSION  f 125) 

TKAP(H,  y 1 , y2)=o-5*H*(n  + Y2) 
on  3 l=t»NsT 
RE AM ( I ) = BEm2 ( I ) 
tFIK.FU.5)  BLAM(i)=hEM3(I) 

3 CONTINUE 

, w2="*w 

g2=Q*Q 

w«2  = **')*2. 
nNL=SIGN l 1 ,0 , H) 

IFIU.EU.I..)  CO  TO  17 
yO=«/Q 
1 7 CONTINUE 

REAMlU)=0tAMFl(I)=BCAMF2I(l)=BEAMK3I(l)=O, 
nO  Ml  1=2, nSI 
u=xs«(l)-xswu-i) 
h = h*  *L 

rF  I =UE AM ( I ) *F ( i ) 

BFIi=HLAM(l-l)*F(I-l) 
r!-  2I=HT  I *F  ( I ) 

Rf2Il=bFIUF  ll-l) 
rFjI=BF21*F( l) 
pF3Il=BF2l 1 *F ( I - 1 ) 

Rl=!MAPl«,RlAM(I),OEAM(I-l)) 
q’=THAM(M,HF  l, HU  1 ) 
b3=TwAP(m,bF  2I,«F?I1 ) 
q.jrTRU'(h,nF  ^1  ,bF  111  ) 
qf  t'il  l )~0r  A.  'i  ( 1-  t )♦«  i/«L«»2 
rlamfui  )=m  Ai«r Ki-i  )m?/»i  **3 
rfamFSKD-OI  AMF2I(I-l)*B3/nL*A<l 
3EAM(  31  ( I) -HI  AMFjI  ( I — l ) + UM/w(_*-*5 
1«  fONtlNUt 

jULAsBLAMI (NST) 

»WFAL  = oCAMKl(f;of , 

jPE AL2=HE  AmF  2 1 ( Nb  T ) 

»kEAL3=mLAmF jI(nST) 

nZ=-C0NN*(*2*A«rA fupxAwr  al2-»iU2*arlal ) 
oM=  CUNM*(,,2XANr  AL  + ‘)2»AREAL3-*02*AREA12) 
fl-lO.EU.C.)  uij  10  12 
fFIXU  + Xl  U ) 12,12,11 
ll  jF(XO-XFG)  13,12,12 
1 3 fALL  GLOMV(XU) 

RZ=-CUNN*(w2*r*Wto2AHHL2-<*«2*BRL) 
q M=-C0NN»  („2'*hKLM)2*bNLi-«  02*  BkL2) 

DNLP  = S10m( | . , — x 0 ) 
o/=  ( ) jtriNCP 

Oii=  ( l)h-Mri«^)  *(>nFP 
lc  oZ=OZ*UNL 
nMSOM*L'NL 
4om=ARHS 

fFlK.L0.3)  AHMsAKMP 
OKS  OZ*AM“ 


4 

■ 
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gUbWOUT INC  GLOmv ( X08 ) 

C 

C vEHTiCAL  CRAG 
C 

COMMON  BFAM(25) ,rtEAMI (2S j ,HtAMFl (2b) ,BLAM12I (25) 

*,PEAMF5H25) 

CflMMUN  /It./  AA,MX,AU,AM,aB,CB,CF,DIR,DXOU»l-U,<;,NST,NVAL, 
*pIfRHO»SH,||U»KL,XLG,XFf.,  COIL,  CONN 
COMMON  /X/  OUMMYC300),X.'><*(3‘j) 

COMMON  !L\/  C.R,URL,BKL2,0HL3 
YOlXO,Xt,X2,Yl,Y2)=Yl+<XO-Xl)MY2-Yt)/(XC-Xl) 

XC'B  = XUb  + AA/wL 
DO  i I s 2 » N S T 
K-l 

lFUOH.G£.xSM(I-l).ANt).XOB.LT.XSH(I))  GO  TO  2 
t CONTINUE 
2 K l =K - 1 

Xl=X3w(Kl  ) 

X2=XS*(K ) 

bH  = Y0(XUO,x1 ,X2,3CAMr ( K 1 ) ,6EAMI c k > j 
gRL  = YU(XUo,X  t ,X2,8EAMFHKI)  ,BEAMFI(K)  ) 
RRL2=YUCXOR,Xl,x2,BEAvF2i(Kn ,HLAMF2I(K)) 

F(RL3  = YO(XOB,Xl,X2,UEAMFiI(Kn  , BL AMK3I ( K ) ) 

»E  TURN 
END 

VV 
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1 1 

T j 


SUBROUTINE  PKL'SSC  T,Y,VUL  , AO,XC,YC,ZC,KC,MC,NC) 

PEAL  KC  /M(,»nC 
llIMENSXUN  y C 12) 

CUM"1).-J/AI)C/  DRAFT (25)#WEIGHT 

COMMON  n M/  AA,»IX,AUf  AM,BO,CD»CF,OfR,DXDU*M),G,NST,NVAL# 

*pl  r RHCJ, SP, ijO»KL»XLG,XFG, COLL,  CONN, fHUUDL/CC* DU,  ANUfALODfCLO 
COMMON  /FR|jS/  CDIS,RHiJHA,PHiri,PHIl,A(M,PMAX»AC,OEM,lPrt 
COMMON  /"G  f/  BijUYAN,  [ N»C T , * hq , -XI) 
fOMMUN  /l-LO*/  PC,UF,QU,VDOTP, AOP,AX 
COMMON  ZbSLl  AK/  BLEAK  , SEE  AK 
HEV=Y (9) 
pH  I = Y CIO) 

THT  = YC1U 
4F=09. 

iFCT.NL.O.)  CO  TO  1 
pH  = 0UUYAN*ol.-M 
PC=(ME1GMT-FB)/AC  + ATM 
pcgage=pl-4tm 

4I=(PHI(I+PhI UPCGAGE)/CCDCS*SQRT(2.*CPCGAGE)XHHOWA) ) 

PM4X=ATM“PhIU/PMII 

g(l  TO  10 

l 4S0=£.*Ob*(bLrAK+SLEAK) 

4 T = Ai  + AU+As8 

C,AM=1.0 

pDXF=PULD-ATM 
TFCPDIF.GT.O.)  GO  TO  20 
OF=PHIU 

OL=A I*CDXS*SGRT(2.*A8S(PUIF)/RHUWA) 

(jO=AU*COIS*SQRT  C2.*A8blPOIF)  /RhU«A) 

TFCXPR.nl.  a)  MRX  IF.C6.202)  f / PULI) 

3|)2  FORMAT ( IX,*PC  l.FSS  THAN  ATMOSPhCHIC  PRESSURE*, 

*r;X,*f=*,F  1 'j  . £ , 5 X , * PC  =*  , F t 0 . 2 ) 

C,n  TO  2 

20  tFCPULC/.LT.PMAX)  GO  TO  3 

OF  =-CJlS*A|-*50R  f ( AOS(POLO-Pm4X)/RHOhA) 
fjl.=-A  r *GDTS*SiJRT  (P.aPOIF  /RhllwA) 

00B”All*CDXS*SlJR  'f  (2,*P0  IF/RHOnA  ) 

TFCXPh.nl. .j)  WRITE(6,2l)3)  T,POLD 

2o3  FORMAT (U,*PC  GHLATER  Than  PMAX *, 5X , *T=* , F 1 0 .2, 5X , *PC=*, F 1 0 .2) 
CO  TO 

3 nF  =PhIU+phI U(POLO-ATM) 

fjL  =-aT*CdIS*5uRT(?. *CPOLO-ATM)/RHURA) 
no  =-ao*cdis*5QRT (?.*(polg-atm)/rhowa) 

£ vOLOllT  = OF  + QL 

v=vnLO+vui.nuf  •*(  r-TnLO)**L/SP 
pC=PULO****GAK/VOL**GAM 
iFCPC.LT. ATM)  PC'.  = ATM 
10  pPIF=PC-ATm 

7ARM=00-CC-Hfc  V**L 
4 N = AC  * T H T 
3 T A M = w L * f A M C fMT) 
p TtiL  = H I AN-gLl  Ak 

TFCTHt.LI  .0.)  btpl  = htan*sleak 

TP19THL.GE.10.)  M0L=t8. 

TFC3T0L.LT. -18.)  0T0L=-18. 

TFC8LLAK  .GT.u.)  ZARM=Ul>-0.5AQTUL 
iFCSLLAK.GT.O.)  ZARM  = L)O  + O.b*(5T0L 
TFC  TbLLAM-SLLAK) .NL.O.)  AN  = 2 . *bB*HTBL 
XC  = AN*PDIF/jJCM 

yc=-phx*ac*pdxf/dem 
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c 


<->  o 


c 


zC  = -ACtPI>lF/Ol'M 
kC=-YC*ZARm/wL 
mC=XC*ZAKM/«L 
NiC  “ 0 • 

pOLOsPC 

TdtO=T 

VUL0=VUL 

tFU.'^.W)  BCTUHN 
vpuTP=vuLonr 
4nP=VUL 
REIUrtN 
p NO 
vv 


r,  i 

I 


I 
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r-  r. 


' 


! ] 

SUBROtn  INE  PLOT  f(Al,  A2,A3»  A*l,  Ab,  A6,NP,NC,NG) 

<?F  AL  NUMt'j) 

CO«MUM  /HRmJ/  0 f,NSTi;P#NPRNT,  IP 

pIMENSION  AUSOO) rA?(500) ,A3(5oO) ,A«(bOO),A5(500},A6(SOO) 
pIMENSlL'N  I<:LK0(o),KANGE(6),  A ( 1 32) ,3  TAR (6)  ,RSfc  T (o)  , R AI4GFI  < 6) 
nATA  IZEK0/21,G0,in,21,6b,Ul/ 

| 0aTa  R3Cl/.l,.2,l.,2.,4.,10./ 

0ATA  s I AW/ 1HA, 1HA, 1H*» IHOf 1H0»  l HO/ 
pATA  oLANK/lrl  / 
pA  T A DASH  / 1 H *•  / 
pATA  EYE  / l H I / 
p A 1 A PLUS  / l H + / 

oa r a rEE/iHtv 

pATA  .MUH/tHl),  nil,  1H2,  1H3,  1H4,  IH'J,  1M6,  IH7, 1H8,  1H9/ 

* jF  IQ)=IZ+G*SCL 

L iNIUAtUE 

c 

NG  = 6 

nSET=6 
nC  Ts 1 31 
MX  = 10 
kX=1u 

c 

t SCALING  FOR  A X I S 

no  2o  i=i,ng 

bO  TU  (21, 22> 23. CM, 25, 27), I 

2 1 CALL  HAXX(nP,A1,VMAX) 
rC'  TU  2b 

22  CALL  HaXX(nP,A2,YHAX) 
rO  TO  26 

23  CALL  NAXX(NP,A3,yMAX) 
r,n  TO  26 

2'l  CALL  MAXX(nP,A<|,YMAX) 
pO  TU  26 

25  CALL  HAXX(nP»A5,VMAX) 

C,0  TO  26 

27  c*LL  HAXX (nP, AO, YMAX ) 

26  r4LL  SCALECNSL1 ,RSET,YHAX,RNG)  ' 

pAilGL(IJ=HNG 

20  RAivGH  ( U=-(?ANG(.  (1  ) 

pPlMT  Y-AXlS 

wRITL(0,20yJ  RANCH ( 1 ) , RANGE < 1 ) , H ANGM ( 2 ) , H AnCE ( 2 ) , R ANGM ( 3 ) , H AhGfc ( 3 ) 

2oO  Fl'HHAU  iHl,Fb.  1 , I2X,7HHEAVL  = *,  1 3X  , Fa  , l , <1 X , f b . 1 , 1 1 X , 7HI’  1 TCH  = * , 

•*  i a x , r o . i / ox  » f b . i , i o x , 1 1 hx a ve  hc  r . =* , n x , r-  a . n 
rBITL (6,20 | ) UANGHIO) .RANGE (6) ,RANGM(5) .RANGE (5) , R ANGM ( 6 ) , H A NGE ( 6 ) 

2(1 1 FORMA  r ( lX.Fb.l,  J3X,6HK0LL=0,  1 3X  , , F 0 . 1 , 0 X , F b . 1 » 1 3X , 5MY  A A=0,  1 fl  X , 

*F«  . I »0A,F5. 1 . 10X, 11HPRISS/ 100=0, UX.F4.1) 

X.  pWEPA«L  PlnTIlMG  A K K A Y 

1 ^ ^ 

oO  I 

OH  2 K s l f M(,  1 


SCI.SKC/KNG 
IFU.NL. 1)  GU  TO  5 
Xl'I  = U*NC 
ri.05.u-NC 
kNT=10 

00  0 K=ILO,lHI 
aOOsHLUS 

xfiknt.ne.nxj  go  ro  6 

a(M=LYE 

K N T = 0 

6 «NT=KNT+1 

5 r,U  TO  17#8,9, lu, 1 1,  14),J 

7 a(IZ)=HLUS 
rcssrcAi t:j; 
r,0  TU  12 

8 a(IZ)=HLUS 
jC=3F  CA21 1 ) ) 

00  TO  12 

9 aUZ)=HLUS 

jc=sr-  ( a3  1 1 ) ) 

r,o  TO  12 
to  jC=SF (A41I) ) 
on  TO  12 
u xr=SFlA51I)) 
o'l  TU  12 
td  Tr.  = SF(A6U)J 
t2  aCIC7=STAP(JJ 
3 CONTINUE 

TF1KX.NE.Nx)  GO  TO  13 
tZl  = ULRU(t) 

TZ2=iZLKU(2) 

TZ3=lZLHU(3) 

a(IZ1)=uash 
aUZ2)=0ASH 
» (IZ3)=0ASH 
XFll.EO.l)  GO  TO  16 
Kl  = U-l)/t'J  + .l 
k2=k 1 + 1 

XFlKl.Gt.lo)  GO  TO  17 
A(1ZJ  + 1)=M.|NIKC) 

A(  IZ2+1)=NU'Ukc) 

a(  1 Z3-1 ) =Nuh IK2 ) 
oO  TO  16 
[7  t2  = NUD(K1 * 10) 
x 1=(K1-I2)/IU+1 
12=12+1 

a(IZJ+1)=NupU1) 
a ( IZ  l+2)=M.|NU2) 
a(IZ2+1  )=NuNU1) 
a(U2*2J*MiHU2) 

a ( 1Z3-2)snijmU  i ) 

A ( IZ3-1)  3MJN  U2) 

(b  kX  = <1 
1 3 kY  = kXU 

wR  I TL  10*202)  Ca(k),K=1»NCT) 
2n2  FOHMA I 1 IX, 1 31  At ) 

1 CONTINUE 

no  13  1=1, NCT 
to  a(I)=MUANk 
XZ1  = ULRU(  t ) 
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TZ2=iZCRU(2) 
tZ jsiZLRuC^J 
4CIZ1)=1EE 
»(IZ2)=TEE 
4CIZ3)=TEE 

wR 1TE (0 f 202 ) C ACK) *K= I »NCT  ) c 

' RETURN 

fNO 

vv 


i 


i 
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SUbROUTINE  MAXXlhP.A,YHAX) 
OlMENSiOn  All) 
yMAX=A  C 1 ) 
oO  1 I=2»NP 

1 yMAX=AMAXt(YMAXrA^S(A(I))) 
r£TIjKN 
£ND 
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r r- 


!• 

t’ 

i 


^UOROU  I I •'I  f PL01XY(XP,YYP,NP) 

DIMENSION  A ( ‘)  I ) * RANGE TO)*AX(lO) 

DIMENSION  xP13o&),YYPCS00),IP(S00) 

C{1MM0N  IP 

data  BLANK  /lH  / 

DATA  ETL/IhU 
1-iaTA  0 ASM  /lH-/ 

IjATA  PLUS  / i N > / 

DATA  STAR/] H*/ 

DATA  RANGE  /a.  ,5.  ,10./ 

data  AX/lHo»  IHI#  l M2,  IMS,  1H0,  1H5*  1H6,  IHY,  IHdt  l H9/ 

L SCALE 

L 

i nR=3 

nri=nm-i 
kXAX1S=0 

yMAX  = ABS(YYPU)J 
no  ii  i=2»np 

nPSAV=1 

|FIXP(I).LT.O.O.O«.ABSCYYP(I)).ST.BO,)  GO  TO  18 
t 1 yMAXsAMAXl ( Y M A X » A8SIYYPCI) ) ) 

1 3 mPsnPSAV 
yHI=80. 
ylo=-yhi 


flKDEH  X 
DO  I I=t»NP 

knt=i 

oCI  2 J=t*NP 
IFU.EO.J)  GO  TO  2 
tFiXPUJ  .GT.XP(J)  ) KNT=KNT+I 
2 CONTINUE 
fPUJsiXNf 
l CONTINUE 

L fNIT IAHZF 

L 

rZ=81 

jFIYYPCNPJ.LT.O.O)  12=11 

kSAVC=0 

(_INE  = 1 
NC  = t)0 
nL=P  1 

SCL=NC/YMJ 

mY=SCL 
I Y = I 0 
kY=LY 

nX=(NY*S./jO.)*2. 

I X = l 0 
KX=LX 

fF  lYYP(NP)  .GL  .0.0)  wRITE (6 *201) 

2ol  FORMAT ( I m!/jOX,7hY  VS.  x/ 

* >0X*3M  + fii;,  l 7X,  3M  + 60,  I /X,  3*<  MO.  17X,3H*20.  IPX*  l HO/  ) 
rFTYYPtNPl.Ll.O.O)  «RITf (6*202) 

2()2  FOR"AUlMt,b6X,7HY  vs.  X/ 

*30 X,  1M0,  l8X*->M-20*  1 rX,SH-Oi),  17X.3M-6I),  I 7X.TH-00/) 

L 

t PREPARE  PLOTTING  ARRAY 


c 


145 


oil  3 K = t/NL 
a(KJ=PLUS 

rFUl.Nt.LY)  CO  TO  14 
A ( K ) =E YE 
KY  = 0 

1'1  *Y=KY+1 

3 CONTINUE 

oO  9 1=1/ Np 
aIIZ)=PLUS 
qO  4 J=l/Np 

4 jFl  I0(  J)  .l-Q.  1)  IO=J 
TY=Ii-iCL*YYP(lO) 
jX=J+(XP(In)*3CL*5./10.)*2, 
fF  IKoAVE.LQ.O)  GO  TO  3 
XFUX.NE.KSAVL  ) GO  TO  7 

3 TFUX.G1.LINL)  GO  TO  7 
a(IY)=5TAK 
ksavl=linf 
lFU.EU.NPj  GO  10  7 
r.n  tu  9 

7 IFIKX.NC.LX)  GO  TO  10 
*XAXiS=KXaxIS+1 
IFIKXAXX3.eQ.1J  GO  TO  17 
IFIYYP(MP) .LT.0.0)  GO  TO  15 
rZUNG=IZ+t 

f ZT  nO-LZ  + 2 
GO  TO  16 
15  i Z0NL=1Z“2 
lZTi“0=XZ-t 

(6  a ( IZUNL)=AX  CKXaXISJ 
a(IZ1*U)=Ax(  l) 
l7  A (IZ)=UASh 
kX  = 0 

ID  kX=KX+1 

IF l A CIZ J .Nf.OASH.AND.A ( IZ J .NL.S1  ARJ  A ( IZ)=PLUS 
WHITE  1 6»  Ei'y)  UlIJ),IJ  = t,NL) 

1)0  fOHMA  r C30X  , 9 l a 1 ) 

fFU.  EG.  up.  ANP.  LINE.  cu.  IX)  CA1L  EXIT 
of)  tj  K=t,Nl_ 

8 ACKJ=6LANK 
L I NE  = L X NL+  1 
GO  TO  5 

9 CONTINUE 
PETUKN 

end 


SUBROUTINE  RUNGS  (X,H,N, Y/YPRIML, INDEX) 

D I PENSION  yU2)  /YPRtMEi  t2),Z<l2),*l  ( 12)  ,W2C12)  »«5(12)/«aua) 
LRUmGS  - WUNGC-MjflA  SOLUTION  OF  SET  UP  F‘IRSI  ORDER  U.O.L.  FORTRAN  OP 
i.  dimensions  must  be  set  for  each  program 
U y INUCREnUI NT  VARIABLE 

l M INCREMENT  DELTA  X,  MAY  BE  CHANGED  IN  VALUE 

L N NUMBER  OF  EQUATIONS 

L Y UEPENDfM  VARIABLE  BLOCK  ONE  DIMENSIONAL  ARRAY 

u YpRIM£  DERIVATIVE  3LJLK  ONE  I) I m(-.nSIlN AL  ARRAY 

L TRE  PROGRAMMER  MUST  SUPPLY  INITIAL  VALUES  OF  YU)  TO  Y ( N ) 

L jUDCX  IS  A VARIABLE  «MR»i  SHOULD  HE  SL T TO  ZERO  BEFORE  EACH 

C INITIAL  ENTRY  TO  THE  SUBROUTINE/  I.G.,  TO  SOLVE  A DIFFERENT 

t SET  UF  CUUaUOnS  OR  IU  START  RITH  NEW  INITIAL  CONDITIONS. 

L THE  programmer  MUST  WRITE  A SUBROUTINE  CALLED  DLRIVt  WHICH  CUM- 

U pUTES  I HE  DERIVATIVES  AND  STORES  THEM 

l the  argument  list  is  subroutine  derive(x,m,y,yprime) 

jF  (INDEX)  5/5,1 
1 dO  2 I=t/N 

M UI)=H*YPRlMfU) 

£ 7(l)sYll)t(M(l)*.5) 
a=x+h/2. ■ 

call  DLKIVeIA/N/Z/YPRIHE) 

D(’  3 I - 1 / N 

w2U)=M*YPRlMt'U) 

3 z(I):mjT.5*W2U) 

iVSX+H/E. 

CALL  OLRIV(:(A,n,Z/YPRIME) 

do  a I-t/M 
w3U)=h*YPr[MC(I) 
a Z(I)=YU)+W3U) 

4 = X + H 

CALL  DERIVE  IA,N,Z,YPPIME) 

1)0  7 I = l , N 

W*H  I ) sh* T PR  I Ml (I) 

7 Y(n  = Y(lJ  + (U2,*(B2(t)+«5(l)))+wi(I)Fwa(I))/i,.) 

X=X+H 

CALL  DERIVE  (X,N,Y,YPR1ME) 
r.O  TO  tj 

5 c 4CL  DERIVE  (x,N,Y,YPRImE) 

iNUEXsl 

6 pETURN 

gND 

vy 
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SUBROUTINE  CUM9(a,N,ND,0,m,NERR,0) 
t SULiiTION  uF  S]MULT.CU.  PORTING  KUlIA  CONO. 
OIMLNblU.M  A(l),  9(1) 
pUUl VALENCE  C I , F I ) , CK,TK) 
n=NEHR=l 


10 

l>n  60  1 = 1,  N 

MISSOObO 

alJMAX  s A ( I J 

MIS-80070 

tjmax  = l 

MISSOOBO 

lFlN.EU.Dr.0  Tr.'  30 

151)  20  J=2,N 

MISS0060 

|J  s I + ( t ) *NI) 

MISSOI00 

TFUBSUt  I.j))-Aas(AIJMAX  ))25, 25,20 

2° 

4 I J MAX  = A(lJ) 

MISS0120 

jJMAX  = IJ 

MISSO130 

CONTINUE 

* MISS 0140 

jF  CAIJMAX)  30,999,30 

MISSOISO 

30 

oil  00  J = l , N 

MISS 0160 

tJ  = I * (J-l)*ND 

•U3S0170 

4 ( I J ) = A(tJ)/AIJMAX 

MISS01  0 

0=0*  A[ JMAX 

MISS 0190 

00  41)  J = l,M 

MISS0200 

jj  = I ♦ (J-l)*NO 

MISS 0210 

a o 

«(IJ)  = 0(IJ)/AIJMAX 

MISS0220 

OP  70  K=1,N 

M I SS0230 

fF  (K-U  50,70,50 

MISS0240 

sO 

jm a x = Umax  ♦ ix-i) 

MISS0250 

4 F A T = -A(XJMAX) 

MISS 0260 

kJ  = K 

MISS 0270 

jj  = I 

MISS02  0 

OG  60  J=1,M 

M1SSO290 

JF  CACXJ))  55,08,55 

MISS0300 

5 5 

4 ( K J ) = AHaT*A(IJJ  t A (K-J ) 

MISS 0310 

s3 

*J  = KJ  ♦ 60 

MISS0320 

6') 

jj  = IJ  ♦ MO 

M1SS0330 

4(KJM4X)  = 9.0 

MISS0340 

KJ  = K 

•IISS0  350 

T J = I 

MlSSOinO 

0*1  09  4“t  , M 

MISS  il  370 

I F («(IJ))  65,68,65 

AISS03  0 

65 

q(K.J)  = APal*B(IJ)  + «<XJ) 

MISS0J90 

6 a 

uj  = KJ  + 60 

M l S3 0 4 00 

69 

r.l  = IJ  + 60 

MISS0410 

7 0 

(■Oi'lT  1 4UC 

MISS0420 

*J  = UMAX  - I + l 

MISS0430 

60 

4 C X J ) = F I 

M I S S 0 4 4 0 

00  100  I*1,N 

M1SS0450 

K = I 

* I SS0460 

63 

Tl  s KAMI}  _ NO  + l 

M I SSOm  70 

FK  = 4(11) 

M 1 S SO  4 0 

rF  (K-l)  93,  100,95 

MISS0490 

60 

TJ  = I 

MISS0500 

JK  = X 

MISS0510 

Oft  99  J = l,M 

m I SSOS20 

4(2)  = DIM) 

U3S0330 

n ( I J ) = a(jKj 

HSS0540 

q(lK)  = A(J) 

M ISSI'550 

|J  = IJ  ♦ 6() 

M 1 SSOSbO 

t)6 

JK  = IK  + UO 

M15S05/0 

loo 

r on  f [t<uc 

MISSUS  0 

nFRR  = -J 

M1SS0590 

SUBKIJUIINC  blMPSU  (N,H,Y,A) 

c routxnl  io  hi  RruK*  s iMpson  integha r i un  fuu  even  mjhreh  of  iNCwtntNrs, 
c 

C tgrNUMBCK  OH  INOFt'CNDC NT  VAKIaBLCS 

t H=1MCHLHLNT 

t yrlNDEHENOENT  VAHIAOLt 

L a=INTLGWAL 

L 

OlHENSlON  y(l) 

A|N  = N 

j F 1**U0  INN , 2 ) • NE • 1 ) 00  TO  10 

15  k0UNT=0 
Nl=N-l 
4 = Y ( 1 ) + Y U) 

THlHDM=H/3, 

00  1 I=2»Nl 

IF  (KOUnT  .E<?.  I ) GO  TO  2 
4=A+«.*YtI) 

KOUN'fsl 

r,o  to  l 

2 A = A + 2.*UI) 
kOUNT=0 
1 CUNT1NUC 
A=  THXHL'h*A 
PE  I0HN 
JO  a =0.0 

DO  20  1=1, N 
a = A + Y C I ) 

>0  CONTINUE 

A=A-0.5*Y( jl-0,5*YCN) 

A =A*H 

PETIIHN 

£N0 


SOHHOUIINE  Ht  S$Fl ( IU>  X » V ) 

OlMtNoIUN  T ( l 0 0 0 ) 

T«=l ./ 12. 

TH=1 ./ j. 

n«=io 

M=3.+S.*Xa*  lH-f«).AX**TH  + A*AXl  (OH»X) 
JF  ,,Nl  .0)  M = N*l 

M l =M- l 

m2=m-2 

T(R)=0 

T ( M 1 ) = 1. 

Z=2./X 

J=H2+t 

mx=m2/3 
SnCRH=0. 
lid  J l = t»Mx 
JsJ-l 

TU)=J*Z*T(J+l)-T(J  + 2) 

J=J-1 

TlJ)=J*Z*T(J+l)-TCJ+2) 
t SNUWM  = tiMjHM+l  ( J) 

SMJKM  = J.*5NUllM-T  ( 1 ) 
V=T(iU+l)/SNORM 
return 
fNO 
vv 
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from  the  cushion  air  dynamics  as  well  as  to  the  environmental  excitations  due  to  waves. 
Non-linear  contributions  including  effects  due  to  large  motions,  viscous  flows,  and 
control  logics  are  considered  and  the  ship  responses  are  solved  numerically  through  time 
domain  integrations.  Sanple  examples  have  been  exercised  using  this  analytic  model  to 
examine  the  effects  of  several  chosen  environments  on  an  SEV  and  the  results  are  pre- 
sented herein. 
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